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ABSTRACT 

• 

A  theory  for  the  analysis  of  stresses  in  laminated  circular 
cylindrical  shells  subjected  to  arbitrary  axl symmetric  mechanical 
and  thermal  loadings  has  been  developed.  This  theory,  specifically 
for  use  with  pyrolytic  graphite  type  materials,  differs  from  the 
classical  thin  shell  theory  In  that  It  Includes  the  effects  of 
transverse  shear  deformation  and  transverse  Isotropy,  as  well  as 
thermal  expansion  through  the  shell  thickness. 

Solutions  In  several  forms  are  developed  for  the  governing 
equations.  The  form  taken  by  the  solution  function  Is  governed  by 
geometric  considerations.  A  range  In  which  the  various  solution 
forms  occur  was  determined  numerically. 

As  a  sample  problem,  the  slow  cooling  of  pyrolytic  graphite 
deposited  onto  a  commercial  graphite  mandrel  was  considered.  In- 
vestlgatlon  of  normal  and  shear  stress  behavior  at  the  pyrolytic 
graphite  -  mandrel  1  iter face  showed  that  these  stresses  decrease 
In  magnitude  with  Increasing  E/Gc  ratio  and  Increasing  deposit  to 
mandrel  thickness  (ha/h^)  ratio.  This  Implies  that  a  tMn  mandrel 
and  a  material  weak  In  shear  are  desirable  to  minimize  the 
possibilities  of  flaking  and  delamination  of  the  pyrolytic  graohlte. 
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NOTATION 


X 


» 


a.  b  subscripts  indicating  upper  or  lower  lamina 
a1j*  bij  constants  defined  by  equations  (C.2)  and  (C.3) 

A,B  constants  defined  by  (C.8) 

c  preferred  direction  in  a  transversely  isotropic  material 


C1 


roots  of  equation  (22)  defined  by  (C.ll)  and  (C.13) 

h<-  (1  ■  a.b) 

1  -  v4 


120  -v 
d/dx 


e  natural  base 

E,  Ec  Young's  Modulus  In  the  plane  of  isotropy  and  "c"  direction 
respectively 

Gc  shear  modulus  relating  stress  and  strain  across  the  plane 

of  isotropy  (=G  =  G„  =  G  =  G  J 

r  \  «xz  «0Z  «zx  «z0/ 

0 

g.j  constants  defined  by  equation  (C.4) 

ha,  hb  Individual  lamina  thicknesses 

h  ha  +  hb 

l,  j  subscripts 

kj  constants  defined  by  equation  (C.6) 

m,  n  constants  defined  by  equation  (C.7) 

m  constant  defined  by  (11a) 

MX,M0  stress  couples 


5 


^Tx*  ^T0  thermal  couples 

M  defined  by  equation  ^41) 

Nv,  Na  stress  resultants 

X  0 

^Tx*  ^Te  dermal resultants  defined  by  equation  (41) 

P11  °z  ^V2)  (*  =  a»b) 

p21  °2  (-hi/2)  (1  »  a.b) 

P(X)  Pia  -  P2b 

H  h 

Pj  joint  normal  stress  (=c2  (~y-)  *  °z  (-Ky) 

Qj  shear  resultant 

^  defined  by  equation  (41) 

R  Radius  to  shell  reference  surface 

T  temperature  measured  from  the  stress-free  temperature 

of  the  material 

Ujj.  u2  deflections  In  the  "x"  and  "z"  directions  respectively 

uoj  axial  deflection  of  lamina  middle  surface 

uf  virtual  displacement  of  shell  middle  surface 

Wj  radial  displacement  of  lamina  middle  surface 

wi  “  jZ  (acT)i  dz  (i  =  a,b) 

wf  virtual  radial  displacement  of  lamina  middle  surface, 
x  axial  coordinate  for  cylindrical  shell 

z  cylindrical  coordinate 

a, a  thermal  expansion  coefficients  in  the  “x"  and  "z"  directions 

respectively 
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3 

1 


r 


c 

tj 

V 


V 

C 

vu 


0  u 

Tj 

e 

0-j  •  02 

(  )' 


rotation  of  the  normal  to  the  undeformed  lamina  middle 

surface  due  to  deformation 

defined  by  equation  (C.  10) 

defined  by  equation  (C.9) 

virtual  rotation 

defined  by  equation  (21) 

strain  component 

Poisson's  ratio  in  the  x  -  9  plane  *  vQX) 

Poisson's  ratio  (^x2  =  vQz) 

Poisson's  ratio  defined  as  the  negative  of  the  ratio  of  the 
strain  In  the  j-direction  to  the  strain  In  the  i -direction 
due  to  a  stress  in  the  i -direction 
stress  component 

joint  shear  stress  (~o^  (-hfl/2)  =  ozx  (+hb/2) 
cylindrical  coordinate  in  the  circumferential  direction 
defined  by  equation  (C.10) 


* 
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I.  INTRODUCTION 


The  ever-  expanding  missile  and  space  technology  continually 
demands  materials  capable  of  maintaining  structural  integrity  at 
very  high  temperatures.  Of  late,  attention  has  been  focused  on 
refractory  materials,  their  anisotropy  in  physical  and  mechanical 
properties  making  them  ideally  suited  for  a  wide  range  of 
Insulation  and/or  structural  applications. 

Of  the  many  refractory  materials  possible,  pyrolytic 
graphite  (PG)  has  probably  received  the  most  attention  of  late 
although  it  was  known  to  Edison  (1)*  in  1883  who  described  methods 
for  its  manufacture,  the  technique  Involving  formation  of  carbon 
deposits  onto  substrates  heated  in  carbon-con;a1ning  gases.  For 
structural  use,  pyrolytic  graphite  Is  generally  deposited  at 
temperature  from  3500°F  to  4000°F  in  a  stream  of  hydrocarbon  gas, 
such  as  methane,  onto  a  substrate  of  commercial  graphite  maintained 
at  temperatures  of  1500°F  to  5000°F.  The  rate  at  which  the  material 
Is  produced  depends  on  a  number  of  factors  which  Include  the 
temperature,  the  reaction  pressure,  the  hydrocarbon  flow  rate  and 

♦Numbers  In  parenthesis  Indicate  corresponding  references  In  the 
bibliography 
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the  surface  to  volume  ratio  of  the  substrate  surface  (2),  (3), 

(4),  (5).  X-ray  analysis  of  the  resultant  deposit  shows  a  well- 
crystallized  structure  having  much  in  common  with  the  single 
graphite  crystal  (6).  Growth  is  always  normal  to  the  substrate 
surface  and  after  a  thickness  of  0.1“  -  0.5"  is  reached,  the 
deposition  process  is  stopped  and  the  deposit  allowed  to  cool  for 
several  weeks. 

The  result  of  such  a  .formation  process  is  a  material 
highly  anisotropic  in  physical  properties.  The  PG  has  one  plane 
of  isotropy  parallel  to  the  mandrel  surface  (x,e  direction  --  see 
figure  1 1-1  for  geometry)  and  a  single  preferred  direction  (z), 
a  state  commonly  referred  to  as  transverse  isotropy.  With  thermal 
conductivity  between  100  and  1000  times  greater  In  the  (x,e) 
direction  than  in  the  (z)  direction,  the  material  acts  as  an 
excellent  conductor  along  its  surface  but  also  as  a  good  insulator 
In  the  thickness  (z)  direction.  The  coefficient  of  thermal 
expansion  in  the  (z)  direction  is  from  10  to  30  times  greater 
than  that  in  the  plane  of  Isotropy  (x,e)  so  that  thermal  expansion 
through- the- thickness  must  be  considered  in  many  analyses  of  the 
material's  thermal  behavior. 

Other  curious  effects  due  to  the  anisotropy  are  manifest 
In  the  Poisson's  ratio  which  Is  negative  in  the  plane  of  Isotropy 
(v  .  *  vft  *  -0.21)  but  large  and  positive  in  the  preferred  direction 

X0 
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(vw,  *  vfl,  *  +0.9).  Furthermore,  the  ratio  of  the  elastic  modulus 

XZ  oZ 

In  tne  isotropic  plant  to  the  shear  modulus  in  the  transverse 
plane  (E  /G  or  E _/G_ _ )  may  range  from  20  to  50,  compared  to  an 
E/G  ratio  of  2.5  for  an  Isotropic  material  with  0.25.  There¬ 
fore,  In  an  analysis  of  a  structure  composed  of  such  material,  1 

transverse  shear  deformation  even  for  thin  cross-sections  must  be 
considered. 

Thermal  and  mechanical  properties  can  be  found  readily 
(2),  (5),  (7),  (8),  (9),  (10),  (11),  (12).  Typical  properties, 
given  In  Table  1-1,  are  taken  from  (10),  which  are  reasonably  close 
to  those  given  In  other  references,  discrepancies  most  probably  being 
due  to  variations  in  the  deposition  process. 

Among  the  earliest  analyses  of  structures  of  pyrolytic 
materials  were  those  of  Garber  (13)  and  Levy  (14)  who  treated 
thermal  stresses  in  cylindrical  and  spherical  shells  and  also 
the  residual  stresses  caused  by  the  general  anisotropy  of 
pyrolytic  graphite,  but  neglected  transverse  shear  deformation 
and  did  not  account  for  the  high  thermal  expansion  coefficient 
In  the  (z)  direction.  McDonough  (15)  has  considered  thermal 
stresses  in  shells  of  revolution  of  pyrotytic  graphite  type 
materials  subjected  to  axially  symmetric  loads,  Including  trans¬ 
verse  shear  deformation  and  thermal  expansion  through  the  thick¬ 
ness.  He  was  able  to  show  that  neglect  of  transverse  shear  de¬ 
formation  would  lead  to  an  over-estimate  of  the  stiffness  coefficient 
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TABLE  1-1  PHYSICAL  PROPERTIES  OF  PYROLYTIC  GRAPHITE 
MECHANICAL  PROPERTIES 

Young's  Modulus  (PSI) 


TEMP  (x,e)  Direction  (z)  Direction 


70°F 

5.4xl06 

1.5  xlO6 

1000°F 

4.3xl06 

1 . 29x1 06 

2000°F 

3.5xl06 

1.05xl06 

3000°F 

2 . 7x1 06 

0.81xl06 

Poisson 

's  Ratio 

70°F 

"xe  =  °'21 

vx2  «  0.90 

THERMAL  PROPERTIES 

Thermal 

Expansion  In/in  -  °F 

70°F 

0.0 

13.1x1 O”6 

1000°F 

0.6xl0"6 

13. IxlO’6 

2000°F 

1.2xlO"6 

13.1x1 O”6 

3000°F 

1 . 7x1 0”6 

13.1x1 0”6 

Conductivity,  BTU/hr-ft-°F 

70°F 

290.0 

1.25 

1000°F 

165.0 

0.82 

2000°F 

100.0 

0.60 

3000° F 

60.0 

0.60 

for  properties  representative  of  PS  while  excluding  thermal 
expansion  in  the  (z)  direction  leads  not  only  to  erroneous 
stress  predictions  but  that  even  the  sign  of  the  stress  (ten¬ 
sion  or  compression)  may  be  wrong.  Kliger  (16),  (17)  has  ex¬ 
tended  McDonough's  work  for  the  case  of  conical  shells  in  that 
he  derives  equations  for  non-axially  symmetric  mechanical  and 
thermal  loadings.  Raju  (18)  studied  the  case  of  shallow  shells 
of  pyrolytic  graphite  type  materials  subjected  to  a  variety  of 
axially  symmetric  and  non-axially  symmetric  loads.  Daugherty 
(19),  (20)  treated  the  case  of  non-circular  cylindrical  shells 
of  pyrolytic  materials. 

The  preceding  deal  with  single-layer  shells.  Anisotropic 
laminated  shells  of  revolution  with  elastic  properties  symmetric 
about  the  middle  surface  of  the  composite  shell  are  extensively 
treated  by  Ambartsumian  (21),  whereas  Dong  (22),  (23),  et  al  (24) 
treat  layered  shells  wherein  the  structure  is  assumed  to  be 
composed  of  an  arbitrary  numbers  of  bonded  layers  each  of  dif¬ 
ferent  constant  thickness,  different  orientation  of  elastic  axes 
and  different  anisotropic  elastic  properties.  Since  Dong  does 
not  assume  elastic  symmetry  about  the  middle  surface,  flexural 
and  extensional  deformations  are  coupled  and  solution  techniques  for 
homogeneous  shells  do  not  carry  over  directly  for  anisotropic  elas¬ 
tic  shells.  Hence,  alternate  methods  of  solution  are  developed. 

Radkowski  et  al  (25)  considered  laminated  isotropic  shells 
of  revolution  with  variable  thickness  using  E.  Reissner's  for- 
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formulation  (26).  Radkowski  extended  this  work  to  include  variable 
laminated  orthotropic  material  properties  (27).  Both  formulations  were 
restricted  to  axisymmetric  loads.  In  Radkowski 's  works  and  that  of 
Sepetoski  (28)*  the  governing  equations  were  cast  in  finite  difference 
form  and  solved  with  the  aid  of  a  digital  computer.  The  introduction 
to  Dong's  paper  (23)  makes  interesting  reading  regarding  the  hazards 
of  this  perfectly  valid  technique. 

Other  treatments  of  laminated  cylinders  have  been  by  Jones 
and  Whittier  (29),  Tsai  and  Azzi  (30),  Paul  (31),  Au  (32),  Keeffe 
and  Windholz  (33).  These,  and  most  other  references  cited  herein 
are  characterized  by  neglect  of  transverse  shear  deformation.  A 
recent  work  of  great  theoretical  elegance,  even  though  it  neglects 
transverse  shear  deformation, is  that  of  Zudans  (34)  which  presents  a  the¬ 
ory  for  arbitrarily  loaded  Mechanically  &  thermally)  shells  of  revolution 
with  internal  masses  and  ring  stiffeners,  derived  under  the  Kirchoff 
Hypothesis  and  consistent  with  balance  of  energy  as  well  as  linear 
and  angular  momentum  and  invariance  under  transformation  of  middle 
surface  coordinate  systems  and  rigid  body  displacements.  The  ele¬ 
gance,  unfortunately,  does  not  carry  over  to  the  computational 
techniques  (35). 

Laminated  isotropic  plates  have  been  considered  by  Vinson  (36) 
who  treated  thermal  stresses  in  circular  plates,  neglecting  transverse 

shear  deformation.  Summers  (3?)  treats  thick  and  thin  isotropic 

and  orthotropic  laminated  plates  including  transverse  shear  deformation. 
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Mehta  (34)  considers  orthotropic  and  Isotropic  laminated  as  well  as 
single-layered  rectangular  plates  of  pyrolytic  qraphite  type  materials 
under  static  mechanical  and  thermal  loads.  Wu  (39),  (40),  and  (41) 
treats  the  lateral  vibrations  of  both  small  and  large  amplitude  for 
rectangular  plates  of  pyrolytic  and  graphite  materials. 
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Except  for  those  references  dealing  with  pyrolytic  graphite 
type  materials,  all  the  works  reviewed  which  analyze  multi-layered 
structures  are  characterized  by  their  neglect  of  either  transverse 
shear  deformation  or  thermal  expansion  through-the-thickness  or 

both. 


Prompted  by  the  absence  of  a  definitive  treatment  of 
laminated  shells  applicable  to  pyrolytic  graphite  type  materials, 
this  thesis  was  undertaken.  It  Is  an  extension  of  McDonough's 
work  In  that  layered  cylindrical  shells,  Including  the  effects  of 
transverse  shear  deformation  and  thermal  expansion  through-the- 
thickness,  subjected  to  arbitrary  axl -symmetric  loading  are 
considered. 
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II.  DERIVATION  OF  GOVERNING  EQUATIONS 


The  coordinate  system  used  is  shown  in  Figure  (II-l).  The 
three-dimensional  equations  of  thermoelasticity  (uncoupled)  for  the 
case  of  axial  symmetry  are  given  by: 

Stress-Strain  Relations  (transversely  isotropic  material) 

ex  =  r(ox  ■  '■vW  +  aT 


ee  =  r  (oe  -  V0X  ’Vz1  +  aT 


ez *  ^  -  r  (ox  +0e>  +  acT 


xz 

xz  ■  257 


(1) 


where  exQ  =  cQz  =oxQ  =  oQz  =  0  by  symmetry 


G1j  ancl  °ij  are  the  Physical  components  of  the  strain  and  stress 
tensors  respectively,  and  for  brevity  a..  -  a.  when  i  =  j. 

‘  J  * 

E»  Ec  ,\>j  vc>  Gc-are  five  independent  elastic  constants  where  again 
for  brevity  v  =  vXQ  =vQx  and  vc  =vxz  =vQz.  a,  ac  are  the  coefficients 
of  thermal  expansion  in  the  x  and  z  directions  of  the  materials.  T 
Is  the  temperature  measured  from  the  stress-free  temperature  of  the 
material  in  units  consistent  with  the  a's. 
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FICURf  n-l:  LATEHEO  CYLINDRICAL  SHELL  COORDINATE  SYSTEM 


6 


K 


The  preferred  direction  for  the  material  is  everywhere 
coincident  with  the  z  coordinate.  This  restriction  is  carried 
throughout  this  work. 

Equilibrium  Equations 

The  equilibrium  equations  when  applied  to  the  present  pro¬ 
blem  become  for  each  lamina: 


R(  1  +  I)  i°x  +  R(1+f) 


3X 


a  a, 


xz 


+  0 


az 


xz 


9c «  y  do 

Z  .  n/i  ,  Z^  XZ 

ax 


RC  +  t>  3T +  R(’ +  f) 


(2) 


The  third  equation  is  identically  zero  ;  ^om  symmetry  considerations. 

Strain-Deformation  Equations 

The  strain-deformation  relations  can  be  written  correspondly  as: 
au„ 


x  ax 


ee  = 


1  u. 


RO+f) 


au. 


(3) 


Gz  =  az 


exz  2  az  ax 


ua  =  0  by  symmetry 
0 
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The  displacements  an;  positive  In  the  direction  of  the  positive 
corresponding  coordinate  (See  Figure  II-l). 

Assumptions; 

1.  Thr.  thickness  of  the  shell  Is  small  compared  with  other 

dimensions,  hence  Love's  First  Approximation  is  applicable: 

l.e.  -  «  1  (4) 

^  min 

2.  The  displacements  are  small  compared  to  the  thickness  of 
the  shell  and  the  angles  of  rotation  are  small  compared  to  unity. 

3.  The  transverse  normal  stress  is  small  compared  with  other 
normal  stress  components  and  is  neglected  in  the  stress-strain 
equations. 

4.  A  linear  element  normal  to  the  undeformed  middle  surface 
undergoes  translation  and  rotation  and  remains  straight,  implying 
deformations  of  the  form 

ux  *  uQ(x)  +  ze(x)  (5) 

5.  Transverse  normal  strain  due  to  thermal  expansion  will 
be  included,  that  due  to  mechanical  (or  isothermal)  loads  will  be 
neglected. 

This  implies  a  deformation  of  the  form 
uz(x,z)  =  w(x)  +  w(x,z)  (6) 

where 

__  z 

w(x,z)  *  /  o  Tdz 
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I 


6. 

With  these 


Material  properties  are  constant  for  each  lamina, 
assumptions,  the  equation  (1)  becomes: 


Ex  =  r  K  '  woe>  +  aT 

e6  T't'  vox>  +  aT' 

(7) 

ez  *  aj 

- 

exz  “  2GC 

Making  use  of  equations  (5)  and  (6)  and  neglecting  z/R  in  comparison 
to  unity,  the  equations  (3)  become: 


e9  *  I  <W  +  W) 

(8> 

_  3W 
eZ  ~  3Z 

CXZ  *  1/2^6  * 

where  (  )'  «o  (  )  =  ^-  (  ) 
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Integrated  Equations 

The  stress  resultants  and  couples  are  defined  as  follows 


'1 


*  / 2  V2 

A 


'1 

~r 


N 


Tx 


'1 


E.ja.jTdz 


17 


T 


h< 


N 


T  0  j 


E^a^Tdz 


T 


M. 


~r 


^o^Tdz 


17 


^a^Tdz 


1  =  a  ,b 
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(9) 


The  equilibrium  equations  ,1n  terms  of  stress  resultants  and 
couples  are  readily  obtained.  Integrating  (2)  directly  yields 


\  +  'll  -  '21  *  ° 


/  h1  hi 

+  t  'n  +  r'21 ' 0 


(10) 


/  N8< 

Q1  ‘  ”R^  +  P11  *  P21  '  0 


1  s  a,b 


Solving  equations  (7)  for  normal  stresses,  first  Integrating 
them  from  -h^/2  to  +  h^/2;  then  multiplying  through  the  equations  by 
z  and  Integrating  once  again  between  the  same  limits,  making  use  of 
the  definitions  of  the  stress  resultants  and  couples  (9)  and 
simplifying,  the  following  stress-strain  relations  result: 
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An  integrated  shear  stress  strain  equation  is  required  in 
addfti  on  to  (11).  The  necessary  expression  is  obtained  using  weighted 
Integration,  the  procedure  being  analogous  to  that  in  reference  15, 
and  for  convenience  given  in  Appendix  D. 


°(  1  T+  thiGc  (Bi  +  V  +  4i 


hi 

where  mi  =  +  t21)  • 


(HA) 


Solving  for  Qg  and  from  the  second  of  the  integrated 

equilibrium  equations  (10)  and  substituting  into  the  third,  making 

use  of  M  and  N  from  (11)  and  the  definitions 
xi  9i 

E.h.  E.hj 3 

C.  ;  D4  (1  -  a.  b) 


1  (1-v,2) 


i 


12(1-Vj  ) 


Tj  *  °zx(-  r  >  *  °zx l+  r  > 


Pj  '  °z  (‘T>  *  °z  (+  T> 


two  useful  relations  are  obtained: 

1  C,V  C  V  h 

DaD  Ba  -  “IF  °Sa  '  T  °Bb  +  T  °'J 


Va  C 
Pj'  Ra  DU°b'2Ra  Wfe 


a  *  v 

17  1 


(12) 
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(13) 


3  "b 

&b  +  2~  ^Tj  +  Pj 


cbvb 


tt2 


Note  that  In  the  above,  use  has  been  made  of  the  conditions 
that  the  laminae  are  bonded  together  and  that  no  slippage  occurs  in 
the  joints  between  laminae.  The  former  is  given  by 


-h  h.  _  -h  _  h. 

s  ^2(2^-)  which  implies  that  wa(x,  -j2-)  +  wfl  =  wb(x,  7^-)  +  wb 

hb  _  ‘  -h 

or  wa  =  wb  +  wb(x,  2~)  -  wa  (x,  — ) 

*  wb  +  E  (14) 


»  _  h  b  _  "h 

where  E  *  wb(x,  -  wfl(x,  — 


1 


The  latter  condition  is  expressed  by 


(15) 


Making  use  of  (11),  (14)  and  (15)  in  the  first  integrated 
equilibrium  equation,  two  further  relations  are  obtained: 
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“37  +  "3 


(16) 


C  h  ry  C  h,  0 

-M-  d2b_  +  -4-4  d2b. 
c  a  i  o 


2  Cava 

t.  +  COY  +  -p-  Dw. 
j  a  oK  R  b 


t.  +  CbD2u  +  Dw 
j  b  °b  Rb  b 


“47  +  n4 


(17) 


Using  ( 1 1  A)  and  the  second  of  the  integrated  equilibrium 
equations  together  with  (14)  and  (15),  two  final  relations  are 
obtained: 


D  D2  6  -  f-Gah  a  +  I,  hi.  -  |  GahaDw. 
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with 


"l  =  "  pla  ‘  2  Dt 1  a 
hb 

*2  =  +p2b  "  ~2  °T2b 


3  =  "Tla 


n4  =  +  T2b 


n5  =  "  12  haTla 


The  governing  equations  may  be  transformed  to  a  more  useful  form  for 
solution.  Since  they  are  cumbersome  for  manipulation  in  their  present 
form,  a  matrix  notation  was  used  to  define  the  coefficients  of  the  un¬ 
knowns.  These  are  tabulated  in  Appendix  C.  After  some  manipulation  of 
the  governing  equations,  the  following  are  obtained: 
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^1°7  +92°5  +93°3  +94°^  Wfc  = 

(22) 

CbllD4  +  b12°2  +  bl3>  6b  ■  *-nCx)  -  (b,/  +  b,50)  wb 

(23) 

a23DSa  +  LIII(x>  '  Ca21D2  +  a22l  V  a24Cwb 

(24) 

“^Oh  *  "-tv  Cxi  -  k1D26s  .  kgD  6b  -  kjDWjj 

(25) 

TJ  '  LV  (*)  ~  V^ob  •  k5Dwb 

(26) 

PJ  *  LVIW  •  k6°38h  '  k7DTJ  +  k8Duob  +  k(,wb 

(27) 
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III.  SOLUTION  OF  GOVERNING  EQUATIONS 


Homogeneous  Solution; 

Consider  first  W^.  Assuming  a  solution  of  the  form  Wb  =  esx 
and  letting  y  =  s^,  (22)  then  takes  the  form 

y3  +  52  y2  +  £3.  y  +  sa  =  0  (28) 

9]  91  9i 

for  the  homogeneous  solution. 

The  solution  of  equations  in  the  form  (28)  is  developed  in 
reference  (42)  and  leads  to  three  possibilities  for  the  roots: 

Case  1 :  there  are  two  conjugate  imaginary  roots  and  one 
real  root 

.  Case  2:  there  are  three  real  and  unequal  roots 
Case  3:  there  are  three  real  roots  of  which  at  least  two 
are  equal. 

Case  1  leads  to  a  solution  of  (22)  which  is  of  the  form 
WbR  =  Vjeci*  +  V2e“clx  +  ec2*(V3cos  c3x  +  V4sin  c3x) 

+  e‘c2x(V5  cos  C3X  +  V g  sin  C3X )  (29) 

where  V-j  -  are  constants  to  be  evaluated  through  boundary  conditions. 
The  constants  c-j  -  c3  are  defined  in  Appendix  C. 

The  Case  2  solution  can  take  on  several  forms  depending  on 
wether  the  roots  of  (28)  are  positive  or  negative.  The  final  forms 
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for  the  case  of  one,  two  and  three  positive  real  roots  respectively, 

are: 


wbh 

=  V1ec4x  + 

V2e“c4x  + 

V3 

cos 

c5x  + 

V4  sin  C5X 

+  Vg  cos 

C6X  +  v6 

sin 

CgX 

(30) 

WbH 

=  V^ec4x  + 

V2e“c4x  + 

V^e1^ 

+  V4e* 

-c5x 

+  V5  cos 

C6X  +  v6 

sin 

C6X 

(31) 

“bH 

*  V^ec4x  + 

y^e-c4x  + 

V  ec5x 
3 

+  V4e' 

■c5x 

+  V5ec6x 

+  V6e“c6x 

(32) 

where,  again,  the  V-j  -  Yg  are  boundary  value-constants  and  c4  -  Cg 
are  defined  in  Appendix  C. 

Case  3  represents  the  degenerate  forms  (30)  -  (32)  where  two 
of  the  roots  are  equal.  The  equations  then  take  the  form: 

WbH*  VleC4X  +  V2e"C4X  +  (V3  +  V5X>  cos  C5X  +  (v4  +  v6x)sin  C5X  <33) 

wbH  *  (Vi  +  V3x)ec4x  +  (Vg  +  V4x)e“c4x  +  V5  cos  Cgx 

♦  V6  sin  c6x  (34) 

WbH  =  ^1  +  V3x)ec4x  +  (V2  +  V4x)e-C4x  +  Vgec6x 

+  V6e’c 6X.  (35) 

In  treating  single-layered  cylinders,  it  is  possible  to 
write  the  c-j  -  Cg  explicitly  in  terms  of  the  physical  quantities 
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Involved  and  thus  have  a  feel  for  the  physical  behavior  of  the  shell. 

In  the  present  derivation,  however,  these  expressions  are  so  lengthy 
and  involved  that,  unless  one  has  a  specific  example  in  mind,  their 
explicit  form  for  the  general  case  is  of  dubious  utility.  A  point  is 
reached  where  one  must  decide  whether  to  obscure  the  physical  situation 
with  mathematics  or  obscure  the  mathematics  with  the  physical  quantities 
involved.  The  former  course  was  chosen  to  show  where  the  mathematical 
formulation  is  analogous  to  that  for  the  single-layered  cylinder 
(Case  1)  and  where  it  diverges  (Cases  2  and  3}.  To  import  greater  phy¬ 
sical  significance  to  the  equations  it  becomes  necessary  to  either  con¬ 
sider  a  specific  case  and  determine  the  form  and  constants  listed  in 
Appendix  C  will  take  or  resort  to  numerical  evaluation  of  the  constants. 


Consider  next  equation  (25).  This  may  be  integrated  directly 

to  give 


(x)  dxdx  -  k^a 


k26b 


-  k, 


ijwi 


dx  + 


V 


+  V 


8 


(36) 


where  and  Vg  are  also  boundary  value  constants.  Since  there  are 
also  eight  boundary  conditions,  four  at  each  edge,  it  follows  that  the 
homogeneous  solutions  for  the  remaining  unknowns  &a  and  B^  are  not 
required.  Their  particular  solutions  will  suffice  to  satisfy  the 
governing  equations.  This  is  equivalent  to  setting  the  boundary  value 
constants  of  the  Ba  and  3^  homogeneous  solutions  to  zero. 


Total  Solution: 


The  total  solution  will  consist  of  the  homogeneous  solutions 
given  in  the  previous  section  together  with  whatever  particular 
solutions  are  called  for  due  to  the  mechanical  and  thermal  loading  for 
any  given  problem.  Total  solutions  for  displacements  will  take  the 
form: 

Wt>(x)  -*  +  Wjj  part 

6t>(x)  *  part 

Ba(x)  =  6a  part 

where  Wb  part  is  the  particular  solution  of  (22)  due  to  Lj(x),  and 

0b  part  Is  the  particular  solution  of  (23)  due  to  LIl(x)  - 

(b1403  +  bi5D)Wb(x),  and 

0a  part  is  the  particular  solution  of  (24)  due  to  Ljjj(x)  - 
(a21D2  +  a22)Sb(x)  -  a24DWb(x). 


(37) 

(38) 

(39) 


Once  (37)  -  (39)  are  known,  uqP  may  be  found  from  (36)  and  the 
joint  shear  and  normal  stresses  from  (26)  and  (27)  respectively. 
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IV.  BOUNDARY  CONDITIONS 


Boundary  conditions  for  plates  and  shells  are  listed 
in  many  sources  (21),  (34),  (37),  (36),  (37).  For  the  axi-sym 
metric  case,  they  are  usually  stated  as: 

At  the  edges  x  •  o  and  x  ■  L 

either  u  or  N  prescribed 

either  w  or  Q  prescribed  (40) 

either  &  or  M  prescribed 

For  multilayered  problems,  the  same  boundary  conditions 
usually  apply  if  N,  M  and  Q  are  interpreted  as  resultants  IT, 

"FT.  Considering  a  two-layered  cylinder  for  example,  laminas  a 
and  b,  the  resultants  would  be: 

■  *,  +  \ 

7  •  qa  ♦  Qb  (4i) 

« •  Ma  +  «b  +  (jji  +  £k)Na 

However,  (40),  (41)  do  not  make  use  of  the  no  slip  -  no 
delamination  conditions.  These  provide  two  constraints  not  only 
on  displacements  but  on  boundary  condtions  as  well. 


Appropriate  boundary  conditions  can  be  derived  using  the 
principle  of  virteal  displacements  of  the  layer  middle  surfaces. 
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(1  *  a,b) 


(42) 


B 


1 


* 


where  the  asterisk  denotes  a  virtual  quantity.  Multiplying 
equations  (10)  by  the  virtual  displacements  (42)  in  the  order 
listed,  adding  the  products,  integrating  over  shell  length  and 
summing  over  the  layers,  we  get: 


DMXi  -  Qi  ♦  1/2  hi  (T11  +  T2<)]  8/  + 

[Q1  -(1/R)N  e1  +  {Pu  -  P21)3  wi*  }  dC]  -  0 

since  the  virtual  work  done  by  a  shell  in  equilibrium  through 
a  vtetual  displacement  is  equal  to  zero. 


After  Integration  by  parts,  the  virtual  work  principle 
for  the  multilayered  cylinder  takes  the  form: 


z 


(  [N 


U  * 

Xi  1 


♦  Mw  8  f  + 


k1 


1«a 


^  C(f  “  \2\  *  V2  hi  (t^  +  t2i )  8^* 
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♦fa  li  "  p2i )  V  3  dx  > 


U. 


* 


+  Hx.  V  +  W +  Wf‘) 


(44) 


The  first  quantity  in  brackets  on  the  left  hand  side  of  (44) 
involves  terms  which  are  specified  at  the  boundaries,  namely: 


QiV] 


L 

o 


Summing  over  layers  gives 


*■.  iL 


EVa*  +  Mx  V  +  Va*  +  VV  +  Hx  •  b*  *  %Wb^ 


S*  v.a  this  derivation  does  not  allow  for  slip  or  delam¬ 
ination,  the  virtual  displacements  must  be  constrained  to: 


ga*  *  u  b*  +  ^ /2  ^ha  6a*  +  ^  /2  ^b  8b* 

wa*  *  V  +  (”b  '  V* 


(  45) 
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Substituting  into  (40)  and  collecting  terms,  the  quantities 
to  be  specified  at  x  =  o  and  x  =  L  are  found  to  be: 

either  Ny  +  N¥  or  uK  specified 
a  xb  b 

either  1/2  haN  +  M  or  B  specified 
xa  a  a 

f  6  ) 

either  1/2  hbNx  +  Mx  or  8b  specified 
a  b 

either  Qa  +  Qb  or  wb  specified 
Qa  (wb  -  wa)  specified 

Note  that  the  first  four  conditions  are  not  a  unique  set. 

Other  oossibilities  are  (ua,  ufal  Ba,  w^)  or  (ua>u  eb,  w  )  soecified, 
but  the  set  (46)  seems  to  be  a  good  choice. 

.  The  last  condition  of  (46)  is  the  result  of  retaining  thermal 
expansion  through-the-thickness  while  dropping  terms  of  order  h/R 
throughout  the  remainder  of  the  derivation.  For  cases  where  w  must  be 
retained,  the  following  procedure  can  be  employed.  For  simplicity, 
free-free  boundary  conditions  are  considered,  though  the  treatment 
Is  analogous  for  any  set  of  conditions  specified. 

Since  the  last  of  (46)  is  a  temperature  dependent,  -  wa 
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Is  a  priori  specified.  Hence,  for  free  boundaries,  either  Qa  =  0 

or  (w!  -  vT  )  s  0.  Qa  =  0  is  acceptable  for  special  cases  but  is  not 
b  a 

generally  true.  (17.  -  w  )  can  be  riqorouslv  satisfied  by  redefining 

d  a 

the  reference  surface  location  in  layer  a,  i.e., 


(I  o 

Wa  (*.  z0>  =  «c  Ta  (x>2)  dz 
Jo  a 

t  hb/2 

Wb  (x,  hb  )  =  I  ac  T  (x,  z)  dz 

2  /  o  b 


Under  these  conditions,  the  no-slip  requirement  becomes 


uo  =  V  +  ('Zo  6a  +  1/2  hb6b) 
a  b 


(•) 


and  all  other  results  remain  the  same  if  1/2  ha  is  replaced  by  -Z0. 

Note  that  when  this  approach  is  used,  the  second  strain 
definition  of  equation  (0)  must  be  used  in  the  form: 

wi  +  "i 

v  jrpnr  (1  =  a-b> 

In  order  to  have  a  zero  stress  state  for  the  case  of  an  isotropic 
material  with  the  same  ac  in  both  layers  and  T  s  constant. 
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V.  SAMPLE  PROBLEM 


The  problem  selected  was  the  slow  cooling  of  pyrolytic  graphite. 
Because  of  the  difference  in  thermal  expansion  coefficients  in  both  "a" 
and  "c"  directions  of  the  pyrolytic  graphite  and  mandrel  material  and 
also  because  of  curvature  effects,  normal  and  shear  stresses  at  the  de¬ 
posit-mandrel  interface  will  be  formed  during  the  cooling  process.  These 
may  be  of  sufficient  magnitude  to  cause  flaking  or  delamination.  An 
investigation  of  their  behavior  with  changing  material  and  geometric 
properties  is  therefore  of  interest. 

The  test  case  considered  a  laminated  cylinder  with  free-free 
edges  and  a  constant  temperature  TQ  =  -1000°F.  The  material  oroperties 
used  were  averaged  in  the  range  3000°F  -  2000°F.  Layer  a,  the  top-most 
layer,  was  taken  to  be  pyrolytic  graphite,  layer  b  commercial  {ATJ) 
graphite.  The  properties  used  are  given  in  Table  CB.l).  The  calculations 
were  performed  with  the  aid  of  a  CDC  6600  computer.  The  program  tab¬ 
ulation  is  given  in  Appendix  A.  Results  showing  the  behavior  of  T.  and 

J 

Pj  due  to  variation  in  lamina  thickness  and/or  E/Gc  ratio  are  in  Ap¬ 
pendix  B.  From  these,  the  following  conclusions  may  be  drawn: 

1.  Figures  (B.l)  and  (B.2)  indicate  a  decrease  in  normal 
and  shear  stresses  at  the  mandrel  -  PG  interface  as  the  E/G  ratio  is 

V 

increases,  implying  that  a  material  weak  in  shear  is  desirable  to  min¬ 
imize  stresses  and  therefore  the  possibility  of  debonding. 
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2.  Figures  (B.3)  and  (B.4)  show  results  when  h  =  0.50, 

E/Cc  *  20  and  5  >_  ha/hb  1  1 •  Figures  (B.5)  and  (B.6)  show  the  case 
where  hb  =  0.25,  E/Gc  *  20,  4  >_  ha/hb  ^  1.  All  the  curves  indicate 
that  a  high  h  /h.  ratio  is  the  desirable,  for  given  material  proper- 
ties,  for  minimization  of  normal  and  shear  stresses  at  the  joint.  This 
implies  that  a  thin  graphite  mandrel  is  preferable  to  a  thick 
mandrel.  Note  that  when  ha/h^  =  1,  the  ultimate  stress  (o^  *  18,000 
psi)  in  tension  for  pyrolytic  graphite  is  exceeded. 

3.  To  determine  the  behavior  of  the  roots  of  equation  (22) 
and  establish  a  range  wherein  the  various  forms  of  solution  ( 29 ) - ( 35 ) 
will  occur,  computations  were  also  made  for  1000  >  h  /h.  >  0.001  with 

—  3  D  — 

E/Gc  at  50,  20  and  2.6.  Results  show  that  the  Case  1  solution  in  the 
form  (29)  occurs  whenever  ha/hb  >_  1  Independent  of  the  E/Gc  ratio. 

The  constants  c^  -  c^  are  affected  by  E/Gc,  the  constant  c^  being 
much  more  sensitive  to  a  change  in  E/Gc  than  Cg  or  c^. 

4.  The  first  two  terms  of  (29)  have  their  principal  effects  at 
the  edge  only.  It  was  also  observed  that  for  hfl/hb  »1,  the  value  c-j 
becomes  so  very  large  that  the  boundary  value  constants  and  tend 
to  become  very  small,  and  as  a  first  appromation  the  terms  containing 
them  can  be  dropped  from  the  solution  functions.  In  this  case,  (29) 
takes  on  a  form  similar  to  that  for  the  solution  for  radial  displacement 
of  a  single-layered  cylinder,  with  suitably  defined  constants. 

5.  It  should  be  kept  in  mind  that  all  the  above  remarks  are 

based  on  the  numerical  results  and  are  valid  in  the  range. 

1000  »  h  /hb*  0.001;  400^L/K^26;  0,05  *  h/R^  0. 0033 
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FORTRAN  IV  PROGRAM  PCE  I  ■  NPUT , OUTPUT . TAPE 5- INPUT , T APE9-0UTPUT »  •••••  | 

OIMENSION  IC0I6I,  IIMI6I.  I  RE  1 6 1  PGE  2 

DIMENSION  OROI 10 1  t  MPllOl,  MQI10I,  COEFIS),  GJAC0<5),  R00TSI2.6I  PGE  3 

COMMON  /SMEAR/  OA.QB.FO  PGE  A 

COMMON  /KIP/  ICI  ,  IC2,IC3,!CA,IC5, ICR, IM1, IM2.IM3, IMA, IMS* (Mb,  IRl,  IPGE  5 

1R2.IRI.IRA, IA5.IR6  PGE  6 

COMMON  01  TEMP  PGE  7 

COMMON  TEMP.EA.E AC, NUA.NU AC ,EB, EMC. NUR, NUBC.HA.HR, R,G)A, Gift, CA.CR. PGE  8 


lOA.Oa.M.CACB.CAPCB.'IUNU.GAMA.GANB.MAOl,  AI55  »,GI  7  I,  ATI  30  )  .  FN  T  OA,  F.*PGc  9 
2TOB.BOEI60OI .lOAOt.AMlBOl ,L0 A02 . ENI6, 7  I ,  F NXA , FNX8 , F  AX  A, FMXB, FNX. FMPGE  10 
SX.FNOA.FNOR.AC () I ,KLlP,CLL,OELPH<  201.VI6I, AL I ISI.MM.OMM, OOnW , MB, OMPGE  t 1 
4B.DOMB,ODOM8,MA,OmA,mU.DMU.  ODMU,OOOmU.TAu.DTAU,MPJ,MmA,Omma, mUA.UwP&E  12 


SUA 

PGE 

13 

EQUIVALENCE  MCl.lCOlin,  IIMl.llMmt.  IIRI.IREIIII 

PGE 

1AE 

REAL  NUA.NUAC.NUR,  VJBC.NUNU.K, LOAOi  .L0AO2 

PGE 

15 

RE AO  <5,121  (OROI  II, 1-1, At 

PGE 

16R 

READ  <5,121  ( HP < 1  1,1-1,61 

PGE 

1  7R 

RE *0  (5,121  <H0<  Ml ,M-1 ,A I 

PGE 

18R 

REA0  <5,111  TEMP.OLTEMP 

PGE 

19R 

RE AO  15, 111  EA.EAC.NUA.NUAC 

PGE 

2  OR 

REAO  <5,111  E6 , E  6C  «NUB , NUBC 

PGE 

2  1R 

REA0  <5,91  BTI1)  ,AT(7)  .BTm.BTlA) 

PGE 

22R 

REAO  <5.101  (DELPHI JI.J-l, 111 

PGE 

23R 

WR1TC  <9,131  ( HQ (Ml, M- 1 , A  1 

PGE 

2AM 

WRITE  <  9, l A  1  TEMP 

PGE 

25w 

WRIIE  <9.151  EA.EAC.NUA.NUAC 

PGE 

26m 

WRITE  (9,161  EB.EBC, NUB, NUBC 

PGE 

2  7m 

COCFIll-50. 

PGE 

28 

COEFI2I.2C. 

PGE 

29 

COEF  (  31  -2.  6 

PGE 

30 

OO  6  llA-1,6 

PGE 

31 

OO  8  H-l,« 

PGE 

32 

OO  8  L-1,3 

PGE 

33 

MA-HP ( | 1 A ) 

PGE 

3A 

MR-HO(M) 

PGE 

35 

R-30. 

PGE 

36 

C3AC0I L l-EA/CCEFILI 

PGE 

37 

G3A-G3AC0I L 1 

PGE 

38 

CjB-c B/ 1  2.  •  < 1,-NueCI  1 

PGE 

39 

CALL  PRELIM  (L.M.IIAI 

PCE 

AO 

CALL  P0LYR  <6, G. ROOTS, 01 

PGE 

A  1 

CALL  RICF  ( KOuTS  , AC.KL 1 P 1 

PCE 

A2 

CALL  THERM 

PGE 

A3 

call  Ml  SC 

PGE 

AA 

OO  7  1  Jr»*l  ,A 

PGE 

A5 

ElL-CRQI 1 JPI 

PGF. 

A6 

IF  (ELL. EO. 0.01  GO  TO  7 

PGE 

A  7 

AMIll-AC(2l«*2-AC(}l*-2 

PGE 

A8 

AMI 2la2.*(AC<2l*AC()l 1 

PGE 

AN 

AM  (  J  1  -AC  (  2  1  •  AM  (  l  )-ACm«AM(2l 

PGE 

50 

AMlAl-ACI  il*AM<  1  1  ♦  A C  (  2  >  •  A  M  2  1 

PGE 

51 

AMI |7I- AMI  |  | -ACI  3  I-3MI2I •*  .<2  1 

PGE 

52 

AMI  I8I-A'M  1  l-ACI  2  l*AM  «ACI!I 

PGE 

53 

AMI  25 1 -AMI  1 1 • AC  1 2I»AMI2»*AC< i» 

PCE 

5A 

AMI  26  1 -AMI l l-ACl  3 1 -AM| 2 1 *»CI 2  1 

PGE 

55 

AM<S>-AC(2>*AM(3>-AC(9>*A«<a) 

PGE 

56 

AMI  6) -AC  111 -AMI  J 1 *ACI2l*AMI4l 

PGE 

57 

IF  IKL1P-5I  1,2.3 

PGE 

5a 

CONTINUE 

PGE 

59 

IF  IRLIP-2I  A, 5, 6 

PGE 

6u 

47 


*  CALL  EKLIP5 

PGE 

61 

CO  TO  A 

PCE 

62 

3  CALL  EXLIP6 

PCE 

63 

CO  TO  6 

PGE 

64 

4  CALL  EKLlPl 

PGE 

45 

CO  TO  6 

PGE 

46 

3  CALL  ERLIP2 

PGE 

47 

*  CONTINUE 

PGE 

48 

CALL  RISC 

PGE 

69 

T  CONTINUE 

PGc 

70 

t  CONTINUE 

PCE 

71 

* 

PCE 

72 

PCE 

71 

PCE 

74 

9  FORMAT  I4E15.7) 

PGE 

75 

10  FORMAT  I6E12.5I 

PGE 

7b 

11  FORMAT  I4E15.0I 

PGE 

77 

12  FORMAT  1  I6F4.01 

PCE 

78 

13  FORMAT  1  AM  MB ■ 1 2F  9. 3 1 

PCE 

79 

14  FORMAT  | 1X.SHTcMP-E12.5I 

PGE 

•  0 

15  format  iix.6maprop-4ei5.ti 

PGE 

81 

16  FORMAT  ( IX.6HBPR0P-4E15.  7/) 

PGE 

a? 

ENO 

PGE 

•  3- 

SUBROUTINE  UC  OFF  1 FNTMN, CONST  .  X) 

1 

COMMON  OLTEMP 

2 

COMMON  TEMP.EA.EAC.NUA.NUAC, EB.EBC.NUB.NUBC.HA.HR. 

R 

.G1A.G1B. 

CA.CB. 

3 

10A,DB.M,CAC0.CAPCB,NUNU,GAMA,GAMB,K(40), AC  551 .Cl  7) 

t 

6T 1 30 I . FN  TOA, FN 

4 

2TOB.OOE 16001 .LUAni.AM(50l .LGA02.EMI6. 7 1 • FNXA, FNXB. 

FMXA.FMXB, 

FNX, FM 

5 

3X.FNOA.FNG0.ACm  ,  KL  1  P  .  E  L  L.  Ot  L  PHI  20  1  .  F 1  6  >  .  AL  (  1  B  1  .  UM 

.GUW.UOriW 

.WR.Dtf 

6 

4B.D0MB , DODMB  .mA.OmA.wU.OkU. UOWU.OOOrfU. TAU.OTAU.MPj 

• 

MWA.UMWA. 

MUA.UM 

7 

5UA 

• 

REAL  HUA.NUAC  .NUO.NIJPC.NUNU.X.LOADI  .L0AO2 

9 

PP4.-IELL/2.  »*ACm 

10 

PP5--IELL/2. 1 *&C ( 2 1 

11 

PP6— IELL/2.  1  •  AC  1  11 

12 

CONST-BUEI 19l«MA«M0E(20l*MB«B0E(21} •FNTNM-BTI20I*X/CAPCB 

13 

RETURN 

14 

ENO 

15- 

SUBROUTINE  PRELIM  IL.M.I14I 

•  •••* 

1 

COMMON  OLTEMP 

PRELI 

2 

COKMuN  TEMI»,EA,EAC,NUA,NUAC,E6.E6C.NUB.NUBC.HA,MB.R.G  3A,G3B,CA,CB,P3EL  I  3 
lOA.OB.H.CACH.CAOCB  ,'luNU.GAMA,GAMft,M*3l ,  Al  55),GC  7 » ,  n  T  I  30  I ,  PNT  G«  ,  FNPAg  1 1  4 
2T0B.BUE  1600)  »  LOAn  I  « A’M  50  I  .L0A02  . E  M  6.  7  I  .  FNX  A ,  FNXO. FMX  A. FM.xB.  FNX,  FM.PTEL  I  5 
3X.FN0A.FN3B,  AC.  (3)  ,  XL  I P  ,  E  L  L ,  DE  L  Ph  I  2  0  J  ,  V  (  6  1  ,  All  l  B  >  ,  w*  ,  OWW ,  00«  *  ,  *B  ,  U«PT  E  L  I  6 
AB.OOmB.ODOwB.mA.OmA, WJ.OWU, OOMU.OOOMU. r AU.Of AU.^PJ.MMA.OrfrfA,  rtUA.OrfPTEL I  7 


5UA 

PAELI  a 

REAL  NUA.NUAC .NUP .NU8C.NUNU.X.LOAOI .L0AD2 

PRELI  9 

CA-EA4HA/I l.-NUA**2l 

PRELI10 

CB«EB*HB/ 1 1. -NUB* *2) 

PREL 1  1 1 

OA*E A* I Mf •• 3 1 /  1 1  2 . • 1 1 • -NUA* *2  1  1 

PRELI  12 

OB *EB*lHil ••31/112. •(  l. -NUB* *21  » 

PREL 1  1  3 

M«HA*MB 

PREL 1  14 

ca:b-ca*cb 

PREL 1 1  5 

CAPC6-C A»CS 

PRELI lb 

NUNU-NUA-NUO 

PREL 1  1  7 

GAMA  *15. / 6.  1 • IG1A«HA 1 

PRELI  lb 

CAMS*  1  5. /b.  )*(G1B*hBI 

PREL  1  19 

XI 15 1 ■ 1 • *CB/C A 

PREL 120 

M 161  -  1. «CA/CB 

PREL 121 

R( 1 71 *HA/HB 

PREL 122 

48 


1 


A I  II  )•!  I./I2.  )  •|MA**2I«M*HA/I6.*KI 16)  I 
A 1 12)>-HA»NUNU/I  I  2.  *R  )*M  16)1 

At  1  ))>(  Me»H6)  •  (*  I  IS)  /  K  I  16  )  )/l2.*H*HR/  |6.*KI 16) > 

A ( 16  I ■MB* A  I 12 )/MA 

Al l 5 ) >- I N/H  A ) »A( 121 

KKJH»(CA*NUA*CA**JUH»  **2 

couh>cap:a>>2 


BOUM> ICA>CAPC8 )•( R**2) 

AU6I  *  (FOU^-LOLU  )  /ft DOM 

A(2H»(H8*>2)/12.» <  5 . /?6.  )  •  (HB**?  )/M IS) 
AI22)>-(5./6. »•«  CGie/6S)*< l.-NU6*»2) ) 


A(7))»(5./26.  )*IHA>HB)/KI  IS) 

AI26)>-|«5./6. I *H8/HAI «A(  121-AI22  > 

AI31 )  >11  1./I2.  )*(hA**2)  )•  II  .♦5./I2.  Ml  16)  )  ) 

Al 12 ) -- I S. / 6. I'llC )4 /£»)•< l.-NUA**2)) 
Ami»IMl5l/miM  )  •  A  (  2  3  I 
A(36)>(S./t2. >»IHA/(R*K( 16)  )  )-CAPA/C A 

A I  6 1  I  «  31  11)  AI51)-BI21> 

AI62I  ■  B(  12)  AI52)  •  Bl 22) 

A I  6  3  )  ■  Bill)  AIM)  *  6(231 

A 1 66  I  >81161  Al 86 )  «  Bl  26) 

A165)  ■  BI18)  AIM)  >  Bl  2D) 


A(61)>A(13)>A(21)“A(  1 1  )  *  A  I  2  1  ) 

Al 6? ) * A  I 16)*A(23)-A(  11)*A(22)~A(21 I *A I  121 
AI63I--AI  1 2 ) *  A  I  2  2  I 
A(66)«A(16)*AI23)-AI  1 1 ) *  A  I  2  6 ) 

AI65)»A(16)>A(23)-AI  1 2  I  *  A ( 26  I 
AI5ll>AI2)l**2”A 111) «A(21  ) 
A(52)»-AI)l)*A<22l-A(32l*A(2lJ 
Al  53)  »- A  (  32  )  *  A  (  2  2  I 
AI56|-A|?)I«A( 36 ) - A l 31  I • A | 26 ) 

A(55I>-A02)«A(26) 

chi  -  oin 

Cl  3 )  -  0121 

Cl  5 )  •  DO) 

CI7)  >  0(6) 

Clll*AISLI«A(66)-A(6l)«AI56) 

C(3)>A01)«A(65)»A(S2)«A(66l-A(6?)*A06)-A(6l)*AI55l 

CID)»AIS2)6A|6S)6A|S3I*A(66I-A162)*A|D5)-A|6))6AI56) 


C(7)»AI53I*A165)-A|6)|*AI  55 ) 

Cl  21 *0.0 
CI6I.CI2) 

C(6»«C(6) 

WRITE  19,3)  L.G3A.G3B 
WRITE  (9,6)  HA.HB.A,  TE"<» 

WRITE  19,5)  I I0P.GI (OP) , I 0P>1, 7) 

WRITE  (9,1)  CA.CB.OA.OB 

WRITE  19,6)  CACR.CAPCe.GAMA.GAHH 

WRITE  19, B)  HA,HB , ELL ,1, Al 1 7) 

WRITE  (9,71 
K 1 1 )  *  1, 0 
KI2)  *GO)/CI  1 ) 

K|3>«G(5)/GI l ) 

R( 6) »G( 7)/GI I ) 

ERUKP-I  ).»II))-*(2)*»2IM. 

HRHP«  I  2  .  •  I  M  2  I  •  •3I-9.«IK(  2)  «KOI  )»27.»K|6)  )/27, 
MUMP  •  DM  ALL  A 
HU)  >  SMALL  B 
RATIO-HA/HB 


PRELI23 
PREL I  26 
PAELI25 
PREL I  in 
PRELI27 
PRELI2B 
PRELI29 
PRELI30 
PRELI31 
PRELI32 
PREL133 
MHt)6 
MEL  I  35 
MEL  I  36 
MEL  I  3  7 
MEL  I  3d 
MEL  I  39 
MEL  16  0 
MELI61 
MEL  I  6  2 
MEL  I  63 
MEL  I  66 
MEL  I  6  6 
MEL  I  66 
PREL 16  7 
MEL  I  6  8 
MEL  I  69 
MEL  I  50 
MEL  I  5  1 
MEL  I  52 
MEL  I  5  3 
ME  L  I  56 
MEL  I  55 
MEL  I  56 
MEL  I  5  7 
MEL  I  58 
MCLI59 
MEL  I  60 
MEL  I  6  l 
MELI62 
MEL  I  63 
MEL  I  66 
M  EL  I  65 
MEL  I  66 
MEL  I  6  7W 
MEL  I  66W 
MELI69W 
MEL  I  70W 
MEL  I  7iw 
MEL  1729 
MEL  I  7  3 w 
MEL  I  76 
MEL  I  75 
MEL  I  76 
MEL  I  7  7 
PRELI 78 
PREL I  79 
PRELI80 
PREL I  8  1 
PAEL I  82 
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OISC-(MAMP**2l/4.*(FRllMP**3l/2  7.  ME l !  63 

WRITE  19,21  Ml, MB, RAT 10, 01SC  PREII84W 

01  SC ’TMINANT.LT*  0  IMPLIES  RLIP  -  2.3,4,  OR  5  OEPENOING  Vi  RPR S L  I  85 
01 SCRlMlNANT.EQ.O  IMPLIES  DEGENERATE  CASE  PREL  1 86 

OISCRIKI NANT.CT.3  KLIP  -  l  MODIFIEO  CLASSICAL  SOLUTION  PRELI8T 

RETURN  PRELI8B 

PREL  189 
PRELI90 

1  FORMAT  I  t*.3HCA»El2.S,l*,lMC6-El2.S,l*.  1H0A-E12.5, IX. 3H0B-E1 2.5/ »  PREL I  91 

2  FORMAT  I  //  l  *  ,  1HHA.F10.  1,  l*,  3MHB-F  1 0.  3 , 5  A  ,  6HH  A/ MO-  F  l  S.  7,6X,  i  3H’J|  SCRPR  EL  I  92 

1ININANT-E16.7//)  PREL  I  93 

3  FORMAT  I2K.4HI.3A  II  2, 2H|  *c  IS.  7.4MG38  E15.7/I  PREL  I  94 

*  FORMAT  I2K.4M  MA«F|3.2,2X,4M  HB-F 1 0, 2 , 2 K  ,  2HR «F 10. 2, 2X , 6H  TEMP-FIO. PREL 195 

12/  I  PREL I  96 

5  FORMAT  (4(2X,2HC(I7,2H)-E12.5I/3(2X,2HG(  I2,2M)-E12.SI)  PREL 197 

*  FORMAT  I2X,5mCACB-E12.5, l  X ,  6HC  APC  B- E 1 2 . 5  ,  IX, SMCAMA.E 12. 5,  IX,  5HGAMBPKEL I  98 

1-E12.5/)  PREL I  99 

7  FORMAT  I20H  All  I  ARE  I  PREL  I 

*  FORMAT  IT 12.5.F12.5.FI2.  5,1  3, IX, 7M  MA/HB-F12.5)  PREL I 

ENO  PREL  I  - 

SUBROUTINE  Ml  SC  *  •' 

COMMON  OLTEMP 

COMMON  TFMPtEA,EAC,NUA,NUAC,EB,EnC,NUB,NUBC,HA,HB,R,C3A,r.3B.CA.CB, 
10A,0B,H,CACB,CAPC!3,NUNU.GAMA,GAMri,K  (401  ,  A I  SS  I ,  Cl  7  )  ,  H  T I  30  I  ,  FN  T  0  A  ,  FN 
2T0B.B0EI 6001  .LOAD l ,AM| SOI ,L 0 A02 , E N I  6 , 7 ) , FNX  A , FNKfl , F MX  A,  FNXtf, FNX , FM 
JX.FNOA.FNOH.ACm  ,  KL  I  P  ,  E  L  L  ,  JE  L  PM|  20  I  ,  V  I  6  )  ,  Al  I  l  8  I  ,  WM  ,  uWW  ,  OJWW  ,  *B  , 
4B,OOWO,OODWB,WA,JwA,  wU  •  0<4U,  OOWU,  OOLMU,  TAU,OTAU,WPJ,WtfA,OWMA,MUA,OM 
3UA 

REAL  NUA,NUAC,NUP,NUBC,NUNU,K,L0A0l.L0AO2 
•OEI  19I«CA*MA/I2.*CAPC0» 

BOE(20)-IHB/MA>*805<  191 
BOEI 21  I •< 1. /R)*ICA*NUA*C8*NUB)/CAPC6 
BOE ( 22 1  -CB 

•OEI23l»"OEI22»PNUB/R 
BOE 124 )«(mO«*2 I* AOEI  221/12. 

B0EI2SI-MM/2. 

B0EI26)»60E(22I/|R**2) 


K(1)«60E( 19) 

Rl 2 1 -BOE I  201 
KI3I-BOEI21) 

R  (4 )  *BOE I  221 
RISI -XI 81 -BOE 123  I 
K I  61 -HOE  I  24) 

R I  7 ) -BOE I  251 
Rl 93 -BOE I  26) 


AMI24)-ACI2)**2*ACI3)**2 
BDEI 37) -CAMA-CAMB 
BOEI3M3-AOE I  3 7|-l M/12. 1*606123  1 
L0A02-D Till) 

LOAOl-O.O 


BOE  I  6 1  ) -6DE  I  I  9| *POE I  27  I 
BOEI62)-M06I20I*BUEI22) 

BOE 163) -BOE I  21 1*6 OEI 22  I -BOE  (21) 

BOE 164) -ODE  I24I*80EI2SI*(BOEI22)*BOEI20)  ) 
BOE  1 65 l-BDEl  75)* I BOE 1 22 1 *8DET  »9)) 


*  l 
2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 
S3 

34 

35 

36 

37 

38 

39 
60 
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■43  €  j  66) «MOE  1251*806(611.  47 

BOEII>7)-BDEI23)«60E(19)  42 

■OE  4  681 »*0E 421)*0OE(2O)  43 

BOE(69t<rAOE(26)-eOE(2))«0DC(21)  44 

ANI32)'»«CA*.NUA»C8*NUBI/A  45 

AHmt-CA«HA/2.  44 

AMI )4)»CA«HB/2.  47 

AM|35)»OA*(M*MA)  «<CA/4.  44 

AM(36)»0«MH*Hfl)  *CA/4.  4, 

)3>«AM(  341  50 

AH I  30  I ■  AMI  171 *NUA/R  31 

AMI39)>GAMA»GAMB  32 

■N(42)-0*«<CA/4. |4<HA«*2I  53 

AM(4})*|MA*MBI«(CA/4. »  34 

B0E(tl)-AMI 31l*NUA/«  35 

BOEI8?)»AM( 34>*NUA/R  54 

AMI44I-Ar|42)-AH(1)I«60E( 191  57 

AMI  43) -Ar|43)-AM( 11) «A0c 4 20 )  39 

AMI  46) -COE  4 «l  1-4*4 J3I«RUE 1211  34 

AM|4T)-0«*ICA/4. I •|HB*«2 )  40 

AM|48)-Ar.<4  t)-AM|  J4MBOE4  191  4| 

AM(49)-A»I4  M-AMl  34)  *BDE  (201  42 

AM(30)-BGE482)-A*M14)»m0E(21)  43 

BOE (631* 1*4/ 12. (•(B0EI22) «B0E(  19))  44 

■DEI84I-IH/12.  >«inOEI22)«BOEI20>)  43 

BOE  483) -AMI  39) «(h/ 12. )*IBOE  4  22  )*BOEI21)-BOEI23)»  46 

BOEI71I-AM()3I>CA«60E(19)  47 

BOE  172)  *AM(  J4)-CA*B0EI?0)  49 

BOE  4  731  ■  I NUA/ R-R0E  121)1  49 


BOE174I-CA*('4U**RTM)/R*RT(20I/CAPCB>-BTCS1/I  l.-NUA) 
BOF I  73  I ■— CB*B0E I  19) 

BOE  I  76  I  ■— CB^BUEI  201 

BOE 177) »CB* I NUB/ 9 -BOE 121 ) ) 


BOE  (91 1 -AM|44MA«(48I  74 

BOE  4  92) -AMI  43  l*AM(  49)  75 

BOE  4  93 ) *  AMI  46 ) ♦AH ( 30 1  74 

77 

RETURN  79 

EMO  79- 

SUBR0UTINE  THERM  MM»  1 

COMMON  OUT EMM  THERM  2 


COMMON  TEMP,EA,EAC.NUA,NUAC.EB,EBC,NUB,NUBC.HA,HB,R,G3A,G1B.CA,CB,THERM  I 
lOAt09,Ht;ACB.CAPCB,HUSU»CAM4,CAM0tM40)f  A  4  33  ),  Cl  7 ) , BT ( 30 ) , FN T0A , FN Tm= R M  4 
2TOB  tSOE  I  (00  )  .L0A01.AMC30)  ,Lj402fEN<6,n,FNXA,FNXB,FMXA,FMXB.F:N<,FMTMERl  «» 
3*»F NO A . F**3B  •  AC  4 1 ) .fLIP.cLL.JELPh(20).V(6),»LI  1 8 ) <  WW . OWM , ODrfrf • MB »  D*  THE  <  M  6 
4B.OOMB  »03©MB.  HA.  DMA.  »U , 0  MU, OUWU , DODmU  ,  TAU.  OT  AU  .  WP  J ,  MM  A  ,  EiRWA.  WU  A ,  CM  Thc  KM  7 


-  THERM  0 
REAL  NUA.NUAC.NUB,NUBC.NUNU.K,lOAOl,LOAD2  THERM  9 
BT(l)  •  AAPH44A)  THERMIC 
tf(2)  -  AAMMA4AC)  THERM  1 1 
•  74  3 1  •  AAMHA4B)  THERM12 
8744)  -  9APHAIBC)  THtRMli 
B7I3)  •  NTXA  -  *4 T OA  THERM  1 4 
BTIbl  •  NfXB  ■  NT08  THERM  1 3 
BT4  7)  -  BAR - 4 )  AT  2  -  -HA/2.  THERM  1 6 
BT ( B I  •  M4BAR-BI  AT  L  •  ♦Hb/2.  THERM l 7 
8*1 9)  •  61  BAR  I  ■  *<4A*-6>  -  M ( BAR- A )  THCRMld 
BT (10)  ■  ALPHA  lit/)  THERM  1  t 
BT  4 1 1  ;  »  ALPHA  4  2.7)  THER*<?0 
B74 121  »  MTXA  TMERM21 
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•  TU3I  -  fcTXS 

BTI14I  -  EXPANSION  IN  THICKNESS  DIRECTION  -  A 

•  T I  1M  •  EXPANSION  IN  THICKNESS  OIRECTICN  -  8 

•  T I  16)  -  MCOIFIEO  AXIAL  THERMAL  RESULTANT  A 


TWERX22 

THEAM23 

THERM24 

TMERM2S 

THERM26 

THERM27 

THERM28 

THERM29 

TMERM30 

IHERM31 

Therm J2 

THERM33 


•Tim  -  AXIAL  THERMAL  RESULTANT  (MODIFIED)  -  • 

•  T I  IS  I  *  8TI14I 

•t 1 19)  -  stiisi 

IT  I  SI ■ IE  A*HA I • I TE MP*BT  1 1 1  I 

•  T|Al»IEN*MBI*ITEMR*BTI2l  I 

•  T I  TI*-BT(2)*ITEMP*HA/2.  I  ♦DLTEMP*! BTI 2  )*HA/B. » 

•TIBI-BT  141  •(  TEMO«hM/2.  I «0LTEMP*I8T |4)*MB/B.  I 

•  TI9I-BTI  8I-6TI7I 

•  Tl  101*8 TISI/IR*! l.-NUAI I -C A*BT I  9 ) / I R *«2 »- 1 OA/ 12. AHA  I )•!  OLIEMP/f R*TH£RN34 

1*2  I  I *B  T I  2  I  THEAM3S 

•  Tl  11I-BTI6I/IR* I l.-NUBI  I  -  I  08/ 1 2. «HB I  I • I  01 T F MP/ ( R*«2 > ) *8T I  4)  TMERM36 

•Tl 121-1 TE  HP*  8  T llll*(EA*HA**3)/( 12* *H) ♦  I E A*MA** 3  I •( OL  TEMP*BT llll/l THERMIT 

I12.4HAI  THEAN38 

BTI  IJI-I TEMP*ST( 3 ))• I EB*HB**3 J /I  12. *R I ♦< EB*H6*« 3) •! 0LTEMP-8T I 3II/ITHERM39 
112. *HBI  THERM40 

•  Tll4l-OA*NUAI*ITEMP*PTl2l  J/R  THERM4  1 

•  TI1SI-ID9*NU0»*ITEMP*8T(4I l/R  THE  RM42 

•Tl  16I-8IISI/I  l.-NUAI-CA*|NUA*BTI9)  |/R-IOA*DLTEMP/l2.*HAH*INUA*8TTrtERM43 

1I2I/RI  THERM44 

BTI17I-BTI6I/II,  -  NUB  J-IDA*DLTEMP/(2,*MB)  )  *1  NUB*RT  I  4  l/R  )  TH  E  RM4  S 

•  Tl  1BI-I0A*NUA/R I *I8T(2I • TEMP  I  THE  RM46 

■  Tl  1 91-1 OB*NUB/R l«(MT|4|*TEMP|  THERM4 7 

BTI20I-6T1  161*81 1  171  THERMSd 

•  TI21I-I UA* IOL  TF  M0*MT(2II/I2.  •HA»I*INUA/RI  THERM49 

■  T I  22  I -I 0B*( 0LTEMP*8T (41  I /  I  2. *MB M • I  NUB/ R I  THERMSO 

PP3«A<23)«A(S)I  THERMS1 

PP10-AI32I  •AI4  3I -All  2  I *4(83  1  THERMS? 

ADUN>IPP3/CAI*IBTI  I'JMRTI  11  I  I  THERM5  3 

•oum-ica*nua*;8*nu8i/ica*ir*;apcbi»  thermsa 

COUH«PP 1*|6T(  161*971  1711  THERMSS 

DOUM-BDUH-COUM  THERMS6 

EOUM-ADUM*OJUM  THERM57 

FDUM-CA*MH/ICb*C APCO »  TMERMS8 

OOUM*  I  F  OUM*  ( S. /  1  2  »  I  I  •  <  PP  l  0*9  T  I  S  I  /  I  l.-NUAI  >  TMERMSN 

H0UM-IS./12.  I  *IHR*«»P10/CAPCB)*BT|  16  1  TMERM60 

PDUM«EOUM-OUUf1*HOUH  THERM6I 

OOUM* I PP 13/LCI*(UM13I-BT(I8I I  THERM62 

ROUM-POUMPQOUM  THERM63 

BOE 1491 -ROUM/OI 7  I  THERM64 

MR  I TE  19,11  THERM 65H 

AOUM.BOUM.COUM, OOUM, E DUM, F OUM.GOUN, HDUM , POUM  ,  RUUM , BDE I  THE  *M66* 

THI RH67W 
THE  RM60W 

IBTII»,I-1,I0I  THERM69H 

TMERM70 
THERM71 

1  FORMAT  I 1X.3SHC0MP0NENTS  OF  B0EI49I  ARE  />  THE RM7? 

2  FORMAT  I1X.1CH  BT  ARE  I  THE  R  M  73 

3  FORMAT  I1X.6E12.SI  THERM  74 

A  FORMAT  I1X.4EIS.7I  THERM  7S 

CNO  THI RM  76- 

SUBROUTINE  tKLlPl  •*••*  l 

01  ME  NS  I ON  J I M I 1 1  I  EKLIP  2 

COMMON  /SHEAR/  QA.OB.FQ  EKLIP  3 

COMMON-  OLTEMP  EKLIP  4 

COMMON  TEMP.E  A  ,E  AC,*<UA,NUAC,EB.C8C,NU6,NU8CfHA,HB,RtC3A,C3B,CA,CB,EKL  IP  S 


WRITE  (9.41 
1491 

WRITE  (9,21 
WRITE  (9,31 
RETURN 
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10A,DB,M,CACB,CAPCft,NUNU,GAMA,GAMH,A(40l,  AI5S),r,(  7  I , B T ( 30 > , FN I OA , FNE K l  IP  6 
2106, BOE  IbOO)  .LOAD  l  ,  AM  I  501  .LOAD2.ESI 6,  7 I  ,  F  NX  A  ,  F  NXB,  fMX  A,  F  M  t  b,  FNX  ,  fhE  4,1  1 1»  7 
3X,F90A,E‘408,AC(l>,Kl!P,ElL'J£LPH<70l,VlM,ALllri),  wm ,  t,-w ,  00  ««w ,  MS  •  0*E  K  L  I  P  8 
AB,DOWB,OJOWB,  MA,  OwA,  HLl.OwU,  DDMU,OOD«J.TAU.OTALI,WPJ,WWA,L)4WA,MtjA,D.<E«llP  9 


5UA  EKLIP1G 

HEAL  NUA.VOAC  ,  NUB  ,  NIIBC  .  NU9U  ,  K  ,  Li  AD  1  ,L0AJ2  EKLIP11 

REAL  JIP  EKLIP12 

P4«ACIll*ELL  EKLIP13 

P5- AC  I  2 ) • £  L  L  EKL I P 1 4 

P6«  AC  I  3  I  •  E  L  L  EKL I P 1 5 

BDE  (  ll«4|44l*IACI  1)**I)4AI45)*AC(1I  EKL I P  1  6 

BDE I  2  )  *  A  I  4  1  I *  (  AC  I  1  )**A|4AI42)*IAC(  1 )**2)*A|43)  EKL I P  1  7 

AMI30I-B0E  (  1  )/BOE  I2»  EKL  I  P  I  b 

B0E(3I*AM(5I*AI4I  >*AMI1)*A|42)«AI4}|  EKLIPI9 

B0£J4|*A*'(6I*AI41  |«4MI2I*AI',2)  EKLIP20 

BDE I 51 « API  3I*A|44I«AC(2I*A(4SI  EKLlP^l 

BOE I  61 -AM <4  l*A(44  |  ,4Ct il »A(451  EKL  IP22 

BDE  I  71 -BOE I  51  ‘BOF  (  31  ♦ *0t I b >  *BOE I M  EKLIP2  3 

BDE  I  HI  -OOF  ( b)  «eJJt  I  W  -POE  I  5i  •P0EI4I  EKL I P24 

BOE I  9 1 «  9PC (  3I**24HDE |4|**2  EKLIP2S 

BDEI  101  *f!OE  I  7)/ACfc(9|  EKL  I  P2b 

BDEIll)*C0£<8l/B0EI9>  EKL I P2  7 

AMI  7)*- At  2  4 | • ACI l  I  *AM( JO  I •(  A(  2 1  I *4  AC  I  I I*  *2  *♦ At  2211  Ekl I P2b 

AM(8)'A(?3I*( ACI l l**2l  EKL I P2  l 

BDEI 121 -A*!  7  I  /  A  h  (  h  I  EKLIPiC 

AH(19|*A(2II«AH(IM«(22|  EKL  I  P  )1 

AM(20)»A(21 I • API  2  I  EKL I p  32 

AM(21I-A(?4)*AC(?I  EKL I P  3  3 

ah(22i*a(24»«ac(  ii  Ekl i p 34 

BOE  mi  «  AH  I  |9I  *BOf  (10)-4M|20>*00E<  III -AM  |  21)  EKLIP3S 

BDE ( 1 4  I  *  AH  I 1 9  I *B0c  IIll4AH(20l  *601 (1JI-AMI22)  EKL IP36 

BDEI IS) -POE  I  1 3)/A 123)  EKL  I  P 3  7 

B0EI16I-CCEI 14I/AI23I  EKLIP38 

AM ( 9 ) -A  M  I  | I  •  •  2 ♦  A*H  2  ) • *2  EKL I P  39 

AMI  101  *AW(  1  |  •BOE  I  151  »AH(  2I*BDE  I  lb  I  EKLIP4  0 

AMI  11 I -AM  I  l >  *bOE (  lb) -AH I  2 1*806 115)  EKL IP4  I 

6DEII7)-AN| 1JI/AHI9I  EKLIP42 

BDE I 1 81  ■ A  M I 11 I/AHI9I  EKLIP43 

AMI  1 3 1 ■ -BDE I  19)*BDEI  1214906 1201 *AH| 30  I -BOE I  2 1 1 /AC I  II  EKL I P44 

BOE  I  271  -00c  I  191*1  AMI  2  I  *BDE  I  IBI-AMI  1 1  *BOE  1171  > *BDE I  20  I •( AM(l) *BOE  I  1EKL I P4  5 
10I-AMI2I  *80c  I  It)  l-«0EI2U  *ACI2  I  E  KL I P4  b 

BOE I  78  I --BOt I  191 •  I  AH|2)*80E  I  1  M  4AMI  1 ) *BOE I  I B 1 1 *BOE( 20 )*|  AH|  1  )*BDEI  EKL |P4 7 


111 |4AMI2l*BOE I  10  I  I  -*0E I  2  I » *AC  I  3  I 
BDE I  291 -AMI  I  I *dOEI 27  l*AMI2)*BOEI2») 
BOE  I  301 -AM( 1  I *BOE I  2ft  I -AM  I 2I*B0EI 27) 
AMI  141-ODE  I  291 /  A  M  (  9 1 
AMI 151-BOE I  30) /AMI  9) 
BOEI31)»EKP|P5)*COS|P6) 

BDEP32I  »cXP|P5)«SIN(Pfel 
BOE I  331 -5  API -P5 1 *l OS  I Pb I 
BOE I 3A)-ExP(-P5l*S!NI*b) 

BOE (51) -ACI 1 1 *(ma*BOE 1 12 )-Hft*AM( 30) I 
BOF I  5?) --BUE I  161 *m A 4  POE  I  1 ll*M8 
BDE I  53) -BOE I l 7) *MA-BOE I IOI*MB 
BOE ISA  I -AC  I  1 1 *00E(S2l4AC(2l  •30EIS3) 
BOE I  551 -ACI 2  I P00EI5?)-ACI 3I-8DEI53I 
AMI  1 2  I -B  r  <  20 1 /CAPCB-6DE 1211*8 f (9) 


BDE 1701 -L0A02*6DE 1 23 )  *BI I  201/C  APCB 


EKL I P46 
EKL I P49 
EKL IP50 
EKLIP51 
EKLIP52 
EKLIP53 
IKLIP54 

Ekl i P55 

Ekl )P5b 
EKIIP5  7 
EKL IP53 
EKL IP59 

Ekl iPbo 

EKLIPbl 
EKL I *62 
EKL I P6  3 
EKL I Pb4 

Ekl ip65 
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BOEI  78I»CB»8TI20I/CAPCB-8TI6I/I1.-NU8I  EKLIP66 

00  l  1*1,6  EKLIP67 

VIII-0.0  EKUP60 

00  1  J*  1 «  7  EKL  I  P69 

EMI  I  «  J )  *0.  0  FKLIP70 

1  CONTINUE  EMIP71 

EMI  1,1 )>4C(  U  >*9DE  I  121-  AMI  451*  AMI  JO  I  I  ♦  AM  14  b)  EKL IP  72 

EMI  l «  2  )  *cN( 1 ,  1  I  EKL I P  7  I 

EMI  1 ,31  *AM44  I  •  I  AC  I  ?  I  •  II  71-ACI  )l  •B0EI18I  I  ♦ AMI  491  •!  ACI 3I«BD£I  1  I  lEMlP  74 
l-ACI2)*mcI  IDM*AM(4M  EMIP79 

EMI  t  ,4)  « AMI  44  »M  »Cm*30E  1 1  7 1  ♦  AC  I  2  I  «80E  I  181) -AMI  451*1  AC  I  3I*80E  I  10IEKLIP76 
l*AC  I  2  I  *BUE  I  1 1 1  I  EM.  I  P7  7 

ENI 1,9I*FN(1,3)  EKLIP7B 

EMI  l,6)*-ENIl,4>  EKLIP7* 

EN(2.1I>a:| l  I  • ( SOE  I  121*4M<48|-AM| 30  I ***149 1  I ♦AMI  SOI  EKLIPbO 

EMI  2  •  2  )  *E  N(  2 , 1 1  EKLIPHl 

EMI2,3l»*Hl<«H)*|-aCIJI*9OEIl8)*AC(2l*B0EI17)l*AMI49|*(ACI3l*B0EIUEKLlPo2 
1)-ACI2I*PUEI  10)  MA“I  SOI  EKLIPHl 


f  M I  2  » a  )  .  A  *•  1 4  n  )•!  AC  13)  *BOE  1 1 7 1 ♦ AC  I  2  )  *PDE I  18) ) -API  49 > • I  AC  I  3 ) *8DE I  1 0) EKL  I  PH4 


1«AC(?I*B0E  III)  I  EKLtP89 

EMI  2  *  9 )  »c  'll  2  »  3  I  EK L  I  P 86 

CN12.6) *-rNI2 ,4|  EKLIPB7 

EMI  1, 1 1  »!•  DE  l  «•>  I*  AC  1 1  l-AHI  10  I  •  I  ODE  I  84l«( ACI 1  I  ••  2  I  ♦  GANB  I  ♦BDE  I  1  2  I  •  I  BOE  Ml  PB8 
IE  I B3 I  *  I ACI  1  I  **2  I  ♦CANAI  EMIP89 

EMIT. 2). -EM  3.11  EKLIP90 

EMI  3.3)  *00S  I  t>9  )  •  AC  I  2  )  -BDE  1 10 1  •  I  C  A*B  ♦ODE  I  34  I  •  A*  III  I  *«UEI  1 1  )•!  BOEI  B«>cKL 
1I*AM(2)  I  ♦HOE  I  l7|*|G4MA*A.Mll)*B0clH3ll-BJEIlrt)*II'0;l83)*Af*l2l  I  EM  IP 7  2 

ENI  3,4|*A:i  31  *036  IB9I-BJF  1 1 1)  *lG4PD*B3El  8414AM  I  )  I-30FI  1  Jl*l  HJEI  84EMIP9I 
lltAMI?)  I*BC- I  t  II  •  I  CANA*  AM  1 1  |«B0c(d3)  M8UEI  l  7  I  •  (  RU€  I  6l)*AfH  21  )  Em  |P7«. 

EMI  3 »  9  )  *  -  £  M  3.31  EKLIP95 

ENI  3,6l*»-M3,4|  EKLIP‘76 

ENI4»l)*EXP(P4)*EN|l  ,1)  EKLIP97 

EMIA,?l»E*PI-P4»l*eNI  1.21  EKHP98 

EMIN, 31 *ENt  1  ,  3)*P JEI  II  I  -  ENI 1,4 ) •BDE I  32)  EKL IP99 

EMI4.4l.rNI  l  ,  !I*PL’EI  32  I  *cNI  1,4  I  *Brit  I  II  )  EKL  IP 

ENl4,5)*rN|  toiepilEI  3II-ENI  1,4  )*eOEH4l  EKLIP 

EM|4,6).C‘I(1<SI*Pl)E(  34  I  ♦  c  N(  l ,  6  I  •  80b  I  111  EKLIP 

ENI5, l)«E*PIP4|*FNI?,tl  EKLIP 

EMI  9,2) »t<PI~P4| •  c  ’ll  2,2  I  EKLIP 

EMIS,3)*ENI2, l»*POE<  U  I -ENI 2,4  I •HOE! 32)  EKLIP 

EN|5,4)«EN(?«S)*PUEi  l?)*tNI2,4|*0DEIill  EKLIP 

ENI  9,91  *t‘ll  2, 9 1  *808133  I -ENI  2,6  I  *4CE  I  34)  EKLIP 

CNI9,6J*:NI2,9>*B0Ef  34  )♦  (ENI  2,6  l*BOEI  II)  EKLIP 

ENI6,l)»t'XPIP4|#t‘Hl,l)  EKLIP 

EMI6,2I*E*PI-P4)*F'H  3,21  EKLIP 

ENI6,3)>*NI  3,  3  I  *  BDE I  It  I  —  £  N I  1,4  I • HOE  I  3? I  EKLIP 

ENI6,4).£M  I,  1 1  •  PDl  I  i?  I  ♦  E  NI  3,4  )  •RUEI  31)  EKLIP 

ENI  6, 9 1  «*NI  I,  9)*  BUS  I  331  -  ENI  1,6  I  «3l)Et  34)  EKLIP 

ENI 6. 6)> ENI 3,3)*80EI  34  )  ♦  £  N<  3, 6  )  *BOE  I  33  >  EKLIP 

EKLIP 

EMI  1.71*1*1  31)«BTt?0)/CAPCB*BOEIB|  I  *0  7 1  9|*BTI  181- I  HA/2.  I*  IBM  91/ I IEKL  IP 
1,-NUA  )  )  —  1  T  I  12  1/11.  -*JUA  )  *  I  HA/?.  |  *HT(  21  )  ♦ AMI  49  I • B3C I  4-11  EKLIP 

ENI  2, 71*  AN  I  34|*B  r  I29I/CA^C8»i30E  142) «3T|9)*mT  I  l9)-|M«/2.  )  •  I  BT  I  9) /I  IEM  IP 
l.-MUAI  1-rtTI  1 1 1/1  l. -NUB  )  ♦  I  M0/2.  )  *8 1  I  221*  AN  I  90  I  *Bi)e  149)  EKLIP 

EMU,?)— ENI1.7I  EKLIP 

ENI?, 7)—  ENI2, 7)  EKLIP 

ENI 3,71 >LGAiH  EKLIP 

ENI4, 7)«FMI  1, 71  EKLIP 

EMI3, 7I*EN(2. 71  EKLIP 

ENI 6,7I*ENI 1, 71  EKLIP 


54 
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(RU» 

MITE  19,11  (IDE  I t ), 1-1, 100)  IM.I*  W 

WAITE  19. M  CKL 1 R  W 

MITE  19.41  I  AA(  JI.J-I.90I  HUE  H 

MITC  (9.101  CKLIR  W 

WRITE  19,4)  (  ICNI  I  ,  j)  ,  J.  1  ,  7  I,  I- l.b)  IUIE  W 

CALL  FMGLE  I  E  N,  V  I  f  U  IT 

WRITE  19, n  I  I ,VI f t, l-l . A)  CKLIR  W 

CKL  ir 

CKL I R 

ARC  12I-BTI20I  /CARCA-BOE 171 1 *8061 49)  CKL  I R 

RR I -AC I  I  I -e'LL/2,  CKL  I R 

RR2-AC(2l«FLL/2.  CKL I R 

RR1-AC (  1 1 *f LL/2.  CKL  I  R 

CALL  CHECK!  I  V  ,AC «  R* l , RR?  «RR)  I  CKL  I  R 

MITE  (9,71  (  I  ,V(  I  I,  1-1 ,6)  CKL  I  R  W 

CKL  I R 

ALIII— AH|10l»vm  CKL  I R 

AL ( 21 -4*1101 • V(2 I  *  CKL  I R 

Allll— BOEI l  0  I *V  t 1 1 -ODE ( l  I  I  «V I  4  )  CKL  IR 

AL (4) -BOE ( I l  I  »VI)  l-KOE ( 10)»V(4 I  CKL IR 

ALII)  »B0E  (  1JI  •  V  ( *»  I  —HOC  (  1  1  I  •  V(  6  I  CKL  I  R 

AL<6)*ao((ui«vm"0Eii0)«v(b}  cklir 

al i  t i *nt>E i  I ? i * v i  i  i  Cklir 

ALIBI— B0SI12I*V(2I  CKLIR 

ALi9i>noc(i/i*v()i*aoE(iai*v(4)  cklir 

AL(  101— AUEI 18I»V()I*60EI17)»V(4)  CKLIR 

ALI  11  I— AOEI l 7I«V(SI-RDE I  1 tf I • V ( b  )  CKLIR 

AL  I  121— BOEI  im»V(5l-RDEI  17KVI6I  CKLIR 

CKLIR 

CKLIR 

ALini*«H(lll«Vin  CKLIR 

ALI14I  —  AHI|1)*VI7>  CKLIR 

AL(li)«A«(l4l*V|ll-4*llll»VI4|  CKLIR 

ALl  1 1* )  —  AH  I  I )  I  *  V  f  )l-AH(l4l«V(4l  CKLIR 

ALII7I«-AW|14)*V(SI«AW(1'»I*VIAI  Cklir 

ALI  l 9  1  — AMI  IS  I -V< 4|-*H(  14 )«V(6)  CKLIR 

MITE  (9,41  CKLIR  W 

MITE  19,91  (ALII),  I  - 1 ,  IS)  CKLIR  W 

CKLIR 

DO  *  KR.l.U  CKLIR 

K-OELRHUFI  cklir 

IF  (X.GT.FLL)  GO  TO  2  CKLIR 

Rl-AC  I  l  I  •*  CKLIR 

R2»  AC ( ? I *K  CKLIR 

R)-AC())*X  CKLIR 

CKLIR 

CKLIR 

CALL  OIFF  |R1,R2 » R 1 1  CKLIR 

CKLIR 

CKLIR 

WW-VI l ) •£  KR( R 1 ) ♦ V  1 7) •EIR(-RI)»EXRIR2)*I9I9)*C0S(R))*VI4|*&I9IR1I  l*EKL  IR 
IEXRI -R7I  •(  VI  ■>  t -COSICD-VI  6) -SlNIH)  l-BOE  149)  EKLIR 

0 WW-VI ll-BOEI10ll*VI2l-aGE(llll«V<ll*BDE 11211 *VI 4 ) *B0E I  111 l»VI*)*BEKl IR 
IDE! 141 l«V(6l*00C  (  111  )  EKLIR 

OOWW-VI  IIMJE  (  102)«VI2I*AOE(  I12)»VI  J 1  *BOE  1 1 27  I  *VI 4  1 •  HOC (  1 121  ♦¥( 5 l*CKL  I R 

iBoci i42i*viA)*bDfii>7t  cklir 

WB-AL lll*EKR|RI)-All2l*clRI-R| I -E XR ( R2 I • I AL 1 1 1 *C0SI R1 1 ♦ AL I  4 1 • S INI REKl | A 
11)I*E*R(-R2)»!AL  (M*C0SIR»)-ALI6)»SIN(R1I  I  CKLIR 


OWB-ALI  1  )4  0CE(  10LI«ALI2I  »BOEI  1 11 1*  ALI  1 1  *BDm  7 1 1 »  AL  1 4 1  *BOE  I  111  I «  AlEKL  I 
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1I3I*B0EI  131I*ALI6)*RDEI131>  EKLIP 

OOWB'ALI  1  ) *ODE 1 1  02 ) ♦ AL I  2  I  •  BOE  I  112  MALI  3)  ‘BOE  I  I  22  )  ♦  AL  (  A  »  •  60E  (  1  32  I  ♦  AE  K  Li  P 
1LI5MBDEI  132  MALI  6  MODE  I  1521  EKLIP 

OOONB>AL  ID  *BOE(  103)  *ALI  2  MBDE  1 11 1  MAL 13)  *BOE<  12  3  I  ♦  AH  4  )  *RuE  I  l  33  I  ♦£  KL  I  P 
1ALISMBDE I 143)«tL(6l *BOE I  151)  EKLIP 

WAa  AL  ITME  XP(  PI  )  vALl  HMEXPI-Pl  MEXPt  P2)  *1  ALI  <>»  *C0SIP3  MALI  10MSINI  EKL  IP 
IP3IMEXPI-P2MIALIU  MCOSIPTMAU  1  2  M  S  I  N  I  P 1 )  )  EKLIP 

DMA- At  l  M  *60E  (  10  l  DAUB)  *B0fcl  l  II  DAL  I  9)  •ttOEC  12  l  MAL  I  1 0  MCDE I  1  J 1 )  ♦  AE  KL  I  P 
lLltlMBDEI  131MALII2MR0EI131)  EKLlP 

DOWA  •  A  L  I  7  I  »BOE  1 102MALHH  *BCEI  1 12  )  ♦  AL  I  9  )  *60E  I  l  22  )  ♦  AL  (  10  )  *BOE  (  1  32  )  ♦  6  KL  I  P 
(All  II  MBOEI  l32MALtl?MB0E(lS2>  EKLIP 

OODMA  ■  AL  I  IMBjEI  1  0  3 )  ♦  Al  IB  )  •  BOE  I  1  11)  ♦  AL  I  1  I  »BOE  I  12  3  MAL  <  10  I  *BOE  I  1 3  3)  E  KL  I  P 
I ♦ AL 1 1 1 ) • BOE 1 13  31 ♦  AL  * 12  MBOEI  131)  EKLIP 

EKLIP 

FNTMWa( 1. /ACI  1)1*1  VI  1  ME  X  PI  PI) -VI  2)  *EXP1»P1 )  Ml  EXP  I  P2>  /  All  I  2  A  )  )*1  VI  CKL  IP 
131*1  AC  12)  *C  OS  IP3D  AC  (!l*SINIPl)l*VH.|*tAC.I7l*SINIP!l- AC  (1)*C0$(P3>£KLIP 
2)>aIEXPI-P2I/A.MI2a))*(VI3)*(-AC(>)»C0SIP1MACI  3MSINIP3I  I“VI6)*I  AC  EKLIP 


H2MSIN1P1MACI1MCQSIP3I  )  I  ♦BDEI  *X  EKLIP 

IF  IKF.E'7.1)  CALL  UCOFF  I  FnTnw .CONS  7*0)  EKLIP 

HU*- BOE  I  191 *wA-BOE I2  0MWR-B0E I  21 l*F NTw*a I 8T I  20  )/CAPCB ) •* ♦CuN  S  T  EKLIP 
DMU*-BD£ I  19 )*D*A- POE  120) aOw«-POE( 2 1  I *w«  *  0 T 1 20 )  /C APCB  EKLIP 

ODWU«-BOC  II  <l*OOMA-*OEI20l*UO<<B“BOEI?l  )  *Dw*  EKLIP 

OOOMU--BUE I  19  1 •O0DW4-P0E I  20 1 •ODO«B- ROE  I  2 l MOO**  EKL IP 

TAU* L0A01 -80c  122  1  *00wl -BO  6 ! 2 1 ) *0* *  EKLIP 

OTAU»-BOE I  22  I ♦OOCwu-BOc ( 231 *DOHM  EKLIP 

MPJ«L0AU2-B0t  (23)  •  0003B-80E(25MOrAU*8DE12>)  *DwuaBOEI  26  MW*  EKLIP 

¥UA«WU«HA*4A/2.  •HB*.B/2.  EKLIP 

DWUA«0WU4MA*0WA/?. ♦mb»9wR/2.  EKLIP 

OOWUAaODMU»NA  *00* A/2. *hB *D0WB/2.  EKLIP 

CALL  RSLT  I x)  EKLIP 

CALL  PRIM  IXI  EKLIP 

2  CONTINUE  EKLIP 

RETURN  EKLIP 

EKLIP 

EKLIP 

3  FORMAT  I1X.9H  0OE  ARE  /I10E12.3II  EKLIP 

R  FORMAT  I  10E 12.3)  EKLIP 

S  FORMAT  I  2 OH  AMU)  ARE  I  EKLIP 

*  FORMAT  (IX.  IE  1 1 • 3  I  EKLIP 

1  FORMAT  13  1 IX, 2HV I  I  1,2H)-E  12.31/21 1X.2HVI  [ l, 2HI-E12. 3) I  EKLIP 

•FORMAT  I1X.2CH  ALII)  ARE  I  EKLIP 

•  FORMAT  I 6E I  2. 3  I  EKLIP 

10  FORMAT  1 2 1 H  ENII.JIBY  RJwS  ARE  )  EKLIP 

END  EKLIP 

SUBROUTINE  EKL I P 2  •••♦•  I 

01  MENS  I  ON  JIPIll  )  EKLIP  2 

COMMON  /SMEAR/  OA.OB.FQ  EKLIP  1 

COMMON  DLTfc.MP  EKLIP  3 


COMMON  TEMP.EA.E  AC  »NUA  »  NUAC , E  B  »  r RC»  NUB.  NU8C*  HA  »  MB* R, C  1A, C  Id. C*.CB, EKLIP  3 
IDA.OR.H.CACB.CAPCO  »NUNU*  0  AMA,GAW8  ,  K  I  30 ) , A  I  33  I  ,  Cl  l)»“.  Tl  30). FNTOA.FN  EKLIP  o 
2T0B, BOE  I  6001 .LOAD l, AM|3J», LwA 02, tNIb, 7), F NX A, FNXB.FMXA.F MX B.F NX, FMFKLIP  I 
31. F NO A, F NOB, AC I  1 ) .KL I P.ELL.LELPMI 2  0 ), VI  b ) , At  I  1 «  )  . * W , 0 , UOWW , *B . OW t KL I P  H 
AB.OOHB.UOOWB.MA.UWA.WU.  UmU.OOMU.OOOmU.T  AU.Or  AU.WMJ,WtoA,ut.rfA,«UA,DtaEKL  IP  9 


SUA  EKL I P 1 0 

REAL  NUA,NUA;,NUB.NUSC.NUNU.K,10ADW10AD2  eklipu 

real  JIP  EKL I P l 2 

P3»AC I l I • c  L  L  EKLlPll 

P3«AC I  2  MF  LL  EKL I P  13 

Pb* AC  I  3 ) •£  LL  EKL I P  1 3 

•OEM  MAI  33  I  •(  AC  I  1  Mal>-AI33l*ACt  D  EKL  IP  16 
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■DEI2I-ACM)*<  ACI  l  I  <*2  >*<ACI  1  >**2I*AI43)  EKLIP17 

AMI  JOI-BDfU  I/80E  (21  IKLIP18 

bde  oi-am hi  •»(* i » *am<  n  •*<  *2 1  ♦*!  4 ji  E<a ip  19 

BDE  m  ■AMIM'AI*  l  I  *AMI2I  •»«*.?»  EKLIP20 

6DE(5)»AM(  3)*At*M»ACI2)*AI<.M  EKIIP21 

BOE  (M«AM(4  l*A(^^  I^aCC  1 1  •  A I  *5  I  EKLIP22 

BO£I7)»BOE<5>*BOEI 3) *BOE ( 61 •BOEI* I  EKL I P2 3 

BOE  I*)  »«0EI6I  »60E  I  II  -BOE  I  5) *BDel*  »  EKL  l  P2<> 

BDE (9>«BDE  I  3  I  ♦  a0E I M  **2  EKL  ?P25 

B0EI10I-P0EI  n /BOE  (9  |  EMIP26 

BOCim»90EI8)/BUc(9|  EKLIP27 

AM{7)>-A(2M»AC(l>MMno)«(Ai2n*<ACtll**2l-A(22>r  EKLIP28 

AN|BI*A| 231*1 ACI  1  »  •*  2  »  EKLIP29 

BOE 1 1? ) *AM 71 /AH  I  8 )  EKLIP10 

AMI  19) »A I  21 > *AM(  1)*A 1221  EKHP31 

AMI  20) -A( 21 )*AM| 2 )  EKLIPJ2 

AMI  21  >"A|  2M  *ACI  2  I  EKLIP33 

AMI22)»AI2<>)*ACI  31  EKLIP34 

BOE 1 13) -AMI  191 *0OE  1 1  0  )  -  AM  ( 20 )  •  BCE  1 1 1)  -  AM  I  2 11  EKL  I P35 

BDEIIM*AM(  lvl*BOE(l  l )  ♦AMI  201*806  I  10) -AMI  22)  EKLIP36 

BDE I15)*80EI13I/A(23J  EKLIP37 

BOE 1161*9061 IAI/AI2I)  EKLIP3U 

AMI9)>AMI|)**2«AMI2) «*2  EKLIP39 

AMI10)*AMIl)  «HOE  I  16)  ♦AMI  2  )  *80£  1161  EKLlP<>0 

AMI  II  I -AM I  1 1 •  BOE  I  16)  -AM|  2  )*8DC< 15)  EKLIPM 

B0EI17I»AMUCI/AM(<*>  EKLIPM 

BOEI  181-AMI  11  )/A.-'<9)  EKLIPM 

AMI13IX-BOE)19|*0OEI  l 2  I ♦ BOE 1 20  I • AMI  »0  I -BOEI  2  U  I ACI  1)  EKLIPM 

BOE  127)  ■BOtl  1 9)  •  I  AMI  2  I  *BUc  (  1 B  I  -  AM  I  1)  *BDE  I  1 7 )  I  »Bi)E  I  20 )  •  I  AM  I  l )  *BDE  1 1 E  KL  1 1»*  ■» 
101-AP.I2l*B0c(ll»)-3r)tl21l*ACI2)  EKLIPM 

BOEI  281  *-306  I  19)*  I  AM|2I**i0EI  1  7  )*AMI  1 1  *BDE  I  IB  )  I  ♦BOE  I  20  )•<  AM  |  1  I *BDE I E*L I P6 7 


111>*AM?)*BUEU3I  )-•()£  <2  l  I* AC  I  3) 

BDE I  29 1 . AMI 1 1 *806  I  2  7  I ♦ AM ( 2 1 *806 I  2 B  ) 
BO£IJC)»AM|1)«buEI23)-AMI2I*BOEI27) 
AMI  1A)-8DEI29)/4M(9| 

AMI  1 S I «BOE I  30 I/AMI9) 

BOB! 3l)»EXPIP5l*C0SI P6) 
B0EI)2)»t*P(P5)*SINIP6) 

BOE  (33  )  *e  XP I -P5 1  *C OS  I  PM 
BOEI 3M »EXP(-P5) *SINI P6) 

BOEI 51 ) «ACI l )*(MA«BDE| 12 )-MP*AM| 30) 
BOEI52)  — 30EI  tai*H**C0?l  l  1 )  *H8 
BOE151)«BOEI17I*ma-BDEI 1 ) ) *MB 
BOE  (5M«ACI3)*80E(57I*AC(2)*80EI5II 
B0EI56)»ACI2l»af)E(52)-4C  i  3)  •BOEI  53) 
AMI12)>6TI20t/CAPCB-80EI21)*BTI9) 


AMI32)-(CA*NUA»CB*NUB)/A 
AMI  33) «CA»HA/2. 

AMI  34  I *C**M8/2» 

AMI 3SI*0AMH*HA)*CA/4. 

AMI  36 1  *0M  I  H*H6)*C  A/4. 

AMI  1 7 )  •  AM |  )  1 1  » AM  i  34) 

AMI  38) »AM< 3  7 1 • NU A / A 
AMI  39) aCAMAtwAMS 
AM|A?)»0A*(CA/9« I *IMA««2  I 
AM(4ll-(HA*H6)*(CA/4. I 
BOE I  M )  “AMI 33I*MU4/R 
B0E(B2)“AM( J4)*NUA/R 


EKI IP4B 
EKLIPM 
EM.  I  P5  0 
EKLIP51 
EKLIP52 
EKLIPS3 
EKLIP64 
IKLIP56 
CKLIPS6 
EMIP57 
EKLIP5B 
CKLIP59 
EKLIP60 
EKLIPM 
EKIIP62 
EKLIP63 
EKL1P64 
EKIIP65 
EKLIP66 
EKLIP67 
IKL I P68 
EKLIP69 
EKLIP73 
EKUP71 
EKLIP72 
EKLIP7  3 
EKLlP/4 

EKL 1 P  75 
EKLIP76 
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A«144|-AH{42I-AM(33I*60EI19I 
AN!  45) -AM  4|J -AH  (11)  *1)06  (  20) 

AMI  46) -BOE <81  )-AM<ll|*80E(2l) 

AHI471-0BMCA/4. I*|HB*»2  I 
AM|4d)-AK(41»-AN( 14 | -BOE  I 19) 

AH|49)-AM(47)-AN(  141*806(20) 

AN (60) -BOE (82  I -AM  14  I  *806  1211 
BDE(H3i>(H/12.  I  *  I  «C*:  1 22  1  ««0£  (  191 1 
60EI84)-(H/12«  )*(Bl)?<22)*B0£(20) ) 

BOE I  85) -AM ( 39  )  ♦ ( H/ 12 1 1 • C  HOE  1 22 ) *80E I  2 1 ) -BOE1 2311 

00  1  1-1.6 
vm-o.o 

00  1  J-1,7 
ENI I, JI-0.0 
CONTINUE 

ENIl.n-4Cm«(AH(44|4B0E(12)-AM(4B|*AM|30l  I4AMI46) 
EN 1 1 .  2)  -C.  0 


ERLIP7? 
ERLIP78 
ERL1P79 
IRLIPBO 
ERLIP81 
ERLIP82 
ERLIP83 
ERL I P84 
ERLIP85 
ERLIP86 
ERLIP87 
ERLIPdB 
ERLIP89 
ERLIP90 
ERLIP91 
ERL1P92 
ERLIP93 
ERLIP94 


ENI1.3)-AM<44>*IAC<2)«80£(17)-AC(3I»BDEI  18 ) > ♦ AMI  45 > •( AC (1) *BDE < 11IERLIP95 
I-ACI2  )*BDE  I  101  )*AM46>  *  ERLIP9* 

ENI1,4>«AM(44I-<AC())*B0E(17>*AC(2)*60EI 18>)-AM145)*I  ACI3IPBDEI IOI6RLIP9  7 
1*AC I  2  I *BOE (11)1  ERLIP98 

EN(1,5)«ENI I, Jl  ERL1P99 

EN(1,6)--EM( l .41  ERL  IP 

ENI2»1)-AC(1(*<BDE(  1 2)  *ANI48  )•  AHOO  l*AN|  491 1 4AHI 501  ERL  IP 

ENI2. 21-0.0  ERLIP 

EN(2«21 -CNI2. 1 )  ERLIP 

EN(2.3)-AN(4b)4(-AC(3|4B0E(  181*ACI 2I*B0E 1 171  )*AHt 491*1 ACI 3)*OOE(11ERLIP 
1  )-AC I  2 1 -FOE ( 10)  ) *AH( 50 1  ERLIP 

EN(2,4|-AMI4d)-(ACIll -BDc 117)4 AC (2) -BOE ( 18 ) I -API 49)*( ACI  3) -BDEI 101  ERL  IP 
|4AC 121 *BOE 111)1  ERLIP 

ENI2.5I-ENI2.3)  ERLIP 

ENI2.6I— ENI2.4)  ERLIP 

ENI3. 11-0.0  ERLIP 

ENI3.2)-BDE(M5)*AC(1)*B0E 1 121*  < BOE 183)*ACI 1 1 ••2-GANA)-AM( 30 )• I BOE I  ERL  IP 
lB4)*AC<l>*-2-GAM8>  ERLIP 

EN13.3>-?0EIB5I-AC<2>-B0E(10I*(GAN8*6DE(84>*AN| 1 1 1480EI 1 1 1  *1  BOE I 84ERL IP 
II4AHI2) 1*8061 I7)-(GAMA*AM<1)*BDE(83))-BDH IB  >*  I  BOE  I  <)1)*AM  2 ) )  ERLIP 

ENI  3*41  -AC  I  )l-80E<84|-ft0E(lll*(GAM6«BDE(»4|*4M(l>|-l)OE(lGI*IBDE(64£RLlP 
1)*AMI2)  l*BOE<  181 -(GAMA-AMm -BDEI  83)  HBOE 117  >*I6DE<  831-AMI  2))  ERLIP 

ENI3.S)  — ENI3.3)  ERLIP 

ENI 3.61-ENI 3.4)  ERLIP 

ENI4, 1 ) --NI l , l I -COSI P6I  ERLIP 

EN|4«2)-ENI 1.  1 1  *  S I N( P6 )  ERLIP 

ENI4.3I -ENI  l, it *006131 1 -ENI  1.4 )*BDEI 32)  ERLIP 

CNI4.4I-ENI  1.  1I4B0EI 121 'ENI 1,4)*B0E( 31)  ERLIP 

ENI 4. 51 -ENI 1,6 ) *006113 1 -ENI 1.4 l-BOE I  34 1  ERLIP 

PN|4V6)-ENI1,3)-B0E<34)*EN<1,6I-I!DE<33)  ERLIP 

fcf)IS.tl-fcNI2,l)*C0SIP6l  ERLIP 

EN(5*2)-£NI2.1I-$IN(P6)  ERLIP 

ENI 5»3)-EN(2f3l-RDE(3l)-6NI2»4)*BDEl32)  ERLIP 

ENI 5*4  > -ENI 2,3)*60E(32)-cN{2,4)-RDE(3I)  ERLIP 

CNI5tSI-ENI2t5)-ODE<31l-eNI2,fc)*BOEI 14)  ERLIP 

ENI$.6)-tNI2,5)-B0E(34)«ENI2,4)*8OEI33>  ERLIP 

ENI6.1)  — ENI3,1I«SIN(P6)  ERLIP 

ENl4*2)-cNI3,l)-C0S(P6l  ERLIP 

ENI6.3)-£N<3.3)-60E( 31)-ENI3,4)*BOEI32)  ERLIP 

ENI6*4)-ENt 3,3)-ROE<32)*ENl J«4)-m>EI31)  ERLIP 

ENI 6,*) -ENI ),5)-B0E(13>-£N(3,4l«6O£t34)  ERLIP 

fNI6.6)*ENI3*5)-B0EI)4|*cNI3*4)-80El»3)  ERLIP 
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EKLIP 

EN(l»7)*AH(33)*BT(?0)/CAPCB*BOEI8l)*0T(9)*BT| |B)- I  HA/2. ) »( BT I  3) / 1 1 EKL IP 
l.-NUA)  )-0T(  121/t  l.-NUAM  (HA/2.  I  *BT  (21  )  *AM ( 46  I  *BOE  I  A  9 1  EKLIP 

CNt2'?l«AM(14)*M (20)/CAPCR«ft0EIS2)  *BT(9I«M(  191-1  HP/ 2.  ) *1  BT ( 3 ) /( If KL  IP 
l.-NUAII-PTI 131/1  1  •  -NUB  )  ♦  (  HB/2.  I  *BT  ( 22  I  *AM(  30  »  *BOE  <  A9I  EKLIP 


EN11,7).-EN(1,7)  EKLIP 

ENI  2 *  7)  —£N(  2  »  7)  EKLIP 

ENI 3,7)*LOAOl  EKLIP 

EN|4,7)*EN( l. 7)  EKLIP 

ENI3,7)*ENI2, 71  EKLIP 

ENI 6. 7) *EN( It  71  EKLIP 

EKLIP 

WRITE  19.1)  (BOEII). 1*1,100)  EKLIP 

WRITE  (9,31  EKLIP 

WRITE  (9,41  I  AMI  J),J*1.30I  EKLIP 

WRITE  (9,121  EKLIP 

WRITE  19,6)  I  IEN(  I ,JI ,J* 1 ,71,1*1,6)  EKLIP 

CALL  FINGLE  (EN.V)  .  EKLIP 

WRITE  (9,71  ( I ,V( I ), 1*1,6)  EKLIP 

EKLIP 

EKLIP 

AMI  12)-er(20l/CAPCB-B0E( 21l*bOE(49)  EKLIP 

PPl-AC(l)«ELL/2.  EKLIP 

PP2«AC(2)*£LL/2.  EKLIP 

PP3«AC(3)*ELL/2.  EKLIP 

CALL  CHCCK2  (V.AC ,PP1 ,PP2,PPJ)  EKLIP 

WRITE  (9,71  (I,V(I),I*1,6)  EKLIP 

EKLIP 

ALI1)*AMI30)*V(1 )  EKLIP 

ALI2I*-AV(3.  I  *V(  2 )  EKLIP 

ALI3I— BDt(10)*V(3)-B0E( 1 ll*V(4)  EKLIP 

ALIA)  *6l)t  ( 1 1 )  *  VI  D-9CFI  10)*VIA  )  EKLIP 

AH3)»BUC(lCI*VO)-RGE(lU*V<6)  EKLIP 

ALI6)«BUE( 11I*V(  3)*P0E(  1 0  I • V(  6 )  EKLIP 

ALI  7)«-0Oc(12)*V( 1 )  EKLIP 

ALIA>*90F (12) *V(2>  EKLIP 

AL I  7 1 *BOE 1 12 1 *V(  I )  EKLIP 

AHH)— Ai)E(12)*V(2)  EKLIP 

AL(9l*OOE(l7>*VI3)«BUE(ld)*Vl4)  EKLIP 

ALI10)*-ADE(1!*I«VI3) *COE I |7)*V(4)  EKLIP 

ALI11I— A0E(17)*VI3MACJ(  l8l*V(ol  EKLIP 

ALI  12I*-BOE(1BI*V(3I-BOE( 17)*V(6)  EKLIP 

EKLIP 

WRITE  1 9, # I  EKLIP 

WRITE  (9,9)  (ALI I  I ,1*1,1?)  EKLIP 

EKLIP 

00  2  KF.1,11  EKLIP 

X«0ELPH(KF)  EKLIP 

IF  (X.GT.ELL)  GO  TO  2  EKLIP 

PI* AC ( 1 ) *X  CKLIP 

P2-AC (2) *X  EKLIP 

PJ*AC( 3 ) • X  EKLIP 

EKLIP 

EKLIP 

CALL  OIFF  (P|,P2,P3)  EUIP 

EKLIP 

EKLIP 

WW*V(l)«C0StPll*V(2t *SIN(P1I»EXP(P2I*(VI 3)*C0SIP3)*V(4)*SIN(P3)I*£EkLIP 
1XP|-P:»*(V(3) *C0S(P3)«V(6)*SIN(P(I I *RUE(49>  EKLIP 

OWN* VI 1 )*BOE  (  191  I *V( 2  I *HOE ( 19  I  I  ♦/( 1  )*nO£ 1 121  I *V I  4  I *BOf I  I II  I *V  1 31 *0*KL  I  P 


AMI  12)«BT(20I/CAPCB-80E( 21 l*B0E(49) 
PP1*AC( 1 )*ELL/2. 

PP2«AC(2)*ELL/2. 

PP3«AC(3l*ELL/2. 

CALL  CHCCK2  (V.AC ,PP1 ,PP2,PP31 
WRITE  (9,71  (l,V(I).l-1.6> 

AL I U  *AM ( 10 ) • VI l ) 

ALI2)— A*(3.I*VI2) 

ALI 31— B0t(10l*V(3!-B0E( 1 ll*V(4) 
ALIA)  *Bl)t  ( 1 1 )  •  VI  D-9UFI  10)*VIA  I 
AL(3)»BUC(lCI*VO)-AGE(lU*V<6) 

AL I  6 ) *BUE ( 11)  *V(  S I ♦POE ( 10I*V(6) 

ALI  7)«-0r>c(12)*V(  l ) 
ALIA)*80F(12)»V(?> 

AL  I  7 1  *6(>E  I  121  *  V  *  1  ) 

ALIH)— Ri)E(12)*VI2) 
ALI9l*OOf(l/>*VI3l*BOf(lS)*VI4| 
ALI10)*-BDE(19I«VI3I  ♦  COE  I  17)*V(4) 
ALI  11 1  — ROE  1 1 7)* VI 3) ♦RCJ( l8l*V(ol 
AL  I  I2I*-6i)E  ( 19 1  *V(  3I-60E  (  1 7 1  •  V  (  fc  I 

WRITE  (9,81 

WRITE  (9,9)  (ALI I  I ,1*1,1?) 

00  2  KF.1,11 
X«DELPH(KF! 

IF  (X.GT.ELL)  GO  TO  2 
PI- AC  1 1 ) *X 
P2*AC(?)*X 
PJ-ACI 3 ) • X 


CALL  OIFF  (P1,P2,P3) 
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10C(14ll»VI*)*60E( 151)  ERLIP 

OOMMaVI II*B01I182>*VI2I*A06(192)*VI1) *906 1 122 1  *V( 4  )*60E ( 132) *V( 5 1*6X1  IP 

iB0€(i42i»v<6i*eo6ii62i  erlip 

M0«AL  111  *C0SI PI)  ♦4l(2)*SlN(Pl)*6XP(P2l*IALI  3)»C0StP3l*AL<4)*SiNtP3ERLIP 
1»I*6XPI-P2)*I AL(5)*C0S(P)I*AII6)*SINIP3I  )  ERLIP 

DKB*ALI  1  1*806  1 181  HAH  21 *606(191  I  * AL I  3 ) *A0C ( 12  11 *Al I A 1*806 1 1 31 )♦ At ERL  I P 
1131*8061  141  HALI4H406  1 161)  ERL  IP 

00MB*  AL I  1 1 *HOE  ll82t*4L(2l*6UEIl92)*ALI 3  I *60t 1 I  22 ) ♦ AL I  4 1 *BO£l  132)*AEKL IP 
Ul  31  *806 1 1421*  ALIA  1*8061 162  I  ERL  IP 

OOOMB-ALI 11*8061 143MALI2I*  606(1911*  ALI3)*6D6I  1 23 1  ♦  Al  1 4  I  *806 1  1 31  HEXl  I P 
1ALI 61*806  I |4 J  |  *41(61 *606 » 15  31  ERL  IP 

MA-ALiri*C0S(Pll*ALI8l*SlN( PI  I *t XP ( P2 I •( AL 1 4  )*COSI P 3)*Al ( 101 *SIN(PERL I P 
11I»*EXP(-P2I *IAL( 1 1I*C3S<P3HALI 12)*StM(P31 >  ERL  IP 

OWA*AL I 7)*B3Ellftl)*4LI8l*B3Ell9l)*ALI9l *MDC I 12ll*ALI 101*6061  lllHAEXLlP 
ll(tll*R06tl41l«AL(l2l*AQEI151)  ERL  I P 

OOWA* AL I TI*60e(lH2HAL(8)*SD£( 192 1 ♦ AL  « 9 1 *BOE ( 1 22 » ♦ AL 1 101 *BOE 1 1 32 l«ERL I P 
IALIUM0O61 142  MALI  12 1  *406 1  1621  ERL  IP 

HOOKA ■ AL I 7 1 *806 1 1 831 « AL( 6 1*606 (1931* AL (91 *4061  123 ) *AL I  10 1 *BOE I  1 311 EKL I P 
1«ALIU)*B0E(143I*ALI  I2I*B0E  (16  31  ERL  IP 

ERLIP 

PNTWM*ll.2AC(lll*<VI ll*SIN<Pl)-V(2)*C3S<Pl>> * I  6 IP ( P 2  I / AH  I  24) I  *  I  V I 36RL I P 
1|9|ACI2I*C0SIP3I*ACI 3 1 *S I Nl P) I > *V 14 ) • ( AC ( 2 1 • S  I  Nl P3 t-ACI 3HC0SIPII (ERL  IP 
*I*IEIPI-P2)2AM(24I >*|VI  61 *I-AC I2)*C JSI  P3HACI  3HSINIP3I ) -VI 6 > • I  AC  I EXL l P 


S2t*SINIPt)*ACI3)*C3SIP3) I  HnOE 1 49 )  *X  ERLIP 

IP  IRE. EO. II  CALL  OCOEE  I  ENT *M .CONST , 0 1  ERLIP 

KO—  8061 191 *WA-R06 120 1 *48-8061  21 1  •E’lT  vv*  I  AT  I  20 1 2C APCB l*R*CONST  ERLIP 
OHUa-8061  19)*utaA-P0£<23l *0*6-806121  I *K**aT I  20) 2CAPC8  ERLIP 

OOKUa-806 1141 *00H A-AOE I  20 ) *0046- BGt (21 )*0M4  ERLIP 

OODWl>a-A0E 1 14) *OOOwA*OOC< 201*300* A-A0E(2l)*00WM  ERLIP 

TAU-IOAUI-BOEI22I *094U-PUE(7II  *044  ERLIP 

OTAUa-806 l22l*00DMO-n0EI21)*U0HM  ERLIP 

KPJ-LOAJ2-BOE(24l*0.9UWH-AOE(26l*orAU»ttOE(23)*OHli*AOEI24)*WM  ERLIP 

MITE  19.101  48,0*6,004(1,  OOOwB , wA«0«A,004 A, DUO* A  ERLIP  4 

MUA-MU*HA*4\22.*K8*4B22.  ERLIP 

04UA-04U*MA*OWA/2.*mh*0*«A/2.  ERLIP 

OOMUA«OOuu*hA*OOkA/?. *HB*30H8/2.  ERLIP 

MITE  19,111  MU.OUU.OOWU.MUA.OHUA.OOWOA  ERLIP  K 

CALL  RSLT  |X|  ERLIP 

CALL  PRINT  |XI  ERLIP 

2  CONTINUE  ERLIP 

RETURN  ERLIP 

ERLIP 

ERLIP 

3  PORMAT  I1X.9H  806  ARE  2110612.61)  ERLIP 

4  PORMAT  113612.6)  ERLIP 

6  PORMAT  1 2 OH  AMI J I  ARE  I  ERLIP 

A  PORMAT  I1X.TEII.4)  ERLIP 

7  PORMAT  l4llX,2HVIIl,2HI*£l2.5l/2ltX,2HV(||.2H|aEI2.6l)  ERLIP 

B  PORMAT  I1X.20H  ALII)  ARE  I  ERLIP 

«  PORMAT  14612.51  ERLIP 

10  PORMAT  I/IX.25H  B6IAIAI0ERIVATIVES  2 El  2. 6,  I  X, E 12. 6, 1 X.E 12. 5, 1EXLI P 

1X.C12.321X.2CH  RETAIL  I  OE Al VAT  I  VC S2 IX , 6 l 2.3,  l X , E 12. 5, 1 1 , E 12. 6 , 1 X , cEX L I P 
21?. 622 )  ERLIP 

11  PORMAT  HI, 23m  Ul  81  DERIVATIVES  2  IX,  E  1 2.6,  1 1,  E 12. 5,  II,  6  12.  52 1X.2EKL  I P 

10MUIA)  JtMVAIIvES  2 1  X  ,6 12. 6. 1 1, 6  1 2. 6  ,  U  ,  5  12.  622  )  ERLIP 

12  PORMAT  I21M  6NII.JIAV  ROWS  ARE  I  ERLIP 

END  ERLIP  - 

SUBROUTINE  6RLIP4  •••••  1 

OlMcNriO'l  JlPItll  ERLIP  2 

COMMON  OL  TEMP  ERLIP  3 
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COMMON  rEMP,EA,EAC.NUA,NUAC.EB,EBC.NUB,9UBC.MA,HBtK,r,lA,r.lB.CA,CB.EKLlP  A 
lDA,08,M,CA:»,CAPC9.NUNU,r,AMA,GAMB,K|AO>  ,  AI55>,G<  7 1  •  3  T  (  10  I ,  FN  COA  ,  FNEKL  I  P  5 
2T08»B0EI600»  tLOAG t ,AM< 50  I  ,LGAD2 , EN ( 6,  7  I  •  F NXA , FNXH, FMXA, FMX*. FNX , FME KL I P  6 
l«,E*iOA,FNOB,ACI>>  tKL IP. ELL. DEL  PHI  20  I  , V 1 6  I  .  AU  1 0 »  ,  WM ,  DWW,  DO-W ,  MB,  OwEKl  IP  7 
6B.OOMB.OOOW8.MA.OmA. MU.OWU.OOMU. OODWU •  T AU.  DT  AU. WP J« MMA.DMMA.MUA. DwEKL I P  B 


SUA 

EKLIP  9 

REAL  NUA.NUAC.NUB.NUBC.NiJNU.K.LOAOl  .L0AD2 

ERLIPIO 

OOil 1 l-A(AAI*IACI 1 1 1 ♦ A (46 1 • ACI 1) 

EKLIP11 

•0EI2>»AIAl>*(AC(  l>*«A|*A(A2)*(ACm»»2t«A(43t 

EKLIP12 

•0Em-B0Elil/B0EI2) 

EKLIPI1 

EKLIPIA 

•OEI A 1  ■ AC  441*1  AC  I  2 1 ••!! -A (A5) • AC  1 2 1 

EKLIPI5 

80E(5)-A(  At  )«(ACt  2I**AI-A(A2I*C  ACI 2  I «*2 1 « AC  All 

IKLIP16 

BOEI Al-BCE 1  A) /BOE (51 

EKL1P17 
EKLIP t 8 

BDEC 7I«ACA6I*(ACC3I**3I-ACA5I*ACI1I 

EKLIPI9 

BOEI0I-AI  At I'CACI  l)*«AI-AIA2>*CACC3t*«2>»AIAl> 

EKLIP20 

•OEC9I-BOEC  n/BOECB) 

EKLIP21 

EKLIP22 

1 

EKLIP21 

80EClOI-BOEm*CA<2ll*IACmA*2l*AC22n-AC2A»*ACm 

EKLIP2A 

•DEC tI)aAI21)*IACIll**2l 

EKLIP25 

■0EC12I-B0EI  IGl/BUEItll 

EKLIP26 

EKLIP27 

•0EC13I-B0EC6I*( A < 211 • ( AC (2 > **2 1 -A  I  22 »» - A 1 2A > • ACI 2 » 

EKLIP28 

BOEC 1  a) ■A(23)*IACIZ)**2) 

EKL IP2-# 

BOEII5I-ROEI lll/COEl  LAI 

EKLIPiO 

EKLIP3I 

B0E<16)-aDEI9)«IA(2t >•( AC  111**21-4(221 1- Af 24 >AACI  3) 

EKLIP32 

•0EI17I-AI21)*IAC(])««2I 

EKLIP31 

B0EClBI-n0Ell6l/A0E( |7| 

EKLIP1A 

P4-AC(tl*ELL 

EKLIP35 

P5-ACI2I*ELL 

EKLIP36 

P6-AC(3)*ELL 

EKLIP!/ 

EKL 1 P 18 

EN(1.1)-I4MI44)*80E( I2I-AM(45>*B0E<3) l*ACI  I l*AMC46l 

EKLIP39 

ENCt.2)>CN(l,ll 

EKL 1 P AO 

INC1»?)-(AM|AAI*B0EI 151-AMI AB)*BOE (61 ) *AC I  2 1 ♦AMI  46 » 

EKLIPAt 

CMC  l«A 1 "0.0 

EKL1PA2 

f Ml l»5l ■ I AMCAAI *80EI tt I-AHI45 I *6D£( 9 1 1 *ACI 11 ♦AMC 661 

EKL IPA 1 

ENCl. 61-0.0 

EKLIPAA 

EN<4,t>-eNll.t>*EXP<P4| 

EKLIPA5 

ENI4»2l-ENIl»2)*FXP|-P4| 

EKLIPAA 

ENt4.1)-eN( t»i)*C0SIP5l 

EKL IPA  7 

CNt4,4)-ENCl.3l*StN<P5) 

CKLIPAB 

ENI4.5l-SN<t.5l*C0S(P6l 

EKLIPAV 

EN(4.6t-ENIt,M*SIN(P4l 

EKLIP50 

• 

EKLIP51 

ENC2»1)bIAMCAMI*B0EI 12I-AMC49) *BDE(  3) t*AC( 1 1 *AM(50) 

EKLIP52 

ENI2.2l-ENI2.ll 

CKLIP51 

2NI2»ll-IAM(4Bl*8DEC I5I-AMI49 J  *BDE(  61 )• ACI 2) ♦AMI  501 

EKLIP5A 

1*12,41-0.0 

EKLIP55 

2NI2»5l-CAMt48l*80EllBI-AM(49l *0DE(9) I* ACC  1) ♦ AMI  50) 

EKLIP56 

CNI2. 61-0.0 

EKLIP57 

ENt).l)-EN<2,l)*FXP1P4| 

EKLIP5S 

ENI5»2I-£NI2.2>*EXP(-P4> 

EKLIP59 

ENI5.3)>eNt2,l)*C0S( P5I 

CKLIP60 

EN<5.4I-EN(2,3>*SIN|P5I 

EKLIPAt 

EN(5,*.»»-ENI2.5I*C0SIP6) 

EKLIP62 

EM(5,6)«EN(2,5I*SIN(  P6I 

EKLIP63 

61 


EKLIP64 

CN(3,l)*P0£(83l«ACtl l-ROEl 3l-(BuC(84l-( ACC l  >**2)*GAHB)*B0EI 12)*(BUEKL IP60 
lf|83)»(ACI 1I--2MGAMA)  EKLIP66 

INC  3, 2) *EN(  3,1)  EKLIP67 

INI  3,  31*3.0  EKUPbo 

INI  3.  4)  «OU6  I  HS»*  AC  I?  I-BOE  46  1*1  60E  (84  >•(  AC  I  2) **2 1  *GAMB  l-BOE  I 101-CBOEKL  I P69 
IE(63)-|AC(2)**2) ♦GA“A I  EKLIP70 

CNI), 01*0.0  EKLIP71 

fN(3,6)*nOE(B6)*ACm»BOEI18)*(ftDE(83)-(  ACC  3  )  P*2  >-GA*<  A  I -BOEI  9  )  •  (  BUEKLl  P  72 


lE(84)*(ACm*»2l-GA*B)  EKLIP73 

|N(6,1)*CN(3,1I-EXP(P4|  EKLIP74 

ENI6,2I*EN< J,2)«EXP(-P4)  EKUP70 

ENI6,3)*-EN(3,4) -SINIPB)  EKLIP76 

EN(6,4)-tN(  3,4I*CGS(P0)  EKLIP77 

|NI6»0}*-ENI3»6I*SI,'C(P6I  EKLIP78 

IN(6,6)*EN(3,6)*C0S( P6  I  EKLIP79 

EKUPBO 

WRITE  10.2)  IBOEin. 1-1,20)  EKLIP81W 

WRITE  (9.))  EKL1P82W 

WRITE  19,4)  MENI I.Jt.J-l.Tf, 1-1,61  EKLIP3)rf 

CALL  BINGO  IEN.V)  EKLIP84 

WRITE  (9,0)  (I, VIII, 1-1,6)  EKLIP8SW 

AMI 12)*BT(2w)/CAPCB-6DE(21l*BDE(49)  EKLIP«6 

RPl-ACI 1 l*ELL/2.  EKLIP87 

PP2  ■  AC ( 2 ) *ELL/2.  EKLIP8B 

PP3  *  AC ( 3 1 -ELL/2.  EKLIPh9 

CALL  CNECK4  <V,a:,PPI,PP2,PP3)  EKLIP90 

WRITE  (9,0)  (I, VII), 1-1,6)  EKIlP-Hrf 

EKLIP92 

ALII) --BCE I  1) -VI 1 )  EKLIP9J 

ALI2) -BOE I 3I-VI2 )  EKL IP94 

ALI3)*BUE(6I-V(3)  EKLIP90 

ALI4)»B0E(6I-VI4)  EKLIP96 

ALI0) -BO* I  9) -Viol  EKLIP97 

ALI6)—  BDE  I  91  •VI  5  1  EKLIP9B 

EKL  IP99 

ALIT)-B0EI12)*VII)  EKL I P 

AL I 8) --80E ( 1 2 ) -V ( 2 1  EKL  I  P 

ALI9)—  BOEU0I-VI  i)  EKL  IP 

ALI10)-B0E(  I0I-VI4)  EKUP 

AUll)--B0E(lb)*V(6)  EKUP 

AL(12l*BuEtl8l-V(0l  EKUP 

WRITE  <9,61  EKUP  W 

WRITE  (9,7)  (ALII ). 1-1*12)  EKLIP  M 

EKUP 

bO  1  KF -1 , l 1  EKLIP 

K-OELPH(KF)  EKLIP 

H-ACI1I-X  EKLIP 

R2-AC I2I-X  EKLIP 

R3-ACI3I-X  EKLIP 

EKLIP 

WW-Vl l)-EXP(Pl)*VI2l-EXP(-Pll-VI 3I-COS(P2I-¥I4|PSIN(P2I-VI0>-COSIPEKLIP 
13)«V(6I-SIN(P3I*80EI49)  EKLIP 

0WW*V(11*«0£(  101 1-VI2I-80EI  III  l«vm-BOE(2Ol>»VI4}PBO£(21l)»V(0)*e£KLlP 
IOE 1 221 ) * VI6I *HOE (231)  EKLIP 

DOWW-VI ll-RDbllC2l-VI2>-nOEIIl2)»V( 3)pBOEI2O2I«VI4I*BOEI212)-V(0)*EKLIP 
1B0E l2?2l*V(ol*B0EI?32)  EKLIP 

WB-AL(1I«EXPIP1I*AL(2)-EXP(-P1 )«AL( 3I-COS(P2I*AI(4|-SIN(P2)*ALI0)*EkLIP 
IC0S(P3J*ALl6l-SIN(P3l  EKUP 

OwB-ALI ll-BOEIIOI ) -AL 1 2 1 *BOE (111) -AL( 3) -ROE! 201 ) ♦ AL 1 4) -BOE 12111- AL EKLIP 
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It5l*BDE 1221 1*41161 *BDF 12 11)  ERLIP 

ODWB-All 1I®BOEI102I*AII2I*80EI1I2I*ALI3!®BDEI202I*AH4 I *BOE (212l«AEKllP 
ILI5l*BDEt222l*4ll6)*PDEI2l2l  ERLIP 

DOOWB-ALC  l  I  •dJEUOD  ♦AU?l*fiOE(ll  3l*All 3l*80f  (  203  I ♦ Al | 4 ) *BUE I  21  3  I  ♦£ KLIP 
IAl(S)*8DE  (  22  3  I  *A  L  (61  *BDE (2131  ERLIP 

Erl  i  p 

WA«AUTI*EXP<Pt)  ♦  4ll8)«£<P(-Pl)«Al(9)»C0SIP2>*AL(  10I*S  INC  P2I  ♦  All  U  EK  L  I P 
ll»C0S(P3l»Al(  121  ®SINCP»I  ERLIP 

OWA«AL  |7l*BDEIl&ll»4LIRI* 80 E(llll*AL 191*8 DE(20ll*AI  >10l*R0fcl211)»AEKLlP 
III  II ) *ODE ( 221  l*AL (  12  I  *00fc (2  311  ERLIP 

OOWA-AH 71 *m)t I102I»ALIHI«OCE(  112)*AL(9I*B0E(202I*AL( 10  I *RuE I  212 » ♦ ERL  I P 
UKU  )«BOE(  222l*4U12)*BOEt232  I  ERLIP 

OOOMA-ALI  ?1*P0EI  1031  ♦Aim*«D6l  l l 3 ) « Al I  9  I  *BOE(  2031*411  10  I  •BOE  I  2 1  3 1  ERL  I P 
14ALI  l  IMBUE  I  22  31  *  Alt 1 2  I • BOE 12 3 i I  ERLIP 

ERLIP 

FNTMW*IV( 1 1 •£  XPI PI  1-VI2I •EXPI-PII I / AC  I  1 1  *  I  V  I  3  I »COSl P2 I -V  I  4 I* S INI P2ERL I P 
in/ACI2l*IVIS»*SlNIP)l-VI6l*:0S<®ll  1/ACI  3I*BUE|49I*X  ERLIP 

IF  IRF.EO.ll  CALL  UCCFF  ( FnTwn  .CONS  f • 0 )  ERLIP 

MU«-BOEI  l9l*w4-B0EI20l*ViB-BC£l21l*FNr«<rf*IRTI20l/CAPCBI®X*CuNST  ERLIP 
OWU«-POEI  191  •UWA-B06I20I  *0we-B0fcl  21  !•*<«♦« M  201/CAPCH  ERLIP 

DOWU«-BCt  I  HI  •  OOMA-HOE  (201  •GJW8-B0H2  H*OWW  ERL  I® 

OOOMU--BOE I 19I*GOOWA-BOE  ( 20  I *000*6-806 12 1 MOD**  ERLIP 

TAU«L0A01-BUE 122  I *00hU-8 JEI 2 1 1 *L «W  ERLIP 

OTAU*-BOE (221 •OOOwu-BDEI 2  3l*0uw*  ERLIP 

WPJ«lCAG2-fli)c  (2*1  *OOUWB-eOE(25l*OTAU*BDEt23)*Owu*BDCI26l*W*  ERLIP 

CALL  RSLf  IXI  ERLIP 

CALL  PRINT  IX)  ERLIP 

IF  IX.EO.  ARSIELL/2.1  I  CALL  VINSON  |  X,  EA,  FR,  H  A,  MB,  G3  A,  C3B ,  F.NXB  )  ERL  I  P 

1  CONTINUE  ERLI® 

RETURN  ERLIP 

ERLIP 

2  FORMAT  IIX.9HB0E  ARE  /I10E12.SII  ERLIP 

3  FORMAT  I1X.23M  ENII.JI  ARE  I  ERLIP 

4  FORMAT  I1X.7C11.4)  ERLIP 

5  FORMAT  (<•(  1X.2MV  I  1 1.2H|«;  12.31/21  1X.2HVI  1 1.2HI-E12.5I  I  ERLIP 

t  FORMAT  I1X.2JH  ALII)  ARE  I  ERLIP 

7  FORMAT  IAE12.il  ERLIP 

END  ERLIP 

SUBROUTINE  EaLIPS  •  ••••  1 

DIMENSION  J1PIIII  ERLIP  2 

COMMON  OlTfcMP  ERLIP  3 


COMMON  TEMP,EA,EAC,NUA,NUAC,E8,=BC.NUB,NURC.MA,MB,R#C3A,G30,CA,CP.F'aiP  * 
lOA.OBtM.CACB.CAPCB.'IUNU.GAMA.GAMB.M  40  I  .  A (  3i  I .  Gl  7  I . I>  1 1  3C  I ,  AN  T CA . F.4tRL  I  P  •» 
2T0B,  BOE  16001  .LOAD  l  .AM  (So  I  .1  „A02.f  N(  6.  7  I  ,  FNX  A  ,  FNXB  .  FfX  4,  F  M  *0.  F  N*  ,  FMCM.  I  P  6 
JX.FNOA.FNCB.AC  «  3  I  «RL  I  P.  E  l  L  .DELPHI  20  I  ,  VI6  I  .  »L  (  t  rt  I  .  MV,  OtfW .  C  Da*  ,  *6,  L  I  P  7 
AB.OOhR.OOOMB, MA.CMA.hU. O4U.GJWU.00CMU, T XU.  l)T AJ,*P J. AtfA.OWMA,  MUA, UhEaL IP  O 


5UA  ERLIP  9 

REAC  NUA,N(IAC,M)B,NUBC,NU*AI,R,l;aO1,L'0AO2  erlipio 

BOE 1 1I«\I44|«|AC 1 1 l**3l*AI4»»*ACI  II  ERL  I  PI  1 

60Et2l>A(MI*<AC(l>**<>l*A(42l*(AC(l>**2l*A(4l|  ERL  IP  t  2 

B0EI3I-H0EI1I/B0EI2I  ERLIP1J 

ERL IPI4 

6DE <4  I *4(44 1*1  AC  I:ia«)|rA(43I«ACI2I  ERL  I  PI  3 

BOE (51*4(41 » • ( AC <2I**4|»a(42I*(ACI2)**2I*A(4II  ERL (p 16 

BOE  l6l»B0E(4l/ei)tl5l  ERL  I  ®  l  7 

ERL  (PH 

bdei7i«a(44i*ia;i  h*»ii-a(4si*aci3i  erl  i p i 9 

BOE  IBI-AHIIM  AC  ())**<•  l-A(<,2l*(XCm**2l*AI431  Erl  11*20 

eOEl9l*BDEI  71/BOE  (81  ERHP2I 


ERL IP22 
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»DE(  tOI»BOECll*C  AC  211*1  ACC  1  >»*2>*AC  221 l-AC 24|»ACC 1 > 
BOE Cill«AI2SI*CAC<ll**2l 
BOE  I 12l»5CEt 10I/60EI 11 ) 

BDEin)-nUE(6)«IAI2n*(AC<2l**2)«AI22l)-*(2«)*AC(2l 

BO£Cl4)«AI2SI*IAC(2)**2) 

BOtCl»l-60EC13)/eOE<  >41 

•OE  CI6l»BOEI9)*CAI2l)*IACl3>**23-AI22l)-A.I24l*ACC31 
BCEI1 7>«AI2J)*IACI3)**2I 
BDEI I « > -noe I  16 > /POE  (  17) 

P4-ACI1I*ELL 
B5* AC  1 21 *ELL 
P»>ACI)I*ELL 

EMC1*1>*IAM|44|«B0E<  12l-Am43)*BOEI  3)  »  *AC  ( l )  ♦  AM)  46 ) 
EMCl,2)«cNll.l> 

ENI 1,3)-UWI44)*»0EC  181-AMUS)  *BDEI 6 )  ) *4CI 2  I  ♦  AMI  46 ) 
ENI1.4)<EN(l.3> 

EMI  l,5)-<  AH|  A<.I*6DEI 18I-ANI45I »BDE(9) )*AC(3)*AM(46V 
EMCl,6)»0.0 

EN  ( 4 ,  1)«£NI1*1)*E*P|P4) 

EMC4*2)*EN(  l,2l*E*PI-P4l 
EMC4,S)«E*Cl.3>*EXPIP3> 

EMC4,4)«ENI1 ,*>*ExPt-P3> 

EM C4»3I«ENI l*3l*C0SIP6> 

EMC  4*61 -EN I I *8) *S I Nl  P6 1 

EMC2, II *14*148 )*BOCI 12I-A«I49>*B0EI II  I  •ACID  *AHI  601 
EMC2*?>-tN(2.1) 

EMC2,3I*IAMI*«)*P0EI  15)-A*U9I*B0EI6I  I  *AC I  2 1  ♦  AH(  50  I 
EMC  2* 4 ) •[ N I  2  • 11 

£NI2,8»«IAMl4d)*P0EI  m-AM|49l*B0EUI I • AC ( 3  I ♦ AMI 50  I 
ENC2. 6)0.0 

EM45, 1)«ENC2.1I*C*PCP4I 
EMC5,2I-ENC2.2I*‘KPI-P4| 

EMI8*3)SENI2* J)*EXP(P5) 

EMI5»4)*ENI2*4l*t XPI-P5I 

EMCS,8)-E;*I2.8I*C0SIP6I 

EMCS»6)-ENI2.0)*SINIP6) 


EKLIP23 

EKLIP24 

EKLIP28 

EKLIP26 

EKLIP27 

EKL  I  P?B 

EKLIP29 

EKL I P  30 

EKLIP31 

EKLIP32 

EKLIP33 

EKLIP34 

EKLIP3S 

EKLIP36 

EKLIP37 

EKLIP3B 

EKLIP39 

EKLIPO 

EKLIP4I 

EKL  IP42 

EKLIP4J 

EKLIP44 

EKLIP43 

EKLIP46 

EKLIP4  t 

EKLIP48 

EKLIP49 

EKLIP50 

EKLIP51 

EKLIP52 

EKLIPS3 

EKL I P54 

EKLIP88 

EKLIP56 

EKLIP37 

EKLIP5B 

EKL I P89 

EKLIP6G 

EKLIPbl 

EKL1P62 

EKLIP63 


EKLIP64 

EMC3.1l-00ElB8>MCtl  I-B0E(3I*CB0EIB4|P|AC(II*P2MCA«BI»B1)ECI2)*CB0EKUP65 
SEC  B3)*l  AC  I  I  »**2I  ♦‘.#AMAI  EKL  I P66 

EMC  1» 2) *lNI 1*1)  EKL1P67 

EMC  3, 3) •"OEI Ail* ACC  2  I -HOE  16 1 • C BOE 1 84 ) •< ACI 2  I *»2 1 ♦ CAMB I ♦BOE 1 1 5 !• I BOEKL  I P6B 
1ECBJ)«(ACI2I*«2>  OAMA)  EKLIP69 

EMC  3,41— £•*!  4,31  EKLIP70 

EMC  I*  5 1 »0. 0  EKLIP71 

|MC3.6)«ROEIB,»)*ACm*BOE(lBI*CBOEIB3l*(ACC3)P*2l-CAPAI-OOEC4)*l80ENLIP72 
|ECB4|*U;<3)*«2)-r,AKB)  EKLIP73 

EMC6,l)«ENC3»ll*t<PIP4|  EKLIP74 

CMC  6. 2 1 •TNI  3,2  l*f <P|-P4|  EKL I P  78 

£NC6|3l«ENCl*3l*c<P(P4|  EKLIP/6 

EM  C  6, 4 ) ■*  Ml  1*4  I *fc  API  -P9 1  EKL  IP  7  7 

ENC6,8)«*E*I|  )»6l*SI*l|PM  IKLIP78 

EMC 6, 6) »f HI  3*61* .0SIP6T  EKL IP 79 

EKLIPBO 

WRITE  19,2)  IBOEI  I  I. 1-1,20)  EKLIPB1M 

WRITE  19*3)  EKLIP82W 


I 


MUTE  (9,41  l(EN(|,J)«J*l«7l, 1*1*61 

IKUPBIw 

CALL  BINGO  IEN.VI 

iKLIPb*. 

MUTE  19,  SI  II, VIII, 1*1. 61 

EKL1PBSW 

AAI12I*8TI20I/CAPCB-BDEI21I*BDE(491 

EKl 1 PH6 

AA1 *AC 1 1  !*Ell/2. 

EKLIP87 

AA2*AC(2l*Ell/2. 

IKlIPBb 

AA3-ACI )>*ELl/2. 

CKL1P89 

CALI  CHECK**  IV,AC,PPl.PP7,PP)l 

EKLIP90 

MUTE  19, SI  II,  VIII*  1*1.  ttl 

EKLIP9IW 

CKLIP92 

ALI 1 B «— BOE  t ))*V( 1 1 

EKLIP9) 

AL  I  2 1 ■ S  0  e I  3  1  *  V  I  2  I 

EKL 1 P94 

ALI  )»  «-P  IE  161 *VI )l 

EMIP9S 

AL  «4I*BI)EI6I*V(4| 

EK11P96 

AL(SI*BOE  (9l*V(M 

EKLIP97 

ALI6I  »-9r>E  1  91  *V(  S  1 

EKUP9C 

6 

EKl I P99 

ALI 7l»BOE(12l*V(ll 

EKIIP 

AL 1 8 1 *-POE I 12I*V(  21 

EKLIP 

ALI 9) «BOC 1 lS»*V! )l 

EKLIP 

ALI  1 0 1  *  -  M.  0  E  <  ISI*V(4I 

EKLIP 

AL llll*-30Ellol*VI6» 

EKIIP 

AL 1 1 2 1  *  B 0 E 1  l  9 1 *V ( S 1 

EKLIP 

MUTE  19,61 

EKLIP  M 

MUTE  (9,71  (All  1  I  ,1-1,121 

EKLIP  - 
EKIIP 

00  1  KF-1,11 

EKLIP 

K—DE l PM ( Kc I 

EKLIP 

A1*AC I  1 1  *  X 

EKLIP 

a2*a:i2i*x 

EKLIP 

A3* AC  I  31 *X 

EKLIP 

EKLIP 

MW* VII 1 •  =  XP I  PI I-VI2I  *EXP(-P|  |*V! )I*EXP(P2I*V(4)*EXP< 

-P2I»VI51*C0S(EMIP 

t*j|*v(6»*$ii«n»i»MD?»*'>i  Em.  i*> 

ovw«vm*eot  i  ion  *vi  ?>*boei  m  i»vm*90Ei  t6t  t«vt4t*BCct  m  i«vt sipbeklip 
IDE <221 >«V(6l*jnE (2)1  I  £KIIP 

OOWW* VI  I  I  *M)b  t 152 I ♦ V ( 7  I *B OE  I  1 1  2  I ♦ V ( ) 1 *BOE ( 1 62  I *V< 4  I *BDE<  172I4VCSI4EMIP 
IBOE (222I*V(6)»,)DE 12121  EMlP 

MS* All  ll*EXPIPll  MUJIU»PI-Pl  »  •  At. «  3I*EXP|P2  I  ♦  AL  I  4 )  *EKP  I -P2 1  ♦  AL  I  V  I  EM  I  P 
t*COSIP)l  »ALI6I«SI’HP)I  EKIIP 

OMft*All l  l*M£  I  t0tl**LI2l*SDE<lll l*AL<  ) I *0051  161 »  **1 «  *  1  *BDE  I  1  7 1 1  ♦  Al  EM  1  P 
l<SI«B0E(72ll»ALI6l*«UE(21ll  EKIIP 

DOWS-Al  (  1  »*PJE<  l  ,2  !♦  All  2  I  «BCE  ( 1 12  I ♦ AL  I  )  l  *60E  1 1  62  1  *61 C  6  I •BDEI  172  I  *AEMIP 
lL(5l*PDE(222l«Ai  IM*R0t(2)2l  EKLIP 

OOOwn-tU  1  IMUEI  1C)I»AI(2I*90EI  11)1* All  »  »  ••SOf «  16)  I  ♦  AM  4 1  •  HOC  C  ITJI^KLIP 
1AI(S)*B0E(22)I«41  (6)  *BDE  (2)31  EKLIP 

EKIIP 

MA-AL  I7»»EXP(PII«AL<AI*EXP<-P1I*AL<9)*EXP(*2I*AI(10I*EXP(-P2>*AI(IEMIP 
I1)*C0S(P)I»4U 12 1 • SI  N  (  P3  I  EKLIP 

DMA*  AH  7I«0j£ ( 1011*4119) -BOEI  l 111 « Al(9l«B0E<  161 1  ♦All  10 l* 60c  I  l  711  ♦  AEM  IP 
1L ( II l*BOE( 221  MALI  17 )*ROE (2)1 1  EKLIP 

OOMA'Al  I  7l*b)E  (lC2l*AUoi  •3UEI  1 12  I  ♦  Al  I  9 1  *B0E  1 1621  ♦  ALI10 1  *AuE  I  l  72  I  aEM  I  P 
1ALI ll)*60E( 2221* M(12I*BU5 (2)2 1  EKIIP 

DOOWA  * Al  (  7  1  •AoEl  l  ))1 -Al (  1 1  •■1UE  (  1 l )» • All  9  I *PDEI  16)  l  ♦  All  10  » *60E  Cl  7  )  I  EM  I  P 
I  PAL  ftl)*AD£(22)!*All  12XCOEI2)))  EM  IP 

EKIIP 


FNTWW-I VI  n*fXP|Pll-VI2l«EXP(-Pm/AC(l)t|VI)l*EXP(P2)-V(4l*EXPI-PEKllP 
121  l/ACI 2>»l  V(  )I*$I’J(  P)I-VI6I  *CGSl  P)  II /AC  I  )  I  ♦  *UEI  <9  I  *X  EM  IP 

If  IX.EU. -Ell/2.  I  CAtl  UC5FF  I FNTrfW.C ONST I  EKIIP 

MW»-BOEI  191  »MA-B0EI20I  •HH-BOE  I  2 1 1 *F NT 4W» | BT |  201/C APC61*X*CWN$T  EKIIP 
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o*u»*bdeii9i*dwa-boei2oi*o*b-boei2i»*w»m6T(2oi/caacb  cku  a 

OOWU«-AOE( 191 *DOwA-BDE  1 20 1 *ODMB-BOE  <  21 l*OMV  CILIA 

OOOMU**BOE  1 191 *OOOWA*BOE<  20 I*003WB-B0£l  2 t 1 *OOMM  CILIA 

TAU'LOAOl'BuE 122 ) *00mU-SOE 1 2 ) t •UHM  CILIA 

OTAU— BOEI 22  >  *OOOMl>-BOE  1 2 )l *00W*  CILIA 

VAJ*L0AD2-B0E 124 I •  00DWB-6DE 1 25  I  *DT AU'BDE 12)1  A0WUABDCI 261 *WM  CILIA 

CALL  ASLf  IX)  CILIA 

CALL  ARlNT  III  CILIA 

OA* I  HA SI 2, I ATAU»C4MB*(MA»0MMI  CILIA 

OB*  I HB  / 1 2> I  * t  AlHG A*rf* I mB*OMM|  CILIA 

FQ*04«QS  CILIA 

MITE  19,11  0* ,08 , AO  CILIA  w 

IA  (X.CQ.ABSIELL/?. I I  CALL  VINSON  <X,€A,£S,HA,MB,G)A,GSB, ANXB  )  CILIA 
I  CONTINUE  CILIA 

IlC  TURN  CILIA 

CILIA 

a  AORMAT  III, 'HIDE  ARC  /(10E12.SII  CILIA 

3  AORMAT  I1X.20H  CN||,JI  ARE  .  I  CILIA 

4  AORMAT  111, IE l l« A)  CILIA 

5  FORMAT  IAIIX,2HV||l,2H|«EI2,SI/2IU,2HVI  1 1, 2H)*C  12.SI )  CILIA 

4  AORMAT  I ll«20H  ALIII  ARE  I  CILIA 

T  AORMAT  I6EI2.SI  CILIA 

0  AORMAT  lll,4HgA*ei2.9,?l,SHG0*C12.3,21,SHQBAR*C12.S/l  CILIA 

CND  CILIA  - 

SUBROUTINE  CILIA4  •••*•  I 

OIMENSION  JIAIUI  CILIA  2 

COMMON  OLTCNA  CILIA  } 


COMMON  TCMA,CA,EAC,NUA,NUAC,CB,CBC,NUB,NUBC,HA,HB,R,G3A,G)B,CA,CB,El  t»  * 
IOA,OB,H,CACB,CAAC8 ,NUNU, GANA.C AMStll 4u),Al)S)«Ct  7)«9Tt St  I ,f N T CA, FN£*Ll A  5 
aTOB.BOCIAOOl.LOAOl.AMISO) ,L0A02,CN(«, ll, ANXA.CNXB.AMXA, AMIS, ANl.FME ALIA  6 
1I,ANCA,FN0B,«CI3> ,ILIA,CLL,OCLAH(20I,VI6I,ALI  tll,MV,OMN,UOMM,WB,OMEKLIA  I 
4B, OOMB, DOOMS, HA, DMA, VU.OVU, OOMU, OOOtfU, TAU, OT AU, WA J , MM A , DMVA, WUA, 0«E ALIA  i 


BOA  CILIA  9 

BEAL  NUA.NUAC ,NUB ,NUSC,NUNU,I, LOAOI.LOA02  CILIA10 

CASE  OF  THREE  REAL  MOOTS .UNEOUAl , AOS  I T I VC  CILIAI1 

BOC I ll*AC44l*|AC I  I  I I • A 1431 • AC I II  CILIA12 

B0EI2I*A«M)«IACI  l)**4|»A<42)*IACUI**?>«AI4)l  ciliais 

BOCISI-BOCID/BOCI2I  CILIA1A 

CILIAIS 

B0CI4I*A(44|A|AC I  2I**))*AI4S)*ACI2)  CILIA  16 

BOC  I  SI *A I A| )•(  AC  I  2 1**4) * A 14?) •  I  AC  121**2  I *A(4)|  CILIA) I 

B0CI6 l«B0CI4l /IOE I  SI  CILIAIS 

CILIA1V 

B0EII>*AI44|6|ACIS)**SI*A|4S)*ACISI  CILIA20 

BOEI t l*AI41 I*IACI)I**4|*A142I*IACI)I**2I*AI43I  CILIA21 

BDCI 9>*B0E I II /BOE I <1  CILIA22 

CILIA?) 

BOEI 101 *BOEI )l*l  A  I ?1 )• ( AC  1 1 1**21  * A(22 1 1 ~A(24  I *AC 1 1 1  CILIA24 

BOC  I ll)*AI2S)*IACI 114*21  CILIAIS 

B0CI12I*B0CI IOI/BOCI 111  CILIA?* 

CILIA?? 

BOC  I ISI*B0C(*I*IA(?1 1*1  AC  121**2) *A 1 22 1 1 "Al 24 1 A AC I ?)  CILIA2* 

BOCII4l*AI2)|A|ACI?IA«2l  CILIA29 

BOCI1SI*BOCI ISI/BUCI 141  C1LIAS0 

CILIAS1 

B0CllAI*BDCI4l*l  AI21l*f ACISI**2I*AI22II"AI?4I*ACISI  CILIA32 

BDCI I 7l*AI2il*< AC  I  31 ••?)  CILIA)) 

BOCIlB)*BOEIl*l/BOCim  HLIAS4 

CILIA)) 


CILIA) 
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P4»AC(l>*Ell 

EK11P37 

P5«AC(2l*ELl 

EKllP3b 

P4>AC(3)*ELl 

EKUP39 

EKIIP40 

EN(UI)-<AM(44l»n0E(12)-AM(4S)*B0E(  >1  )*AC(  1  >  ♦AMI  46) 

EKL1P41 

ENtlt2)*ENU,l) 

EKLIP42 

EN( l,3)-<  AM(44I«B0E(  m-AN(4SI*B0EU>>*AC(2>»AHt46) 

EKLIP4I 

INI 1*4)>ENI l»  31 

EKLIP44 

EN( l ,5)»(AM(44»*B0EI l6)-AM(45)»B0Et9) )*AC(3I«AM(44I 

EKLIP43 

EN(l,6»-EN(l,5) 

EKLIP46 

EKIIP47 

ENI4,l)*EN(l,l)*tXP<P4) 

EKUP4B 

EN(4.2t-EN(l«2l«EXPt-P4| 

EKIIP49 

EN(4,1)«EN(l,)l*EXP(PSt 

EKLIP30 

EN(4.4l«ENt l.4|»CxP{-P5) 

EKLIP51 

ENI4*S)-ENI 1.5)*CXP(P6I 

EKUP32 

EN<4,6I«ENU,6)*EXP(-P61 

EKIIP53 
Ext <P34 

ENI2,ll«UH(4b)*DOE(12l-AM(49|40DEmi*AC(lUAMI50) 

EKUPbi 

EN( 2«  21 *EN( 2 , 1 1 

EKLIP36 

EN(?,3)-< AM(4ul*B0E< lBI-AH(49l*8DEt6l >*4CI2I«AM()0) 

EKLIPS7 

EN(2*4)-EN(2,3I 

EKLIPSB 

ENI2,3»-«AH<4a)*rOEI  lB)-AM|49)»nDE(9) )«AC(3I»AM(30I 

EKIIP59 

EN«2t6)«EN(2,M 

ERLIPbb 

EKLIP61 

EM(5il)»EMf2»ll*EXP|P4| 

EKl  1 P62 

EN(9,2)«EN<2,2>*EXP<-P4) 

IKIIP43 

CN(4«3)«EN(2'i>*EXPtPS) 

EKLIP44 

£N( 5( 4) •EN(2»4I*EXPJ-P5I 

EKLIP6S 

ENI5.5)*EN<2,4I«EXP(P6I 

EKL IP66 

EMS.6I*EN(2*6I«EXP(-P6I 

EKLIP67 

EKLIPbb 

EN(3.1)«P0EI33l«ACm-e0E(3l*(B0E(a4)«l  AC( 1>*'2I ♦GAMB?»BOE(  1  2  »•  <  60EKU  P69 

IE  « S31  •  (AC  I  U**2)«UA*A) 

EKIIP7G 

EN(3.2)*EN(3,1I 

EKUP71 

EN(3t3»*tOE(HD»*ACI2 t-BOE <6 )• < 6UE 1 841 »( AC(  2 > **2 > «GAMB 1 »B0E  1 131*1 BOEKl IP  72 

1E(83)*(AC(2>**2I ♦GAMA! 

EKIIP73 

ENI3.4)—  EN(  § •  3> 

EKIIP74 

EN(3.5)*ROEl*3>*AC<3)-6DE(9)*<BOE<84}*(ACm**2}«CAHOI*BOEi  16l*(nOEKi:P7S 

IE  (  •  31  >  •  1  AC  I  3»**2I*GA*A1 

EKl 1 P  76 

ENC3,6»  — ENU.51 

EKUP77 

EKIIP76 

INI 4il)*ES(3fl)*EKP(P4| 

EKIIP79 

lNIA.2l-ENII.2l«EXP|-P4l 

IKIIPHO 

EN(4.3I«ENI j, J|*EXP(P4I 

EKLtPBi 

IN  1  4  •  4  )  >E  N  (  )t4in  XPI-P3I 

EKIIPB2 

CNI4,9»«E9< 1,)I*EXP(P6I 

EKlIPd) 

EN(4«4)*E9( 4i4I*EXPI~P4I 

EKl 1 P  64 
EKIIP63 
EKl tPub 

AK1I— BDEf  >l*Vt  |» 

EKl  (**87 

Alt 2) *B9E  (3t*V(2) 

IKIIPBB 

EKIIPB) 

Al(3l— 60E(6I*V<  11 

EKIIP90 

Al ( 4 1 ■ BDE (A1*VI4I 

2KIIP91 
EKl IP92 

Al(SI»C0EI9l*V(5l 

EKIIP93 

AU6I»60E<9>*V<6) 

EKIIP9* 
EKl lP9b 
EKIIP96 

ALI7)-BDCI12I*VI  I)  KLIMT 

ALUI»B0EI12)*VI2)  EKLIP98 

EKLIP99 

ALI9|*B0Ef 1SI*VI  3)  KLIP 

ALIIOI  —  B0EI1»)PVI4|  EKLIP 

KLIP 

ALf 1 1 1 "IDE  II #1 *9151  ‘  EKLIP 

ALIl2)»BDEtl8l»V|6|  EKLIP 

CKLIP 

Ntl»miAEXPIPllAVI2l*£XPI-Pi;AVm9£XPIP2)*mi4EXPI-P2lAV.5>4EXPIEKLtP 
1P1)»V(4I«EXPI-P9>«60EI49)  EKLIP 

DWW*VI1I*60EI101  MV(2)«Bt)Ellll>«Vt3)»60Etl41IMt4)*B0Et  171  )«m>*BEKLlP 
10EI 24 1 1  *  V ( 61 *BOE I  2b 1 1  EKLIP 

DOtfN-Vtl)«HJE(192)«VI2l*BOEI112l»VI  3I*BOEI162|aVI4)*BDEI  172)»VI5)*EKLIP 
18DE 1 242) ♦ V (6) *BOE (252)  EKLIP 

WB*AH1)*EXPIP1>*ALI2)*EXPI-P1)*AH1)*EXP|P21*ALI4)*EXP<-P2>»ALC5)EKLIP 
1*EXP|P3)»AL<6I»EXP|-P3)  EKLIP 

OMB* AL( I )*ROE ( 191 ) «AL 1 2 ) *BDEt 1 1 1  )  +  ALt 3 1 *6DEI  141 UAH 4) *BDE1 1 71 1 ♦ ALEKL 1 P 
ll*UftOEI741MALt6)«*OE(291)  EKLIP 

OOMS- AL  ( 1 )  «BJE  1 1 02  >  ♦  AH  2 )  *80£  1 1 12  )  ♦  AL  111  *BOE  ( 142 1  *AL  1 4 1 48DEI 1 72 1  ♦AEKL  I P 
lLm*B0E(S0fc(242l«4lt6l*80E<252)  EKLIP 

ODDWB-ALt  1U6UEI 1 031 ♦ALt 2 l*BOE 1 1 13 > Mil 3 UBOE C 14 J) AAL 1 4) MoE 1 1 7J) *EKL IP 
1ALIS)*BDE<241UALI6UB0EI253)  EKLIP 

EKLIP 

NA*ALITUEXP(P1MALI8UEXPI-P1MAH9U£XP|P2UALI10MEXPI-P2UALIIEKLIP 
UUEXPIP3UALI  12UEXPI-P3I  EKLIP 

OWA-AL17UH0EI1OI  UAH8MB0SI 1 11  UAL!  9 )  «8DEI  141.  «AL  I  10UBDEI 171MAEKLIP 
1 L  <  1 1 )  *8DE  {241H4Ltl2)*B9E(251)  EKLIP 

00WA-ALI7MBDEI1O2UAH8H8DEI  U2U4LI9UBDEI142UALI  lOUBuEl  172UEKUP 
1ALIU  UBUEI2*2UAL(12)  •B9EI262  I  EKLIP 

ODOtfAaAL  (  7) *BUE  I  103)  ♦AttA  I  *8uE I  113)*  AL  (9)  *8DE  (  163UAL 1 10  UBOE  1 173  )EKL  IP 
UAL ' 1 1  UBOE (2*3) ♦  AL  ( 17UPUEI253)  EKLIP 

FN’MM«|VI  lUsXP|Pll-V(2)«EXP|-Pl))/ACI  1  UIVI  3MEXPIP2I-VI4UEXPI-PEKLIP 


12)1/ACI2UIVI5UE<P|P3)-VJ6UEXP|-P6)  l/ACI  3UB0EI  49MX 
IP  IX. EO. -ELL/2. I  CALL  UCOFf  I FNTWN .CONS T I 

MU*-BOE  ll9l*WA-8OEC2O)*WB-8OSI21)*F‘4Twy4(0n2OI  /CAPCB  UA+CQNST 
OWU  —  RDEI  INU0WA-RDM20)  •0*B-B0EI21  UWW*ttM20)/CAPCB 
00MU--80EI  IVU00*A-R0FI23UD0WB-A0EI21  UO*W 
OOOt/U-  -BOE  1 19 1  ♦OOnwA-RQE  I  20  )  »OOCriB-BOE  I  2  1 )  •OOWM 
IAU«LOAO)-BOE(22)*OOhu-BOEI23)*Ohm 
OTAU»-Bf)EI22)*OOOWO-BOEI2  1U0UWW 

HPJ-LCA02-B0EI24 )  •900W0-B0EI2S  )  *OTAU«BOEI  23U0NU*DDEI  26UWtf 

CALL  tESULT  IX) 

CALL  MINT  IX) 

ketukn 

END 

SUBROUTINE  OIPP  IP1.P2.P3I 


EKLIP 

EKLIP 

EKLIP 

EKLIP 

EKLIP 

EKLIP 

EKLIP 

EKLIP 

EKLIP 

EKLIP 

EKLIP 

EKLIP 

EKLIP 

EKLIP 


COMMON  OLTEMP 

COMMON  IEMP.EA.CAC.NUA.NUAC.EB.EBC.NU6. NU8C.HA.HBtR.03A.C3B.CA.CB. 
10A.UM,H,CAC6.CAPCP*NUNU.CAMA.CAMn,K(AOI,AI4S).CI 7I.1TI 3CI.FNTCA.PN 
2T08.90EI 6901 .LOAOl , AM  I  40 ) .L DAU2 . ENI 0, 7 ) , PNX A, PNX8 , FMX A, FMXB. PNX. FM 
3X.FN0A .PNOB. AC  I  3  I .KL IP . ELL.OELPHI 20 ) . VI 4 ) , AL  I  1 8 1 . MM , OWW , 00«H .  ,8 . On 
46.00MB,0UOMBiMA.OWA.MU«OHU.UUMU(000«U.T AU.OTAU.hPJ.hpA.OmMA.IiUA.Om 
5UA 

REAL  NUA, NUAC, NUB, MUBC.NUNU.K. LOAOl. L0A02 
•OE II01)*ACI1)*EXP|PII 
•0EI1C2I-M0EI101 )«ACI1I 
•DEI  1031 *8011102) •ACIll 
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0DEI104) *60E  1 103 1  *  AC  111 


BOEum—  ac<ii*exp<-pi» 
BOE(im«-ftUEIlll>*AC(ll 
•0Eimi*-euc  1 1 1 2 1  *aci  it 
BOEI114I«-60E(U3I*AC<1) 


BOE(121)'EXPIP2I*(AC<2)«COS(P9  l-ACI  3 1 *S 1N(P3I I 
BOEI122 I *E XP«P2> • t A« (1) *C0S(P3 1-ANl 2  )*Sl N(P1 1 > 
BOE  ll23)*EXHlP2l*iAM(3l*C0S(P3  |-AH(<*I*$INIP3H 
BOE (124)  »EXP(P2)  •( Am(5)*C0S«P3  l-AB< 6)*S(N(P3  J  ) 


BOE (131 >«EXPIP2)«( AC(3I«C0S<P1 1 *AC ( 2 1 *S I N| P3 > I 
bde(13?)*e*p(P2>*(4<m?)*cosipii»4pi  1)*SInimi» 
BOE(lS3)«eXt»(P2>*<  A4(4I*C0S(P3 MAPI  ))*$ IN(P1) I 
BDEI134MEXPIP2)  •  <  AM6MC0SIP3  )  ♦  4* (  >  I  •  $ I  N(  P3  I  » 


BDEI141)— EXPI>P2)*IACI2I«C0S(  P3MAC13)«SlN(P3n 
B06(l42l«»£XPI-P2l*U*l  l  »«C3SI  PJI*4.M<  ?l«>niP3)l 
BOE(l43)*-ExP(-»2»»l API 3l«C0$IP3M4>>|4)«SINt  P3I) 
BOE 1 1441 ■♦£  XP I -P2 ) •( API  5 1 *CvS ( P3I ♦AM(4l*$ INI P3 1 1 


BOE(1511»*EXP|-P2l*( AC(3IAC0S(P3I-ACI2>*SIN|P3>I 
BDEI152I  —  SXPC-»2)*IAP(2I*C0S(PI»-AP|  ll*Sl*J»  “» *  * » 
BOE  (lS3)**E4P(-P2)*IAPI4|*CJSI^3l-API3l*Sl*MP3ll 
BDEI154)»*EXPI-P2)*( AP(6> *C  OS  « PJ ) -AH ( 5  I • S I N( P 3 >» 


BOE(l6l)«ACm*EXPI»2) 
BOE(162l*4C(2l*OOc(161l 
BOE I 163)*4CI2I*90E(162I 
B0E(^64)«4C(2)*P0EU63I 

BOE( 1 71 l»-AC(2)*E*P(-P2) 
BDEIl72l‘-AC(2)*P0Ell7n 
80EI173)— AC(2)*4l>E(172) 
BOE( 17*1 »'4C(2I*B0E( 1731 


'  BOE (1611s- AC  I l  l*S I N( PI ) 

B0E11B2I—  (4Cin**2»*C0S(  PI  1 
BOE(18  3l-<ACm**3)«SIN(Pl) 
BDEI1S4I-IACI  M**41*C0S(P1) 


BOE(I9II*ACI1I«COS(PII 
BDC<192M-IAC<1>**2>»SINIP1I 
BOE (I93l«-IAC( ll**3l*C3S(Pll 
BOE < 1941*1  At  I 1)»*4I«SINIP1) 


BOEI20D—  AC(2I«SIN(P2I 
BOE 1202 » *- 1  AC (2I**2I*C0S(P2) 
BOEI203) «( ACI 2)** J  )*SlN( P2) 
BOE(20A,«(AC(2)**A>*C0S|P2) 


14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 
27 
20 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 


69 


74 

75 

ROE  1 2 1 1 ) ■ AC  I  2 > *COS (P2 )  76 

BOEI2l2)a-(ACI2l  •Si'll  1*2 »  11 

BO£(2l3>a-(ACI2>»*)l«C0S<P2>  78 

BOE l214)«(AC(2)**4)aSIN(P2)  T« 

10 

•1 

BOEI221)— ACm*5IN(P3)  B2 

BDEI222la-IAC(3M*2»*C0SIP)l  83 

BDE<223l«IACm**>l«StNtPll  84 

80EI224|a|ACI3)**4l*C0S<P3)  85 

86 

87 

BDE(7)U»ACm«COS<l»3l  18 

RDEI232I*- (AC(3)*«2>*SINIP3I  89 

80E(23))--IAC(3I «*1)»C0S(P3I  90 

BOE  ( 2  34)a|AC(3)**4)*SINIPU  *  9i 

92 

93 

•OE(241)»ACm*£XP(P3l  94 

BOE ( 242  taAC(3)*ftOE(241)  93 

B0E1243I>ACI3»*90EI242)  96 

BOE( 244 ) «AC I  3 ) *6UE (243 1  97 

98 

99 

BOE(251)»-AC(3)«EXP<-P)l  100 

ROE(252)a-A:<  JI*flOEI25U  101 

8DEI2S3I— ACm*00EI252l  102 

ODE I  254)*- AC I 3  3*  BOE (2591  103 

104 

RETURN  105 

ENO  106- 

SUBROUTINE  POLYA  IN.COEFF, ROOTS, 01  •••••  1 

DIMENSION  AIM, 31,  IA(51,3),  R00TSI2.N),  0(1),  C0EFFI1I  2 

INTEGER  DEGREE  3 

OECREE-N  4 

NlaUEGREE ♦ 1  5 

R-10  6 

MMAX-15  7 

OELTA-0.0001  8 

IPS 1 LON ■ 0, 000001  9 

00  I  I • 1 • N1  10 

All, 1 laCGEFF 1 1 )  11 

|AII,l)aO  12 

CALL  SCALE  IAII,l),IA||,ll)  13 

1  CONTINUE  14 

CALL  RSSR  ( A, I  A. ROOTS, OEGREE, A, UMAX, DELTA, EPSILON, 0)  15 

IF  (Nl-IUEGMEEM  I  )  3,3,2  1* 

2  RETURN  17 

2  PRINT  4  ISP 

RETURN  19 

20 

A  FORMAT  ( 21HOSOME  ROOTS  NOT  FOUNO)  21 

I NO  22- 

SUBROUTINE  RSSR  M,  IA,  ROOTS,  DEGREE, R,MMAX,OELTA,FPSILON,  01  •••••  I 

DIMENSION  A  1 5 1 , 3  )  ,  I  A  I  5 1 , 33 ,  R0CTSI2.50I,  0(51),  RONG0(5O),  MROMGO  2 

1150),  N0NRTI50),  NNONRT 150)  3 

INTEGER  DECREE  4 

NaDCGREE  » 
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IF  (HI  1,1.1  * 

1  DCCREE-NCUR  T 

2  RETURN  * 

1  HI*NM  9 

M2«Nl*l  10 

DO  S  I-l.N  11 

A-N2-1  12 

IF  (A(K, 1 1 1  6,9, 6  11 

4  J-Nl- 1  19 

ROOTS! 1 • J) a0. 0  IS 

ROOTS  ( 2  •  J I  B0*  0  U 

$  CONTINUE  17 

0E69EE-O  It 

CO  TO  2  19 

4  N1«K  20 

N-K-l  21 

NCUft-N  22 

Ml-N  21 


7  CALL  A00TSQ  I A , t A , NCUR*H I  29 

CALL  REALRT  I  A, I  A ,M,NCUR , OELT A ,£RS1 LON,ROMOD,HROHOD,NONRT, HNONRT ,N  25 

ICO, ROOTS  I  2* 


IF  (NCOI  12,12,8 

27 

• 

WI-NCUAM 

28 

CALL  CONPRT  ( A,  I  A , ROKOO, 

ROOTS,  M.HNONRT.NONRT, NR 0*00, NCO,OELTA,ERSI 

29 

1LON.NCURI 

10 

IF  (NCURI  12,12,9 

11 

9 

IF  1NL-NCUR 1  11,11,10 

12 

10 

WL-NCUR 

11 

CO  TO  7 

19 

11 

M-N«l 

15 

IF  (HMAX-Ml  1,7,  7 

16 

12 

CALL  RECON  ( ROOT S , At  1 , 1 1 

•  IAI1,  ll.O.OEGREE) 

17 

CO  TO  1 

18 

CKO 

19- 

SUBROUTINE  KOOTSQ  (A,!A, 

NCUR.NNI  •  •••« 

1 

OIRENSION  A ( 5 1 , 1 1 •  IAIS1 

,11  RTSQ 

2 

N1-NCUR«1 

RTSO 

3 

DO  1  J-l.Nl 

RTSQ 

9 

AIJ.2 1 aA ( J ,  1 1 

RTSO 

5 

IAI J, ? 1 ■ 1  A  (  J,  1 1 

RTSQ 

6 

AIJ* J)«0.Q 

RTSQ 

7 

IAU.11-0 

RTSQ 

a 

1 

CONTINUE 

RTSQ 

9 

DO  10  M«1,NN 

RTSQ 

10 

OO  7  J-l.Nl 

RTSQ 

11 

ci-hi- j 

RTSO 

12 

K2-J-1 

RTSQ 

1) 

RM«XN 1  NOE  1 K 1 , K2I 

RTSO 

|9 

IF  (RMI  2,5,2 

RTSQ 

15 

2 

DO  5  L-l.KM 

RTSO 

lb 

CA-XHOOE (L,2> 

RTSQ 

17 

JC-J-L 

RTSQ 

ID 

JCF»J»L 

RTSQ 

19 

IF  ILRi  3,1,9 

RTSO 

20 

1 

S-AI  •’  ,2  1  *  A 1  JLR,  2  1 

RTSQ 

21 

■Z-IAI JL.2I *1 A(JLP,2I 

RTSO 

22 

CALL  SCALE  IX.IXI 

RTSO 

2  i 

CALL  AOO  IAI J,ll ,IAIJ,1) 

,X,IX,AIJ,3),IA(J«11)  RTSO 

29 

CO  TO  5 

RTSO 

2> 

4 

X-AIJL,2)«A(  JLR.2! 

RTSQ 

26 

71 


IX- IAIJL.2I-I A(JLF,2t 
CALL  SCALE  IX, IX) 

CALL  S8TRT  I A 1 J, 3 ) ,  I  Al  J, 3 1 , X,  1  X, Al J, 3 ) , I A<  J,  3)1 
3  CONTINUE 

Al J,3)-2.0*AI J . 3 1 

CALL  SCALE  U(  J,3I,IAI  J,3t) 

X- AI J«2) -*2 

IX- IAt  J,2)*IAU,2> 

CALL  SCALE  IX, IX) 

CALL  ADO  IAI J.3) ,  I  Af  J,3I  , X,  IX  ,  Al J , 3 1 , I AIJ,  3 1) 

JR-XHOOFI J,2I 
IF  IJRI  6,6,7 
A  AIJ, 31— AIJ, 31 
7  CONTINUE 

IF  IMM-HI  10,10,8 

•  00  9  J-I.Nl 
AIJ,2)-A( J,3> 

IAIJ.2I-IAI  J.3I 
AIJ, 31-0.0 

IAI J, 31-0 
9  CONTINUE 
10  CONTINUE 
RETURN 
END 

SUBROUTINE  REALM  ( A , I  A, K ,NCUR , OELT A, EFS I LON.ROHOO, MROMOO,  NONRT, 
1 ONR  I ,  NC  0  ,  R  OU  T  S I 

DIMENSION  AI51,}),  lAISl.l),  ROOT S I  2 , 30 1  .  ROM0DI30I,  NROMOu(SO), 
10NRTI SO) ,  NNONRT 1501,  RATI0IS1),  I P IV I  5  1 1 •  AREDI30),  IAREDI30) 
RATIOI ll-l.o 
DO  6  I-2.NCUR 
I i-XKGOF 11,21 
IF  I A 1 1 , 3  )  I  2,1,2 

1  RATI0in«0.0 
CO  TO  6 

2  T-Atl, 21-All, 21 
IT-IAII .2I-IAII.2I 
CALL  SCALE  IT, IT! 

T-T/AI I ,3) 

IT- IT-IAII, 31 
IF  tIT-21  3,3,1 

3  IF  1IT-21  1,6,6 

6  CALL  UNSCALc  IT, IT) 

RATIOI I  I - T 
IF  111)  3,3,6 
3  RATIOI I  I --RATIOI  I  I 

6  CONTINUc 
RATIOINCUR-D-UO 
IFIVID-I 

IPI VI NCUR- 1 ) ■ l 
00  9  I • 2.NCUR 

X- AASFIRATI Jl  I  1-1.0) 

IF  IX-OELTAI  7,8,8 

7  IFIVI  M-l 
CO  TO  9 

•  IFIVin-0 
9  CONTINUE 

NCURl-NCUR-1 

11-0 

MULT-0 

1-1 


RTSO  27 
RTSO  28 
RTSO  29 
RTSQ  30 
RTSQ  31 
RTSQ  32 
RTSQ  33 
RTSQ  36 
RTSQ  33 
RTSQ  36 
RTSQ  37 
RTSQ  38 
RTSQ  39 
RTSQ  60 
RTSQ  61 
RTSQ  62 
RTSQ  63 
RTSO  66 
RTSO  63 
RTSQ  66 
RTSQ  67 
RTSO  68 
RTSQ  69 
RTSQ  30- 
MN—  •••  1 

NRLRT 
ALM 
RLRT 
RLRT 
RLRT 
RLRT 
RLRT 
RLRT 
RLRT 
RLRT 
RLRT  13 
RLRT  16 
RLRT  13 
RLRT  16 
RLRT  17 
RLRT  18 
RLRT  19 
RLRT  20 
RLRT  21 
RLRT  22 
RLRT  23 
RLRT  26 
RLRT  23 
RLRT  2b 
•  RLRT  27 
RLRT  28 
RLRT  29 
RLRT  30 
RLRT  31 
RLRT  32 
RLRT  33 
RLRT  36 
RLRT  33 
RLRT  36 


2 

3 

6 

3 

6 

7 

8 
9 

10 

11 

12 
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I 


14-  1 

10 

I2-H*I 
MULT *PULT* 1 

IF  1 1 P I VI t  2  I )  10,10,11 

11  R0M00I 14  I  ■  A I 12,11/At  1,31 
I ROM 90* I •( 12, 31—1*11 ,31 

CALL  SCALE  (nOMOOU«l,IAOMOO) 

If  IROMOlH  141  )  12,13,13 

12  R0M0DU41--R0M90U41 

13  CALL  00UPL06  IROMOOC I4> , IROMODtXN,  IXN1 

T«2**MI1  » 

XN-XN/T 

CALL  SCALE  (XN.IXN) 

CALL  OLEXP  I XN,( XN,R0M00( 14), IROMOD) 

IF  IIR0M0D-74I  1  A, |A , 15 
1A  IF  ( I R  0*100*  74  I  15,15,16 

15  R0M09I 141-0.0 
1R0M00*9 

60  TO  1 7 

16  CALL  UNSCALE  (ROMOOI IAI.1ROMOO) 

IT  MROMOOH Al-MULT 

IF  I NCUR* 1-121  19,16,18 
18  1-12 

I  A* I  A ♦  1 
MULT-0 
11-0 

60  TO  19 

19  0-0.0 
NC0-9 

00  28  1*1, 1A 
ML- 1  A* 1-1 

M  — ROPOO(KL) 

15- MR0M00IKLI 
00  26  J* l , I  5 
J-J 

20  CALL  TEST  IA, I  A, M ,0, NCUR i ROMOO ( KL )«fcPSILON,K) 

IF  (Kl  23,23.21 

21  ROOTS! l.NCURi *>M 

roots«2,n:ur»*o.o 

ARE0(l)*4tt,t) 

IAREOI 1 1 « ( A| 1 , 1 1 
00  22  L * 2 « NC U R 
V-AREOIL-l >•* 

JY-IAREOU-ll 
CALL  SCALE  IY,IYI 

CALL  S8TKT  I A  1 L, 1 ) ,1  A ( L , 1 1 , Y, I Y, AAEO(L » , I AREOC L )» 
AIL , 1 ) *AREO< L  ) 

IAIL, 11*1 AREO ( L I 

22  CONTINUE 
60  TO  25 

21  IF  IWI  24,27,27 
2A  X»-X 

CO  TO  20 

25  NCUR-NCUR-l 

26  CONTINUE 
60  TO  28 

27  NCO*NCO* 1 
NONR T INCO I *RL 
MN0NRTINC01-I5M-J 


RLRT  37 
ML AT  38 
RLRT  39 
RLRT  AO 
RLRT  A 1 
RLRT  A2 
RLRT  A3 
RLRT  AA 
RLRT  AS 
RLRT  A6 
RLRT  A 7 
RLRT  48 
RLRT  49 
RLRT  50 
RLRT  51 
RLRT  52 
RLRT  53 
RLRT  54 
RLRT  55 
RLRT  56 
RLRT  57 
RLRT  58 
RLRT  59 
RLRT  60 
RLRT  61 
RLRT  62 
RLRT  63 
RLRT  64 
RLRT  65 
RLRT  66 
RLRT  67 
RLRT  68 
RLRT  69 
RLRT  70 
RLRT  71 
RLRT  72 
RLRT  73 
RLRT  74 
RLRT  75 
RLRT  76 
RLRT  7/ 
RLRT  78 
RLRT  79 
RLRT  89 
RLRT  81 
RLRT  82 
RLRT  83 
RLRT  84 
RLRT  85 
RLRT  86 
RLRT  87 
RLRT  88 
RLRT  89 
RLRT  90 
RLRT  91 
RLRT  92 
RLRT  93 
RLRT  94 
RLRT  95 
RLRT  96 


i 
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28  CONTINUE 

RLR.T  97 

HE  HUN 

RLRT  98 

INO 

RLRT  99 

SUBROUTINE  CGHFAT  (A , I  A* RONOD, ROOTS .M.NNONRT , NONRT , MRONOC.NCO,  DEL !•*•••  1 

1A, EPSILON, NCURI 

»***♦  2 

OiMENSION  A  t  b  1  •  S  I  •  1  A I  5 1  ,  3 )  ,  RONODI  SOI  i  KOOTSI2.SO).  SRI9I.3I,  I  SR  3 

1191, 11,  SACNOOISO),  SR00TSI2.S0I ,  KNONATI50), 

NONRT ISO),  MSROHOOI 5  4 

201  •  NS0NAT<<<4U  MSN0RTI49I,  HROMOOOO),  0(2). 

RI2I,  8(49)  5 

00  2R  1*1. NCO 

6 

JA-NONRT 1 1  I 

7 

ll-MNONRTI 1 1 

8 

11*11/2 

9 

IF  111)  10.2 

10 

1  11-1 

11 

2  IF  IROMOO(JA))  3.28.3 

12 

S  Q-ROMOOIJA) 

13 

00  27  J-l.tl 

14 

CALL  SUBRES  1  A. 1  A .NCUR, SR. 1  SR.  01 

15 

IF  INCUR-41  5,4,6 

16 

A  NSCUR-2 

17 

CO  TO  7 

18 

9  NSCUR-l 

19 

J2-1 

20 

CO  TO  • 

21 

A  NSCUR-NCUR-3 

22 

T  J2-NSCUH 

23 

•  LL-NSCURU 

24 

IF  INSCUR-ll  9,9,11 

25 

9  IF  (SAll.tl)  10,12,10 

26 

10  X-SRIl.D 

27 

IX-ISKfl.il 

28 

V*SRI2*  1 1 

29 

IY-ISRI 2. 1 ) 

30 

CALL  UNSCALE  IX, 1X1 

31 

CALL  UNSCALc  (Y, TV) 

32 

SROOTSI 1.1) *“  Y/X 

33 

NSC UR-0 

34 

CO  TO  13 

35 

11  CALL  ROOTSO  (SR, ISR.NSCUR.M) 

36 

CALL  REALRT  ( SR.  1  SR. M.NSCUR. DELTA. EPS  ILON. SRONOO.RSRONOO.NSONRT.HS  37 

1N0R7.NSC0, SCOOTS) 

38 

IF  U2-NSCUA)  12,12.13 

39 

12  SROOTSI 1 ,  J2I  0,0 

40 

IS  SROOTSI 1.J2I" SROOTSI 1, J? I*ROHOOI JA) 

41 

T-ROMOOI JAItRUNQOUA) 

42 

IF  1 SROOTSI 1 , J2) *  T )  1A, 21.21 

43 

14  W-SROOTSI 1 , J2 ) 

44 

WE-AONOOI  JAIMOHOOI  JA) 

45 

CAM.  TEST  (A.IA.W.WE , NCUR .RONODI JA 1 , ERS ILON.R I 

46 

IF  IK*  20.20,15 

47 

19  ROOTS!  l.NCUAI—uiO.O 

48 

1-4.0-WE 

49 

U-W*W 

50 

T-T-U 

91 

IF  IT*  14,16,17 

1*  T— T 

53 

U-SORTIT) 

94 

ROOTS  11, NCURI -ROOTSO, NCURI -U/2.0 

55 

ROOT  SI  1.NCUH-D--IW-UI/2.0 

56 

ROOT SI  2, NCUR 1-0*  0 

57 

74 


■.4 


ROOTS! 2.NCUR-1 1*0*0 

56 

CO  TO  16 

59 

17 

U-SQRTITI 

60 

ROOTS  ( 2  »NCUR1 -U/2. 0 

61 

ROOTS!  l,NCUR-ll-ROOTS(l,NCURI 

62 

ROOTS!  2.NCUR- 11  --ROOTS!  2  .NCURI 

61 

It 

Dlll-W 

64 

0121-wE 

65 

CALL  QUAD  1 V  INCUR, At  IA,R,0,6) 

66 

JX-NCUR-1 

67 

00  19  JY- 1 i JX 

66 

A(JY,11-B( JYI 

69 

IA( JY.11-0 

70 

CALL  SCALE  (A! JY  ,  1 1 ,  1  A(JY , 1 1 1 

71 

19 

CONTINUE 

72 

NCUR-NCUR-2 

71 

CO  TO  27 

74 

20 

W— W 

75 

CALL  TEST  (A.1A.W, WE. NCUA.ROMOOIJAl, EPSILON, R> 

76 

IF  <KI  21,21.15 

77 

21 

IF  ( J2-(NSCUR,lll  26,22,24 

76 

22 

IF  ( J2- 1 1  26,28,21 

79 

21 

J2-J2-1 

60 

SR00TSI1, J21-0.0 

61 

CO  TO  14 

62 

24 

IF  (SR00TSII.J2I -SROOTSI 1 « J2-1 1 1  25,26,25 

61 

25 

J2-J2-1 

64 

CO  TO  11 

65 

26 

J2-J2-1 

66 

CO  TO  21 

67 

27 

CONTINUE 

66 

2* 

CONTINUE 

69 

RETURN 

90 

ENO 

91- 

SUBROUTINE  TEST  1  A,I  A,  W,  0  .N.RON'jO,  EPSIL  ON.K  1 

•  •••« 

►  1 

DIMENSION  AI51.1I,  IAI5l.ll,  6(11,  10111,  TI2),  EI2I,  CI5I) 

TEST 

2 

6(11-0.0 

TEST 

1 

IX-0 

TEST 

6 

IW-O 

TEST 

5 

16(11-0 

TEST 

6 

6(21  *A( 1 , 1 1 

TEST 

7 

I6I2I-IAI 1,11 

TEST 

6 

00  2  l-l.N 

TEST 

9 

X-W»6(2t 

TEST 

10 

IX- IB  121 

TEST 

11 

CALL  SCALE  (X.IXI 

TEST 

12 

V-0*B ( 1 1 

TEST 

11 

IY-IMII 

TEST 

1 4 

CALL  SCALE  IY.IYI 

TEST 

15 

CALL  ADO  IX,IX,V,IY,2,IZI 

TEST 

16 

CALL  S6TRT  1 A( I* 1 , 11 , 1  At  1 «1 , 1 1 ,2 ,  IZ  ,61  1 1 , 1 6(  11 1 

TEST 

1  T 

IF  (N-ll  2,2,1 

TEST 

Id 

1 

6(11-6(21 

TEST 

19 

16(11-16(2) 

.  TEST 

20 

6121-6(11 

TEST 

21 

16(21-16(11 

TEST 

22 

2 

CONTINUE 

TEST 

21 

ROUNT-1 

TEST 

24 

CEPSIL-EPSILON 

TEST 

25 

Tlll-O. OVTI2I-0.0 

TEST 

26 

75 


N1*N* | 

TEST 

27 

X*2«C*Ef SILON 

TEST 

28 

V*X*AOMOO 

TEST 

29 

Cm«K0MC0»T 

TEST 

SO 

1 121 -ft0N00»C£fSIL 4X0*00 

TEST 

il 

00  4  1*1,41 

TEST 

12 

if  uii, m  s,4,* 

|,‘  ST 

3) 

3 

cm->4<  i.iuicmaii  ,n 

TE.'T 

)4 

SO  TO  1 

TEST 

15 

4 

C(l>*A||,|>llC*tAII,l> 

TEST 

It 

S 

CALL  ON SCALE  ICIII.ICI 

TEST 

17 

mi*Tiii*ein»cm 

TEST 

If 

Tt 2I*TI 2  )*E 12 l«C 1 1 > 

TEST 

34 

6 

CONTINUE 

TEST 

40 

Olf*Tll)-TI2l 

TEST 

41 

if  ioi  la.r.ii 

TEST 

42 

y 

If  18(11  1  1.4,9 

TEST 

4) 

• 

Bl)l--B<)> 

TEST 

44 

« 

If  11801-74)  10,10,12 

TEST 

45 

10 

If  1 1 8 1 1  •  ♦  74 )  12.11.11 

TEST 

4* 

n 

CALL  UNSCALE  181)1, 141111 

TEST 

47 

if  ioif-emi  i2. 12, ir 

TEST 

4  H 

12 

(•0 

TEST 

49 

If  IKOUNT-2)  1), 14,14 

TEST 

50 

IS 

If  III  14,14,15 

TEST 

SI 

14 

If  INI  15,16,1* 

TEST 

S2 

IS 

SENSE  LIGHT  2 

TEST 

S> 

KOUN7-KOUNT*! 

TEST 

54 

14 

AETUAN 

TEST 

SS 

l» 

8*1 

TEST 

Sfc 

00  TO  14 

TEST 

S  7 

14 

If  IIBI2I-T4I  14,14,12 

TEST 

5* 

14 

If  (18121*74)  12,20.20 

TEST 

59 

20 

CALL  UN SCALE  (8121,181211 

TEST 

40 

If  IIBI II* 74|  21.21.12 

TEST 

41 

21 

If  llf(SI*74l  12.22.22 

TEST 

42 

22 

CALL  UNSCALE  181)1,181)11 

TEST 

4) 

*■0*8 121*812) 

TEST 

44 

v-NMi2i*em 

TEST 

45 

2*01 )l*OI )) 

TEST 

44 

M-r»i 

TEST 

47 

If  ITI  21,17,24 

TEST 

44 

21 

Y*-T 

TEST 

49 

24 

oif-oif •oif 

TEST 

70 

If  IDIf>VI  12,17,17 

TEST 

71 

(NO 

TEST 

72- 

1U8A0UT |KE  SU8AES  1 A , 1  A, N, SR, 1 SA , AONOD 1 

01  ME  NS  ION  All! ,) 1 ,  1  A I  SI , )) ,  SR(41,)I,  1SA(5I,)I,  CI51I.  AI50,)) 

2 

Rl-N*l 

1 

T*1,0 

4 

00  1  1*1, N 

5 

J*Nl-l 

4 

T*T*XONOO 

7 

CUI*AIJ,II*T 

• 

IC-IAIJ.ll 

4 

CALL  UNSCAU  ICUI.ICI 

10 

1 

CONTINUE 

11 

CINlt-AINl.l) 

12 

1C*  1  Al N 1 ,  1 ) 

11 

CALL  UNSCALE  ICINII.ICI 

14 

76 


t 


If  14-21  17,17,2  IS 

2  N2-N-2  16 

DO  3  l-l,N2  17 

•(1,11-0.0  16 

6(1. 21-0.0  19 

•  CONTINUE  20 

1-2  21 

6(1, 2) -C  111  22 

•  61 1,31-C  I  11-61 1,1 1  23 

00  9  J-2.N2  24 

8(  J,3I--6U-1,2I-A(J,1I  29 

9  CONTINUE  26 

If  (N-l 3*1 1 1  6,6,6  27 

6  1-1*1  26 

DO  7  J-1.N2  29 

6( J, 1 1 -6 ( J, 2 )  30 

.  61 J.21-HI  J,  3 1  31 

T  CONTINUE  32 

CO  TO  4  33 

•  IF  IN-41  19,9,13  34 

9  IF  (N-I2*l)»  12,10,10  39 

10  1-1*1  36 

DO  11  J-1,2  37 

6IJ,II-B( J.2I  36 

6(J, 21-61 J.JI  39 

11  CONTINUE  40 

CO  TO  4  61 

12  6(3, 3I--6 (2,21  62 

SR  (3. II- '~C  (91*811, 31  63 

I  SRI  3, 1 1 *0  66 

SRI2, 11-6(2,31  69 

ISRI2.il -0  66 

SR(1. 11-0(3, 31  67 

1  SR (1,11-0  66 

CALL  SCALE  (SR(l.l>,ISR(l,lll  69 

CALL  SCALE  ( SA( 2 , 1 1 , I S*( 2 , 1 1 1  90 

CALL  SCALE  (SA(3,1I,ISA(3,UI  91 

CO  TO  16  92 

13  SRIN2,ll-C(Nl-8(l,3l  93 

1  SR ( N2 , 1  I  -0  96 

SRIN2-1,  11—  CINll-612,31  99 

I  SRI N2- 1*11*0  96 

CALL  SCALE  (SA(N?,11 ,ISA(N2,1) 1  97 

.CALL  SCALE  (SA(N2-l,l),ISR(N2-l,lll  96 

IF  IN2-21  16,16,14  99 

16  OO  19  J-3.N2  60 

R-N2* l- J  61 

SR(«,1I--6U,1I  62 

I  SR (K, 11-0  63 

CALL  SCALE  (SR(K,II,ISR(K,1II  64 

19  CONTINUE  69 

16  RETURN  66 

17  SRIl.ll-CIll  67 

I  SR (1,11*0  66 

SR12, 1I--CI2I  69 

I  SR ( 2 , 1 1 -0  70 

16  CALL  SCALE  ( SR  1 1 , 1 1 . I  SR  1 1 , 1  1 1  71 

CALL  SCALE  ISRI2,1I,ISRI2,1I1  72 

CO  TO  16  73 

19  SR f 1 , 1 1 --C (4 1  74 
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SMa.n-ctsi-cit  i 
I  tMi.no 

ISM2.1IO 

60  to  to 

CNO 

SUBROUTINE  RECON  (ROOTS. 4.I4.0.NI 
DIMENSION  ROOTS! 2  .501 »  DISH 
!•« 

IKU 

CRU  UNSC4LE  (1*1X1 
00  I  l-I.N 
0(11-0.0 
1  CONTINUE 
OIN-ll-l.O 
l-l 

Nl-N-l 

a  IE  I  ROOT  S I  2 . 1 ) I  3,7.) 

3  T-ROOTSII .1 l-ROOTSll.l  I 
U-R00TSI2. 1  I -ROOTS (2. I  I 
T-T-U 

U-2.0-RO0TSIl.il 
00  S  J-l.NL 
IE  ( I ♦  J-NI  0.4.4 

4  OIJI-OI J-2I-T-0I  Jl 
OIJI-DI  JI-U-OIJ-D 

5  CONTINUE 
OINI -T-UINI 
OINI-DINI-U-OIN-tl 
OIN-I l-T-DIN- 1 1 
1-1-2 

4  IE  IN-13  10.2.2 
700  I  J-I.N 

IE  IJ-I-NI  4,9,4 
4  D(JI-0(J-ll-0(JI-ROOTS(l«ll 
9  CONTINUE 

OIN-ll—  OIN-U-ROOTSI I  .11 
1-1*1 
CO  TO  4 
to  NS-N-l 

00  II  ll-I.NS 

oiin-onii-x 

rt  CONTINUE 
RETURN 
(NO 

SUBROUTINE  QUAOIV  IN, 4,  I4,R,D,  B> 
DIMENSION  MSI. 31.  14(31.11,  *121, 
4(11-4(1.11 
14-14(1,11 

C41L  UNSC4LE  14(11,141 
IE  IN-21  4,4,1 

1  44-4(2.11 
144-14(2,11 

C4U  UNSC4LE  1*4,1441 
4121-4**4(11-0111 
IE  IN-31  4,4,2 

2  NT-N-I 

00  3  1-1.  T 
XN-RI l-ll-OI I  I 
TN-BI I -2  I -U 1 2  I 
44-411,11 


TS 

TO 

TT 

T4 

79- 

•••••  t 

a 

3 

4 

5 
0 
7 
• 
4 

10 

tt 

12 

13 

14 

15 
10 
IT 

14 

19 

ao 

at 

22 

23 

24 

as 

20 
27 
26 

29 

30 

31 

32 

33 

34 

35 
30 
37 
34 
39- 

•  MM  1 

*  4(491  2 

3 

4 

5 
0 
7 
4 
9 

10 

11 

12 
13 
1* 

15 
10 
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i 


IAA-IAII  ,11 

CALL  UNSCAlE  IAA.IAAI 

All l*AA-l XN»VN> 

J  CONTINUE 
4  XN*B IN-  11*0(11 
VN*6 1 N-2 1 *01 2  I 
AA*AIN,tl 
| AA ■ I A I N 1 1 1 
CALL  UNSCALE  IAA.IAAI 
Rill • A A- I »N*YN) 

AA* A  I N* 1.11 

|AA*IAI'I»1.1I 

CALL  UNSCALE  IAA.IAAI 

RI2I*AA-BIN-1I*0I2I 

RETURN 

END 

SUBROUTINE  OOUBLOG  IX.IX.Y.IVI 

1*44.0 

IF  1X1  1.2,1 

1  MINI  5 

2  Y-0.0 
I  Y*0 

CO  TO  4 
J  TO* I  X 

V* ALOG I  X I  *T0* ALOC I Tl 
I  Y*0 

call  scale  iv.iri 

4  RETURN 

5  FORMAT  I46H0THE  LOG  OF  A  NON-ROSITIVE  NUMBER  IS  REOUESTEOI 
END 

SUBROUTINE  OLEXR  I*. IX, 2,121 
2-EXP  I  XI 
12*0 

IF  11X1  1.11,4 

1  I  — IX 
11*4*1 

00  4  J-l.lt 
RaXMODF | 12,21 
IF  IKI  2,1,2 

2  12*12-1 

2*44.0*2 

S  12*12/2 
2-SQRII2I 
CALL  SCALE  12,12) 

4  CONTINUE 

5  RETURN 
*  I *4* I  X 

oo  r  j.i.i 
2*2*2 
12*12*12 

CALL  SCALE  12,121 

?  continue 
CO  TO  ) 

I  CALL  SCALE  12,121 
CO  TO  S 
INO 

SUBROUTINE  ADO  I  X, IX, Y, IV, 2*12  I 
IF  IXI  3,1,1 
I  2*Y 


17 

10 

IN 

20 

21 

22 

23 

24 

25 
24 
27 
20 
21 
30 

11 
32- 

•#•••  | 

2 

3 

4P 

3 

4 
7 
0 
» 

10 

It 

12 

13 

14 
14- 

•••••  l 
2 

3 

4 
3 
* 

7 

• 

N 

10 

11 

12 

13 

14 

13 

14 
17 
10 
1* 
20 
21 
22 

23 

24 
23 
24- 

••••#  I 
2 
3 
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IX- 1 V 
2  RETURN 
$  IF  m  5,4,5 

4  t>X 
I2-1X 
60  TO  2 

5  10IFF-1X-IV 

IF  (lOIFF)  6,1,7 

4  UMY 
A-V 
B-X 

ID1FF  — 10IFF 
60  TO  • 

T  I A— I X 
A-X 

b*v 

•  IF  116-lUlfFI  4,0,10 

•  2- A 
12*1* 

60  TO  2 

10  IF  IIOIFF)  11,11,11 

11  00  12  1*1,101 FF 
0-A/64.0 

12  CONTINUE 
11  CONTI NUc 

M»1 
12-1 A 

CALL  SCALE  (2,121 

CO  10  2 

ENO 

SUBROUTINE  SBTRT  ( X, I  X  ,  V  ,  IV ,2 , 12 ) 
N—V 

CALL  ADO  (X. IX, W, IV, 2. 121 

BE TUAN 

ENO 

SUBROUTINE  SCALE  IX, I XI 

KEC64-1.G/64.0 

IF  (XI  1,11,2 

1  V—  X 
60  TO  1 

2  V*X 

1  IF  (64.0-Y)  4,5,1 
A  V*Y/64, 0 
U*IX«1 
60  TO  1 

5  IF  I V-REC44 }  6,7,7 

6  V>V»64.0 
IX-IX-l 
60  TO  5 

7  IF  (XI  8,4,9 

•  x—v 

60  TO  10 
4  X*V 

10  RETURN 

11  lx*0 

60  TO  10 
ENO 

SUBROUTINE  UNSCALE  IX, Ul 
IF  (IX4R4I  1,2,2 
I  X-O.O 


4 

5 

6 

7 

B 

4 

10 

11 

12 

11 

14 

15 

16 

1 7 

18 

19 

20 
21 
22 
21 

24 

25 

26 

27 

28 

29 

30 

31 
12 
11- 
'  1 

2 

1 

4 

5- 
1  1 
2 

3 

4 

5 

6 
7 
B 
9 

10 

11 

12 

11 

14 

15 

16 

17 

18 

19 

20 
21 
22- 

1 

2 

1 
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IK-0 

60  ro  * 

IF  IIX-B4I  4.4,3 
X-1.0E*t33 
IX»0 
PRINT  T 
60  10  6 
IF  1111  9,6,3 
X>X»64.0**IX 
!X«0 
RETURN 

FORMAT  I25HOEXP.  OVERFLOW  IN  UNSCALE  I 
END 

SUBROUTINE  CHECK  1  <V,AC,PP1 ,PP2,PP3> 
DIMENSION  Vlftl  •  ACm 

viimviimexpippii 

VI2MVI2MExi><-oPl  1 

AV3*£XP|PP2MIVI  3MC0SIPP31  ♦  V|4MSIN|PPm 

AV4«EXP|PP2MI  VI4MC0SIPP3I-VI  3  l*S  INI  PPTM 

AVS*f  XPI-PP2I  •!  VI5MC0SIPPlMVl6MSINIP»m 

AV6-EXPI-PP2MI  VlftMCOSIPPlI-VMMSIMPPm 

VC  JMAVJ 

VI4MAV4 

VI5MAVS 

VI4MAV6 


RETURN 

END 

SUBROUTINE 


CMECK2  CV,AC,PPl.PP2,PPJ) 


4 

5 

4 

7 

• 

9P 

10 

11 

12 

13 

14 
19 
16 


CHECK 
CHECK 
CHECK 
CHECK 
CHECK 
CHECK 
CHECK 
CHECK 
CHECK  1 0 
CHECK1I 
CHECK  12 
CHECK1 1 

Check i 4 

CHECK  1 3- 

•••••  1 


1 

2 

3 

4 
9 
6 
T 
6 
9 


DIMENSION  V  1  6  1  •  ACHI 

2 

AVMVI IMC0SIPP1MVI2MSINI  PPl  1 

3 

AV2«-VIIMSIMPPIMVI2MC0SIPPI> 

4 

VIIMAVl 

5 

VI2MAV2 

6 

AV3«fc‘XP|PP2MI  VI  9MC0SIPP3MVC4MSINIPP3II 

7 

AVMEXPIPP2MI  VI4MC0SIPPSI-VI  3MSINIPPM  1 

8 

AV5-EXPI-PP2MIV1SMC0SI  PPl  MV  1  6 MS  INI  P»  ill 

V 

RV6-EXPt-PP2MIV|6MC0S(PP3MVI9MSIN(PPm 

10 

VI 3) »AV3 

11 

VI4MAV4 

12 

VlSMAVS 

13 

VI6MAV6 

14 

RETURN 

19 

ENO 

16- 

SUBROUTINE  CMECK4  IV, AC, PPl .PP2.PP3I 

••••i 

►  1 

OINcNSI  ON  VI  6  1  ,  ACID 

CHK4 

2 

VIIMVIIMEAPIPPII 

CHK4 

3 

VI2MVI2MEXPI-PPH 

CHK4 

4 

AV1»VI3MC0SIPP2MVI4MS1N<PP2I 

CHK4 

5 

AVM-Vt  IMSIMPP2MVI4MC0SIPP2I 

CHK4 

6 

AV5»V  IS  MCOS 1  PPl  MVI6MSINI  PP3  1 

CMK4 

7 

AV6  — VI3MSINIPP3MVI6MC0SIPP3I 

CHK4 

V 

VI3MAV3 

CHK4 

9 

VI4MAV4 

CHK4 

10 

VlSMAVS 

CMK4 

11 

VI6MAV6 

CHK4 

12 

RETURN 

CHK4 

13 

ENO 

CMK4 

14- 

SUBROUTINE  CHECKS  IV,AC,PP1,PP2,PP3I 

•  »••• 

1 
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01 MCNSI ON  VI4I,  ACM) 

CHKS 

2 

CHKS 

1 

CHKS 

A 

V 0 1  •  V (  1>*£KMFF2I 

CHKS 

S 

V|4»*V|4I*CKFI-NPM 

CHKS 

6 

AVS*VISI*COSI  AM  »  ♦  VI 4  (•SINtERlI 

CHKS 

T 

AV4*-V(*I«SIN(FFM»VI4I«C0S(FF)I 

CHKS 

6 

MSI-AVS 

CHKS 

9 

vuiovi 

CHKS 

10 

AIT URN 

CHKS 

1  1 

I  NO 

CHKS 

12- 

SUBROUTINE  mica  iaoot.ac.rlifi 

1 

OIMENStON  IC9I6I.  IIMI6I,  IAEI6I 

RICE 

2 

DIMENSION  ROOT  12.41.  ACIlt 

RICE 

1 

COMMON  /KIF/  ICl  .  IC2, IC1, IC4, ICS, IC4,  |Ml,  IH2, 

IMS, IMA. IMS.IH6, 1 A 1 , 1  A  1 CE 

A 

|R2,IR1,M4,IRS,<R6 

RICE 

S 

EQUIVALENCE  MCI.ICOMII,  I|M1,||M«  lit.  IIM. 

IREIMI  RICE 

6E 

MM  IE  19.291 

RICE 

rw 

MM  IE  14.211  MA00TM,JI.I*l,2l,J*l.4t 

RICF 

SM 

00  1  J*l  ,4 

RICE 

9 

ICOUI-O 

RICF 

10 

IIHUt-0 

RICE 

1  1 

IRE  1 J 1  *0 

RICE 

12 

1 

CONTINUE 

RICE 

1 1 

IRO 

RICE 

14 

IC-0 

RICE 

IS 

11-0 

RICF 

1  6 

00  S  0*1.6 

RICE 

1  T 

IF  IA6SIH00TI 1 .JII.lT. 1.E-12I  ROOT  1 1 . J 1 *C. 0 

RICE 

Id 

IF  <ARS<k00T(2,J) I.LT.l.E-12)  ROOT  1 2 . J 1 *0. 3 

RICE 

19 

IF  IR00M2. JI.EQ.0.91  CO  TO  4 

RICT 

20 

IF  IROOTII.JM  2.1.2 

RICF 

21 

X 

IC-ICM 

RICE 

22 

ICOIICI-J 

RICE 

21 

CO  TO  S 

RICE 

24 

) 

Iv-IMI 

RICE 

2S 

MMM  II  ■  J 

RICE 

26 

CO  TO  S 

RICE 

2  T 

4 

IR-IRM 

RICE 

2a 

IRE  1 IR 1 • J 

RICE 

29 

s 

CONTINUE 

RICE 

10 

MRITE  19,221  IR. 1 1 . 1C 

RICE 

liw 

IF  IIR.EC.2.AN0. IC.EQ.4I  CO  TO  4 

RICE 

12 

IF  IIR.EQ.2.ANU. IM.EU.4I  CO  TO  1 

Ricr 

11 

IF  |  IN.E0.6I  00  TO  9 

RICE 

14 

IF  IIM.EQ.2.ANO. IC.EQ.4I  CO  TO  16 

RICE 

IS 

IF  MA.EC.4.  AND.  Tm.Eu.2I  CO  TO  IT 

RICF 

16 

IF  MR. ED. 61  00  TO  19 

RICE 

it 

CO  TO  19 

RICE 

Id 

RICE 

19 

RICE 

40 

* 

KLltM 

RICE 

41 

RICE 

42 

ACM  1  •  AC  SI  ROOT  1 1 ,  INI  M 

RICE 

4  1 

AC  1 21 *ARS 1  ROOT  1 1 , IC1 1 1 

RICE 

46 

AC  1  1I*ABSIR00TI2. ICl 1 1 

RICE 

AS 

RE  TUAN 

RICF 

46 

1 

KMF*2 

RICE 

4  1 

RICE 

4d 

ACMMAISIR09TII.IRMI 

RICE 

49 

1 


ACI2I-A6SIA00TI2, INI  >1 

RICF 

SO 

IF  I  AC  1 2  1  •  NC  •  ABS ( ROOT 1 2  *  I  H2  )  1 ) 

CO 

TO 

• 

RICF 

SI 

ACT  31 *AR SIROO T (2 , 1 H3 1 1 

RICF 

12 

RETURN 

RICF 

S3 

• 

AC(3)*ABS(ROOT(2,IM?l) 

RICF 

S6 

RETURN 

RICF 

S5 

RICF 

S6 

9 

RUP-3 

« 

RICF 

5  T 

RICF 

SB 

AC(tl-RBS(R00T(2,IHin 

RICF 

S9 

IF  IABSI*OOTI2.INlM.NE.ABS(ROOT(2,IN2m  GO  TO  10 

RICF 

60 

AC(2l«AnS(R00T(2, 1  M3  1 1 

RICF 

61 

IF  IABS(*00T(2,IN6n.NE.  ACill. 

ANO. ABS 1  ROOT  12, IN6) I.NE.ACI 211 

CO  TOR  1 CF 

62 

1R12 

RICF 

63 

CO  TO  l A 

RICF 

66 

10 

ACC2l«ABSIRDOTI2.IH2»l 

RICF 

66 

IF  (AnS(ROOTI2,lH3tl.EO.  ACmi 

CO 

TO 

11 

RICF 

66 

IF  (AOSIROOTI2, 1*3)1. fO. ACI2I ) 

CO 

TO 

13 

RICF 

67 

ACI 3)«AOS(ROOT(2, 1  H3  1  1 

RICF 

68 

RETURN 

RICF 

69 

11 

IF  (AeS(R00M2,IM4l|.E0.ACI2n 

CO 

TO 

U 

RICF 

70 

12 

ACI3l*ABS(RuOT(2, 1*6  1 1 

RICF 

71 

RETURN 

RICF 

72 

13 

IF  (ArS(K00T(2. 1*611. NE.ACI1II 

CO 

TO 

IS 

RICF 

73 

19 

AC(3)*A£S(R30T(2, IMS  II 

RICF 

76 

RETURN 

RICF 

76 

15 

ACC  3I-ABS(R COT (2,1*61 » 

RICF 

76 

RETURN 

RICF 

7  r 

RICF 

78 

10 

RLIP-6 

RICF 

79 

RICF 

BO 

AC(l>-ABS(ROOT(?, INI }) 

RICF 

8  1 

AC(2l*ABS(RujT ( l ,  IC1 1 ) 

RICF 

82 

AC(3)»ABS(Rk/0T  (2,  ICl  )  1 

RICF 

B  3 

RETURN 

RICF 

66 

RICF 

86 

IF 

BLIP'S 

RICF 

86 

RICF 

87 

AC(1»«ABS(R0CT (l , tRIII 

RICF 

6d 

ACI 3)«ABS(R00T (2 , 1*2 I ) 

RICF 

89 

IF  (ABSIROOH  l,l«2)).EP, ACI  111 

CO 

TO 

18 

RICF 

90 

AC(2l»ABS<R00l<l,|R2H 

RICF 

91 

RETURN 

RICF 

92 

18 

ACI7I-ABSIR00TI1, IR3II 

RICF 

93 

19 

RE  1 P  * l 7  76 

RICF 

96 

RETURN 

RICF 

96 

RICF 

96 

ricf  9r 
RICF 
RICF  99 

ricfioo- 

1 


2 

Ol*ENSION  FENI 3, 31  ,  VPI1I  3 

OiPENSION  OEM  3,  3)  ,  £9(6,  M,  VI61,  OUOUI 6  )  •  C «  3  •  1  I  * 

00  I  I  •  l ,  3  S 

OO  I  J- I , 3  6 

OFM  I  *  J )  *0,0  7 

1  CONTINUE  • 

00  2  I >1,3  9 


tO  FORMAT  I l OX , I  OH  ROOTS  I 
tl  FORMAT  (IX, £12, 6, 12X.E12.6I 

22  FORMAT  U0X,3H|R«|N,3X,)HM«|4,3I,1HIC>IA/I 
CNO 

SUBROUTINE  FINCLE  IEN.VI 
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DEN(I,1)-ENII*3,1>  10 

D£N(1,2l-ENf  10,3)  11 

OENII,3)«EN<iO,A)  12 

CII'l)-EN(lO,7>  13 

2  CONTINUE  15 

WHITE  (9,91  I (DENI  I, Jl.J-l, 31, 1-1,31  15W 

CAU  LEO  I0EN,C,3,1,  3,3,OET)  16 

Vtl)-Cll,l)  12 

VI 31 aC (2, t I  10 

V(  A I *C( 3,11  19 

WRITE  (9,10)  Vm,V(3),V(A|  20W 

00  3  (•1,3  21 

00  3  J-1,3  22 

EENI I ,  J) *0,0  23 

3  CONTINUE  2A 

00  A  1-1,3  25 

PEN! 1,1) -ENI 1,2)  26 

PEN! 1 ,2) -ENI I «  5)  27 

PEN( 1,31 *EN( I ,61  20 

A  CONTINUE  29 

C< 1« 1 1-cNI 1,71-ENl l,l)*Vll)-ENll,3)*V( 31-ENI  1, A)*V( A)  JO 

CI2»l)-cN(2,T)-ENC2,ll*VCll-ENC2»3l*V(3l-£N( 2, A) -VC  A)  Jl 

Cl  3, 1  l«EN( 3, 7) -ENI  ), 1 >*VI U-ENI J,3)«VI  3 1 -ENI 3,41 *Vl A)  J2 

WHITE  (9,9)  I (P£NI I , Jl , J>  l«  J ) ,  1*1,3)  33H 

WRITE  (9,61  C(l,l),C(2,lliC(3,ll  JAW 

CALL  LEQ  (FEN,C,3,  l,3,3,OET)  J5 

VI 2 ) -C I  1  ,  1 )  36 

VI 5 1 -C 1 2 , 1 )  37 

VI6I-C I  3, 1 )  36 

WRITE  (9.111  V(2)»V(5I*V(6)  39W 

00  5  1-1,6  AO 

OUOUI 1 1 -f  Nl 1 , 1 1 • V 1 1 ) *ENI I »2 1 *V 1 2 1 aENI 1 , 3 l*V I  31 aENI 1 1 A |AVI AIaENI I ,$  A1 

1IYV<M<EN<I  ,6l«Vl6)-ENll  ,7|  A2 

J  CONTINUE  A3 

WRITE  (9,71  lUUOU(l), 1-1,61  AAW 

RETURN  A 5 

A6 

6  FORMAT  (IX, 25H  FOR  THE  FEN  MATRIX  / 10X , 7HC( 1,1)-E12.5,5X, THCI  A7 

I2ill-E12,6,6X,7HC(3,1)-El2*5/)  A8 

7  FORMAT  IIX.36H  SOLUTION  CHECK  BY  SUBSTITUTION  IS  /I6EI2.5I)  A9 

■  FORMAT  I  I X,2aH  FENll.JI  BY  ROWS  IS/IJE1S.7))  50 

9  FORMAT  (IX.2CH  OENtt.JI  BY  HOWS  IS/I3E15.7I)  51 

10  FORMAT  I IX.SMVI 1  I *el2.5,2X,SHV< 31-c 12 . S , 2 X , 5HV ( A 1 -E 1 2. 5 1  52 

11  FORMAT  (iX,)HV<2>«E12.S,2X,SHVm-EI2.5,2X,*HV(6)-E12.5l  53 

END  56- 

SUBROUTINE  LEQ  ( A.B.NEOS.NSOLNS.IA. IO.OET I  •••••  1 

LINEAR  EQUATIONS  SOLUTIONS  FORTRAN  II  VERSION  2 

SOL VF  A  SYSTEM  OF  LINFAA  EQUATIONS  OF  THE  FORM  AX-6  BY  A  MODIFIED  3 

GAUSS  ELIMINATION  SCHEME  A 

5 

NEQS  -  NUMBER  OF  EQUATIONS  ANO  UNKNOWNS  A 

NSOLNS  •  NUMBER  OF  VECTOR  SOLUTIONS  DESIRED  7 

U  •  NUMOER  OF  ROwS  OF  A  AS  OEFINcO  BY  DIMENSION  STATEMENT  ENTRY  • 

IB  •  NUMBER  J F  RO»S  OF  B  AS  OEFINfcU  BY  UIMFNSION  STATEMENT  ENTRY  9 

AOET  -  DETERMINANT  OF  A,  AFTER  EXIT  FROM  LEO  10 

11 

oinensi on  aiu.iai,  bub.ibi  12 

NSI2-NEQS  13 

NBSI2-NSCLNS  |A 

NORMALIZE  EACH  row  BY  ITS  LARGEST  ELEMENT.  FORM  MARTIAL  CETERNT  15 
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OET-l.O  16 

00  6  1*1, NSt l  17 

8IG*AI I • 1 )  10 

If  INSIZ-ll  17,17,1  19 

1  00  3  J*2 ,NS1  Z  2° 

If  tAPSf  (BIOI-AASf  ( A  «  1*  JIM  2,3,3  21 

2  BIG** 1 1 ? J)  22 

3  CONTINUE  23 

8C>1. O/BIG  29 

00  9  J* 1  *  NS  I Z  25 

4  At  I.JI-AI I ,JI»BG  26 

00  5  J*1,NRSIZ  27 

3  6II,J>*B(I,JI *6G  28 

0ET«DET*BIG  29 

6  CONTINUE  30 

START  SYSTEM  REDUCTION  31 

NUMSYS*NSIZ-l  32 

00  16  l-l.NUHSrS  33 

SCAN  HAST  COLUMN  Of  CUAAENT  SYSTEM  fOR  LARGEST  ELEMENT  34 

CALL  THE  ROW  CONTAINING  THIS  ELEMENT,  ROW  NBGRW  35 

NN*IM  36 

•  IG«AM,M  37 

NBGRW* I  38 

00  8  J*NN, NS  I L  39 

If  lABSf IBIGI-A6SE IA(J,I 111  7,8,8  40 

7  BIG** I J,  I  I  41 

N8GRW-J  42 

•  CONTINUE  43 

•GH.U/BIG  49 

SWAP  ROW  I  mITH  ROW  NBGRW  UNLESS  1-NBGRW  45 

If  INBGRW-ll  9,12,7  46 

SWAP  A-MATRIX  ROWS  47 

9  00  10  J* I , NS  I Z  48 

TEMP*AINBGRw, Jl  49 

AINRGRW,JI*AII.J)  50 

10  Al I, JI-TE-P  51 

OET.-OET  52 

SWAP  e-NATRIX  ROWS  33 

DO  II  JH.NOSIZ  54 

TEMP*enOGRw.  J»  •  55 

8IN8GRW.JI *81  I ,JI  56 

11  81 1 , J I ■ TEMP  57 

ELIMINATE  UNKNOWNS  fROH  FIRST  COLUMN  Of  CURRENT  SYSTEM  58 

12  00  15  K*NN,r.SIZ  59 

•COMPUTE  PIVvTAL  MULTIPLIER  60 

PMULT*-AI  K,  I  I  *8G  61 

APPLY  PMuLI  Tj  ALL  COLUMNS  Of  THE  CURRENT  A-MATRIX  ROW  62 

00  13  J**.N,NS1  2  63 

13  Al  R  ,  Jl  *PMUL  T  •  1 1  *  J I  ♦  A I  K,  J  I  64 

APPLY  PMULT  TO  ALL  COLUMNS  Of  MATRIX  B  45 

00  14  L-l.NOSI*  66 

14  II K.L I *PMUL  T*0tI,LI»6(K«L)  47 

13  CONTINUE  68 

16  CONTINUE  69 

00  BACK  SUBSTITUTION  70 

WITH  R-NATA I  X  COLUMN  *  NCOLB  71 

17  00  22  n:glb*i,nbsiz  72 

00  FOR  ROW  ■  NROM  73 

00  21  l»l,NSU  24 

RAOWNSl/M-t  75 
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M  .M-IOUSl.  C„»VU»  •  ««« 

Sg'S'MlSfl*  .0..."  >M 

IF  fNXSI  IB.20.1B 
NO 

II  00  14  fc-t.NXS 

II  T|NF«TE«B*MXX,HCOLM**«NWN.KIl»  |/mn)10m,nR0W1 

20  KNROW.NCOtBI-IBINROU.NCOlBI-TE  tOLUHN  •  NCOLB 

HAVE  VE  FINISHED  AU  BOVS  FOR  B-HATHI* 

21  CONTINUE 

”«  .1  JUST  TIN!  SUED  -It"  •-«'*•«  COlum  -MS.-SI1 
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APPENDIX  B 
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TABLE  B.l 


AVERAGE  MATERIAL  PROPERTIES  FOR  PG  AND  ATJ  GRAPHITE 


PG 

ATJ 
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3.1  x  1 0G  psi 

2.26  x 

106  nsi 

V 

-0.21 

+  0.30 

vc 

+  0.90 

+  0.25 

R 

varies  with  case 

30  in. 

L 

40  in. 

• 

40  in. 

a 

1.43  x  10‘6  in_ 

*v. 25  x 

10"6  jn__ 

Or 

in  -  F 

in  ■ 

ac 

13.1  x  10“G  in 

4.25  x 

10’6  in_ 

u_ 

o 

1 

c 

•r- 

in 
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Figure  B  I  Ef’ect  of  E/Gc  on  Joint  Shear  Stress 
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APPENDIX  C 


DEFINITION  OF  TERMS 


i 


^•j  -  (Cg  +  Cb)  D  {(ai^  +  *  tt-j  +  tt g )  ”  2  ®  ( °47  ^  fl4) } 


+  {£  CbD2  +  (■  R--dR”“b) }  (a37  +  a47  +  "3  +  V 


'a  b 


Z2  =  ^Ca  +  Cb^  {a67  +  *6  "  12  hb  ^a47  + 


+  J2  hbCb  (a37  +  a47  +  tt3  +  *4) 


Z3  =  (Ca  +  Cb^  {a57  +  *5  “  1*2  ha  ^a47  + 


+  12  haCb  (a37  +  a47  +  *3  +  Tr4) 


Z4  =  a31Zl  '  (allD  +  a15^  Z2 


Z5  =  a23°  Z3  “  ^a31°2  +  a32^  Z2 


(C.l  ) 


C  C.h  h 

aU  =Da  <Ca  +  V  ♦  -S-JJ- 


.  .  Wa  Z2_  V 

a12--1  VRb’ 


CAM 

a13  ■  °b<Ca  +  Cb>  +  -X1 
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14  "  TT  12 

O 


W  ,va  V 
ai5  ‘  ~r-(Ra  -r„’ 


/  CaWb\  -/VM 

'16  'U,  RJ  KRbJ 


a21  =  °b  ^Ca  +  Gb^  +  24  CaCbhb 


a22  ’  -  <Ca  +  Cb>  I  Gchb 


a23  =  TS  CaCbhahb 


-  L.  r  C  h  +  (C  +  CK)  |  GDh. 

12  C_CA  MC  Ru  a  b  6  c  b 


5  A 


24  12  abb  "a  "b 


a31  '  °a  (Ca  +  Cb>  +  If  CaCbhaZ 


a32  *  -  <Ca  +  Cb>  I  G>. 


a33  =  a23 


a34  =  T7  CaCbha(R! 


^>-K\  <Ca  +  Cb> 


(ci 
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J 


bll  =  a13a23  "  alla21 

b12  =  a14a23  ’  alla22  “  a 21 a 1 2 

b13  =  "  a12a22 

b14  =  al 5a 23  "  alla24 

b15  =  a16a23  “  a12a24 

b21  =  a23a23  “  a31 a21 

b22  =  “  ^a31a22  +  a32a21  ^ 

b23  =  “  a32a22 

b24  =  a23a34  "  a31 a 24 

b25  =  '  a32a24 
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g,  =  b  b  -  b  b 
1  21  14  11  24 


g?  =  b  b  +  b  b  -bb-bb 
L  21  15  22  14  12  24  11  25 


go  =  b  b  +  b  b  -  b  b  -bb 
J  22  15  23  14  12  25  13  24 


g.  =  b  b  -  b  b  (C.4) 

4  23  16  13  25 


L T ( x )  =  (b  D4  +  b  D2  +  b  )  Z  -  (b  D4  +  b  D2  +  b  )  Z 
1  21  22  23  4  11  12  13  5 


1 1  (x )  *  24 


LIIl(x)  =  Z2 


LjV(x)  =  V(Ca  +  Cb)  (a37  +a47  +tt3  +k4) 


Ly(x)  =  a^7  +tt4 


LVI(x)  =  a27  +  n2  (C .5) 


kl  '  Caha/2(Ca  +  Cb> 
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k2  =  Cahb/2(Ca  +  Cb) 


k3  =  1/R  (C,v  +  Cbvb)/(Ca  +  Cb) 


k4  ■  Cb 


k5  ~  k8  "  WRb 


k6  -  Cbhb  /12 


k7  *  V2 


k9  =  Cb/Rb‘ 


{3 


m  =  l/3j3{~)  -  A2 
91  91 


] 


!  2  (%3- 


n  *  1/27  /  2  (-^)J  -  9  (-^)  (-i)  +  27  (-*) 


l  9> 


V  -9, 


91 


A  *  n/2  +  ,i2/4  +  m3/27  j 
B  =J-  n/2  -  n2/4  +  m3/27^ 
X1  =  A  +  B  -  1/3  (g2/g, ) 


1/3 

1/3 


(C.6) 


(C.7) 


(C.8) 


X2  =  +1/2  (A  +  B)  +  1/3  (ga/g1) 


(C.9) 
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X3  (a  -  B) 

r  *  (x2)2  +  (x3)2 


(C.10) 


(C.U) 


(C.12) 


(C.13) 
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APPENDIX  D 


Derivation  of  Integrated  Shear  Stress  Resultant 
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The  required  shear  stress-strain  relation  can  be  developed  by 

•  • 

weighted  integration  in  order  to  obtain  the  factor  5/6  that  is  gen¬ 
erally  accepted  for  isotropic  plates  and  shells.  The  method  follows 
that  in  reference  (15).  For  convenience,  the  procedure  for  cylindrical 
shells  is  reproduced  here. 


Expressions  for  nonnal  stress  distributions  with  z  can  be  ob¬ 
tained  by  replacing  the  strains  in  the  stress-strain  relations  (1) 
by  the  approximate  forms  (8)  and  neglecting  ( z/R ) ^  in  comparison  to 
unity  in  the  same  terms.  It  is  to  be  understood  that  the  following 
hold  for  each  lamina. 


o 


x 


+ 


/ 

Z8)  - 


EaT  +  vEw 
!  v  R(l-v2) 


a 


6= 


EaT  +  Ew 
’'v  R(l-v2) 


(D.l) 

(D.2) 


Expressions  for  the  stress  distributions  in  terms  of  stress 
resultants  and  couples  can  be  obtained  by  using  equations  (11)  in 
(D.l)  and  (D.2). 


a 


X 


NTx 

TT^ 


vE 

Rh(l-v2) 


vEw 

R(l-v2) 


wdz 


_ E__ 

R(l-v2) 


(D.3) 


f 
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.c  |c\j  jc  |c\j 


Nu  NTo_  F 

°°  "  V 


h 

2  +  1?  z 
wdz  + 

-h  hJ 


MT 

M  +  Jl'- 

0  I  -V 


Rh(l-v  ) 


wzdz 


EaT  ,  _  E _ w 

1'v  R(l-v?) 


(OJ) 


The  shear  stress  is  related  to  the  normal  stress  by  the  equi¬ 
librium  equations  (2).  If  the  first  of  these  is  integrated  with  respect 
to  z,  the  following  is  obtained. 


CTxz  "  T2i 


a 


9X  (R’x>  dZ 


(D.5) 


Using  (D.3)  in  (D.5) 


o..  _  -  X 


xz  2i 


1 

Rh 


‘'Hz.; 


/  1?  f  /  1 

RN  +  4  2  RM  fdz  +  i 

x  x\  K7-h 


n  dz 


(D.6) 


where 


n  = 


R  a 
T-v  ax 


N 


L 


Tx  ,  12z  M 

ir +  :t  mt 
h  x 


EaT 


vE 

tt+^tr 


1 

h 


-h 


wdz  + 

h  / 


,h 

f  jz 

wzdz 

■h 

7 


-  w 


Refering  to  the  integrated  equilibrium  equations  (10)  and  re¬ 
placing  the  normal  stress  resultants  and  couples  in  (D.5)  by  their 
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equivalents  in  terms  of  the  shear  stress  resultant. 


zi  ■  ft j  rhii  -  <’  - +  ^  (n-B) 

Performing  the  indicated  integration 


ftdz 


(H.9) 


From  (D.6)  and  equations  (9),  it  can  readily  be  shown  that 


ndz  =  0 


nzdz  =  0 


CD. 10) 


The  shear  stress  distribution  (D.9)  satisfies  shear  stress 
boundary  conditions  and  the  definition  of  the  shear  stress  resultant. 
To  prove  the  latter  it  is  necessary  to  make  use  of  (D.10)  and 
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/-h 

2 


/  2 


-h 

2 


f(y)  dydz  = 


h 

2 


■h 

2 


(il 

v2 


-  z)  f(z)  dz 


(D.ll) 


McDonough  (15)  points  out  that  the  effect  of  the  last  term  in 
(D.9)  is  to  modify  the  classical  quadratic  shear  stress  distribution 
but  not  the  maqnitude  of  the  shear  stress  resultant.  The  modification 
is  due  to  the  nonlinearity  of  the  normal  stress  distribution,  primarily 
its  distribution  with  x. 

Proceeding  with  the  weiqhted  integration  as  in  (15),  the  shear 
stress-strain  relation  of  the  set  (1)  is  multiplied  by  the  weiqhtinq 
function  D  ■  <  *0  and  then  integrated  through  the  thickness  of  the 
shell  to.yield: 

2 

Ihe  integral  on  the  left  hand  side  is  evaluated  usinq  (D.9)  with 
the  aid  of  (D.10)  and  (D.ll).  The  integral  on  the  riqht  hand  side 
is  evaluated  by  usinq  the  shear  strain  given  in  equation  (8).  Aftpr 
Integration,  rearrangement  and  simplification,  equation  (11A)  results: 
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(D.13) 
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A  theory  for  the  analysis  of  stresses  in  laminated  circular  cylindrical  shells 
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taken  by  the  solution  function  is  governed  by  geometric  considerations.  A  range  in 
which  the  various  solution  forms  occur  was  determined  numerically. 

As  a  sample  problem,  the  slow  cooling  of  pyrolytic  graphite  deposited  onto  a 
commercial  graphite  mandrel  was  considered.  Investigation  of  normal  and  shear  stress 
behavior  at  the  pyrolytic  graphite  -  mandrel  interface  showed  that  these  stresses 
decrease  in  magnitude  with  increasing  E/G  ratio  and  increasing  deposit  to  mandrel 
thickness  (h^/h.) ratio.  This  implies  that  a  thin  mandrel  and  a  material  weak  in 
shear  are  desirable  to  minimize  the  possibilities  of  flaking  and  delamination  of  the 
pyrolytic  graphite. 
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